Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,452,273 G→A 5.2% A457V (GCA→GTA)  paaZ ← oxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,452,2730GA 5.2% 254.3 / 2.5 77A457V (GCA→GTA) paaZoxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (41/32);  new base A (2/2);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

ATCGTCGGACTACCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGCCAGTTGCAGAGCATGTCGCT  >  NC_000913/1452129‑1452417
                                                                                                                                                |                                                                                                                                                
aTCGTCGGACTACCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCg                                                                                                                                               <  1:359998/148‑1 (MQ=255)
 tcgtcgGACTACCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCTGGTGGATTCTTTTGCCGa                                                                                                                                              >  7:128305/1‑148 (MQ=255)
           aCCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCAtt                                                                                                                                     >  1:125827/1‑147 (MQ=255)
           aCCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCAtt                                                                                                                                     >  5:343743/1‑147 (MQ=255)
           aCCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTACCGACTCTTCAtt                                                                                                                                     >  8:26985/1‑147 (MQ=255)
           aCCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTACCGACTCTTCAtt                                                                                                                                     >  8:20558/1‑147 (MQ=255)
            cccTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGa                                                                                                                                   <  5:32510/148‑1 (MQ=255)
             ccTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAg                                                                                                                                  >  2:332196/1‑148 (MQ=255)
              cTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAgg                                                                                                                                 >  8:315140/1‑148 (MQ=255)
                gAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGAt                                                                                                                               >  4:354560/1‑148 (MQ=255)
                     gCGGTTCGCTGCATGTAATGTTTCACCGCGCGTAAACCGCCTAATTCTTCACCGCCGCCTGCGCCCCCGGGCCCCCCATGTACCAGTTGTGGCAGTGGGCAGCCATGCCTGGGGGATAGTTTTGCGGACTCTTCATGGCGGAGCGGaa                                                                                                                          >  6:374076/1‑148 (MQ=255)
                      cGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAAt                                                                                                                         <  1:332196/148‑1 (MQ=255)
                         ttCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCg                                                                                                                      <  8:58575/148‑1 (MQ=255)
                          tCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATc                                                                                                                              <  4:81928/139‑1 (MQ=255)
                          tCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATc                                                                                                                              >  3:81928/1‑139 (MQ=255)
                           cGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGc                                                                                                                     <  6:74848/147‑1 (MQ=255)
                           cGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGc                                                                                                                     <  2:125827/147‑1 (MQ=255)
                               gCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTGATTGAGGATCTGAATTCGCCCATg                                                                                                                >  2:341683/1‑148 (MQ=255)
                                  tGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCg                                                                                                              >  2:45640/1‑147 (MQ=255)
                                  tGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCg                                                                                                              >  5:406096/1‑147 (MQ=255)
                                    tAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTACCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTa                                                                                                            <  7:26985/147‑1 (MQ=255)
                                    tAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTACCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTa                                                                                                            <  7:20558/147‑1 (MQ=255)
                                     aaTGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGAt                                                                                                                               <  2:102337/127‑1 (MQ=255)
                                     aaTGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGAt                                                                                                                               >  1:102337/1‑127 (MQ=255)
                                     aaTGTATCACCGCTCGTAAACCGCCTAATTCTCCACCGCCTCCTGCGCGACAACGCCCACCATGTACCAGTTGTGCCAGTCGGGAGCCAGGCCCGGTGGATTCTTTAGCCGACTCTTCATTGAGGATCGGAATTCGCGCAAGCGTACg                                                                                                          >  7:189141/1‑148 (MQ=255)
                                           tCACCGCTCGTAAACCGCCTATTTCTTCACCGCCTCATGCGCGACAAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCcg                                                                                                     <  6:343743/147‑1 (MQ=255)
                                            cACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCcgcg                                                                                                   >  2:153769/1‑148 (MQ=255)
                                             aCCGCTCGTAATCCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGGCTCGTCATTGAGGAGCTGAATTCGCCCATGCGGACGTGACGCGt                                                                                                  >  2:127424/1‑148 (MQ=255)
                                              ccGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTc                                                                                                 >  1:271335/1‑148 (MQ=255)
                                                gCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCgg                                                                                               >  5:83805/1‑148 (MQ=255)
                                                gCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCgg                                                                                               >  5:293437/1‑148 (MQ=255)
                                                 cTCGTAAACCGCCTTATTCTTCACGGCCTCATGCGCGACCTGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCcgc                                                                                                    <  6:28860/142‑1 (MQ=255)
                                                  acGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCcgcg                                                                                                   >  5:28860/2‑142 (MQ=255)
                                                          cGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCg                                                                                                                      >  3:188736/1‑115 (MQ=255)
                                                          cGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCg                                                                                                                      <  4:188736/115‑1 (MQ=255)
                                                          cGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTg                                                                                     <  1:153769/148‑1 (MQ=255)
                                                           gCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGa                                                                                    >  8:363491/1‑148 (MQ=255)
                                                             cTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCAAGTACCACTTGTGGCAGAGGGGACCCATGCCCGGTGGAGTCTTTTGCCGACTATTCATTGAGCATCTGAATTCGCCAATGCGTACGTGCCGCGTCGGCAATAAACTGCcg                                                                                  >  2:406293/1‑148 (MQ=255)
                                                              tAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATTCCCGGTGGATTCTTTTGc                                                                                                                                                 <  3:61601/84‑1 (MQ=255)
                                                               aaTTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGTGAGCCATGCCCCGTGGATTCTTTTgac                                                                                                                                                >  4:61601/1‑82 (MQ=255)
                                                               aaTTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATATACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGAcgcg                                                                                >  3:132264/1‑148 (MQ=255)
                                                                 ttcttcACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGTAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCa                                                                              <  8:104808/148‑1 (MQ=255)
                                                                       aCCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCTGTGGGGAGCCATGCCTGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTg                                                                         <  5:374076/147‑1 (MQ=255)
                                                                            ctcctGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGa                                                                     <  5:200241/146‑1 (MQ=255)
                                                                            ctcctGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGa                                                                     >  6:200241/1‑146 (MQ=255)
                                                                              cctGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAg                                                                 >  6:248758/1‑148 (MQ=255)
                                                                               ctGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGc                                                                <  6:293437/148‑1 (MQ=255)
                                                                                     gACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCAc                                                          >  2:4456/1‑148 (MQ=255)
                                                                                            cccACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCCTGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGtt                                                   >  7:104372/1‑148 (MQ=255)
                                                                                                 caTGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCg                                               >  2:228526/1‑147 (MQ=255)
                                                                                                   tGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCaa                                            <  3:170016/148‑1 (MQ=255)
                                                                                                     tACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAg                                           <  5:248758/147‑1 (MQ=255)
                                                                                                        cAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTa                                       >  5:53821/1‑148 (MQ=255)
                                                                                                            tgtgGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGt                                                    >  7:193377/1‑131 (MQ=255)
                                                                                                            tgtgGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGt                                                    <  8:193377/131‑1 (MQ=255)
                                                                                                                    tGGGGAGCCATGCCCGGTGGATTCTTTTGCCTACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGc                                                                        <  6:136668/103‑1 (MQ=255)
                                                                                                                    tGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGc                                                                        >  5:136668/1‑103 (MQ=255)
                                                                                                                     ggggAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGca                           >  8:260081/1‑147 (MQ=255)
                                                                                                                      gggAGCCATGCCTGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGcaca                         <  8:128305/148‑1 (MQ=255)
                                                                                                                      gggAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGcac                          >  5:196436/1‑147 (MQ=255)
                                                                                                                            cATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGCCAg                   >  7:366646/1‑148 (MQ=255)
                                                                                                                             aTGCGCGGTGGCTTCTTTTGCCGACTGTTCATTCTGGCTCTCATGTCGCCAAGGCGAACGTGCCGCGGCGGCCATAAACTGCCGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCGGCAAGGCTACCGCCGCCTGCACAAGCCAGt                  <  3:269578/148‑1 (MQ=255)
                                                                                                                              tGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTcc                                                 <  4:134642/116‑1 (MQ=255)
                                                                                                                              tGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTcc                                                 >  3:134642/1‑116 (MQ=255)
                                                                                                                              tGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGCCAGtt                 >  5:185429/1‑148 (MQ=255)
                                                                                                                                    gTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGCCAGTTGCagag           <  7:363491/148‑1 (MQ=255)
                                                                                                                                     tGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTcc                                                 >  8:78201/1‑109 (MQ=255)
                                                                                                                                     tGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTcc                                                 <  7:78201/109‑1 (MQ=255)
                                                                                                                                     tGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGcc                     >  3:41311/1‑137 (MQ=255)
                                                                                                                                     tGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGcc                     <  4:41311/137‑1 (MQ=255)
                                                                                                                                      ggATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGCCAGTTGCAGAGCa         >  5:414799/1‑148 (MQ=255)
                                                                                                                                           cTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCa                                             <  8:13684/107‑1 (MQ=255)
                                                                                                                                           cTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCa                                             >  7:13684/1‑107 (MQ=255)
                                                                                                                                             tttGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCATCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGc                                         <  3:1458/109‑1 (MQ=255)
                                                                                                                                             tttGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCATCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGc                                         >  4:1458/1‑109 (MQ=255)
                                                                                                                                              ttGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGCCAGTTGCAGAGCATGTCGc   <  1:45640/146‑1 (MQ=255)
                                                                                                                                              ttGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGCCAGTTGCAGAGCATGTCGCt  <  7:260081/147‑1 (MQ=255)
                                                                                                                                              ttGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGCCAGTTGCAGAGCATGTCGCt  <  6:83805/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
ATCGTCGGACTACCCTGAACAGCGGTTCGCTGCATGTAATGTTTCACCGCTCGTAAACCGCCTAATTCTTCACCGCCTCCTGCGCGACCAGGCCCACCATGTACCAGTTGTGGCAGTGGGGAGCCATGCCCGGTGGATTCTTTTGCCGACTCTTCATTGAGGATCTGAATTCGCCCATGCGTACGTGCCGCGTCGGCAATAAACTGACGCGCAATTTGCGGATCAGCCGTCACCAGCGTTCCCGCAAGGCTACCGCCGCCTGCACAAGCCAGTTGCAGAGCATGTCGCT  >  NC_000913/1452129‑1452417

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: