Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,739,698 T→C 5.5% intergenic (+171/‑7) malS → / → avtA alpha‑amylase/valine‑pyruvate aminotransferase; transaminase C; alanine‑valine transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,739,6980TC 5.5% 2974.6 / 13.3 990intergenic (+171/‑7)malS/avtAalpha‑amylase/valine‑pyruvate aminotransferase; transaminase C; alanine‑valine transaminase
Reads supporting (aligned to +/- strand):  ref base T (532/370);  new base C (4/77);  total (536/454)
Fisher's exact test for biased strand distribution p-value = 1.89e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1036 total aligned reads displayed.

AGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGG  >  NC_000913/3739550‑3739846
                                                                                                                                                    |                                                                                                                                                    
    ccTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCg                                                                                                                                                  >  3:408257/1‑149 (MQ=255)
       cGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTa                                                                                                                                               <  4:285740/149‑1 (MQ=255)
        gggAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTAt                                                                                                                                              <  4:430152/149‑1 (MQ=255)
         ggAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTaag                                                                                                                                             >  1:2166277/1‑147 (MQ=255)
                  cAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTcc                                                                                                                                    <  2:2107607/149‑1 (MQ=255)
                      cAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCt                                                                                                                                                >  1:1019616/1‑133 (MQ=255)
                      cAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCttt                                                                                                                                >  3:1659875/1‑149 (MQ=255)
                      cAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCttt                                                                                                                                <  2:188772/149‑1 (MQ=255)
                       agtacaaaaTGTCAGCCCTACTCAGCGTTGCCGGATAAATCACCGCTTTCTCTTCAACAGTCCGATTTCCACCCCAGCAATTGTAGCGTTATTGTTACGTTCCCTGCTACAGATTTCGACGGATATCCCGCTATGACATTCTCCCtttt                                                                                                                               <  2:1359187/144‑1 (MQ=255)
                       aGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCtttt                                                                                                                               >  2:1741646/1‑149 (MQ=255)
                        ggATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCttttt                                                                                                                              <  1:1842282/149‑1 (MQ=255)
                        ggATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCttttt                                                                                                                              >  1:651521/1‑149 (MQ=255)
                        ggATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCttttt                                                                                                                              <  2:2308463/149‑1 (MQ=255)
                         gATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTg                                                                                                                             <  4:1038764/149‑1 (MQ=255)
                          aTAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTgg                                                                                                                            <  1:867886/149‑1 (MQ=255)
                           tAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGt                                                                                                                           <  1:1403950/149‑1 (MQ=255)
                           tAAAATGGCAGCCGTAAGCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCTACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACGTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGt                                                                                                                           <  4:2068264/149‑1 (MQ=255)
                              aaTGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGAc                                                                                                                        >  3:251968/1‑149 (MQ=255)
                                     gCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTa                                                                                                                 >  1:175319/1‑149 (MQ=255)
                                       ccTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTAcc                                                                                                               >  3:1392405/1‑149 (MQ=255)
                                           aTCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGcc                                                                                                           >  3:2175889/1‑149 (MQ=255)
                                              aGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACt                                                                                                        >  2:1830975/1‑149 (MQ=255)
                                                       ggATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTa                                                                                               >  4:1529929/1‑149 (MQ=255)
                                                          tAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGc                                                                                            >  3:2233304/1‑149 (MQ=255)
                                                            aaGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTg                                                                                          >  1:1647207/1‑149 (MQ=255)
                                                              gCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGAAAAATTTACCCGCCACTCCGGCATTACGCTGtt                                                                                        <  3:107481/149‑1 (MQ=255)
                                                               cACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATGTCCCGCTATGACATTCTGCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTg                                                                                       >  4:1458416/1‑149 (MQ=255)
                                                               cACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGt                                                                                                                           <  2:1226101/113‑1 (MQ=255)
                                                               cACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTg                                                                                       >  1:2400291/1‑149 (MQ=255)
                                                                 ccGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATACCCCGCTATGACATTCCCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGAt                                                                                     <  3:682921/149‑1 (MQ=255)
                                                                 ccGCTCACCATCCACCAGCCCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTAATTTCCTGGTTACGAAGTTCTACAAATACCCCGTTATCACTTCCTCTCTTTTTGGTCACAAATTTACCCGCCACTCCGGCATTACGCTGTTGAt                                                                                     <  4:871395/149‑1 (MQ=255)
                                                                       aCTCTTCAACGGACCGATTTGCACCCCAGAAAAGGTATCGTTATTGTTACCTTCCTGGCTACAGAGTCCGCGAATTTTCCCGCTAAAACATTCCCCCTTTTTGGTGAAAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGa                                                                               <  1:1988251/149‑1 (MQ=255)
                                                                        ctctCCAACAGACCGATTTTCACCCCAGAAAATGTAGCGTTATTGTTACCTCCCTTGCTCCAGATTTCGACAGACACCCCGCTATGACTTTCTCCTTTTTTGGTGACAATTTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGAT                                                                              <  2:2171767‑M1/149‑1 (MQ=255)
                                                                         tctTCACAACACCTATTTGCCCCCCAGAAAATATGCCGTTATTTTTACTTTCCTTGCAAGAAGATTTGACAGTTTTCCCGCAAGGACATCCTCCCTTTTTGGTAACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATc                                                                             <  4:773782/149‑1 (MQ=255)
                                                                         tctTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATAACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATC                                                                             <  3:2398470‑M1/149‑1 (MQ=255)
                                                                           ttCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGAt                                                                                     >  4:364899/1‑139 (MQ=255)
                                                                           ttCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTG                                                                           >  4:671290‑M1/1‑149 (MQ=255)
                                                                           ttCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTG                                                                           >  3:49849‑M1/1‑149 (MQ=255)
                                                                           ttCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTG                                                                           >  1:1114006‑M1/1‑149 (MQ=255)
                                                                            tCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGGTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCTGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGa                                                                          >  3:1945496/1‑149 (MQ=255)
                                                                            acAAAAGACCTATTGCCACCCCAGCAATTGTAGCTTTATTTTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCAAGAACTTTCTCCCTTTTTGGTCACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGa                                                                          <  1:1237520/148‑1 (MQ=255)
                                                                             cAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTCACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGAGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGaa                                                                         >  4:1821535/1‑149 (MQ=255)
                                                                                      gttttGCACCCCAGCAATTGTTGCGTTATTCTTATCTTCTTTGCTACGGAGTTCGACGGATGTCCCCCTATAACATCCTCCATTTTGGGTGCCAATTTTCCCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGttt                                                                <  4:2343649/147‑1 (MQ=255)
                                                                                      gATTTCCACCCCAAAAAATGTAGCGTTATTGTTACCTTCCTTCCTAAAGCTTCCGACAGATCCCCCGCTATAAAATTCTCCTTTTTTGGTGACAAATTTCCCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGttt                                                                <  3:766748/149‑1 (MQ=255)
                                                                                       aTTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTa                                                               >  1:1599194/1‑149 (MQ=255)
                                                                                       aTTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTa                                                               <  1:81654/149‑1 (MQ=255)
                                                                                        tttcccccccAGCAAATGTAGGGTTATTTTTACTTTCTTTGCAAAGGAGTTCGACAGATATCCCGCTTTAACATTCTCCCTTTTTGGTGAAAAATTTCCCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTac                                                              <  3:861215/143‑1 (MQ=255)
                                                                                        tttGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTac                                                              >  2:1006904/1‑149 (MQ=255)
                                                                                        tttGCACCCCAGAAAATGTAGCTTTATTTTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCTTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTac                                                              <  4:1981041/149‑1 (MQ=255)
                                                                                        tttGCACCCCAGAAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTac                                                              <  3:1227559/149‑1 (MQ=255)
                                                                                          tGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTacgc                                                            <  2:1538940/149‑1 (MQ=255)
                                                                                             aCCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACCAATTTACCCGCCACTCCGTCATTACGCTGACGATGTAAGATCTGAACGACGGTTTacgcacg                                                         >  2:841760/1‑149 (MQ=255)
                                                                                               cccAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGTACGACGGTTTACGCACGcc                                                       >  4:1673950/1‑149 (MQ=255)
                                                                                               cccAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGcc                                                       >  3:1021660/1‑149 (MQ=255)
                                                                                               cccAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGcc                                                       >  2:1109605/1‑149 (MQ=255)
                                                                                                  agaaaaTGTAGCGTTATTTTTACCTTCCTTGCTACAGAGTCCGACAGATATCCCGCTATGACTTTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTgg                                                    <  3:563902/146‑1 (MQ=255)
                                                                                                     aaaTGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACTTTCTCCCTTTTTGGTGAAAATTTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGgcgc                                                 <  2:924983/149‑1 (MQ=255)
                                                                                                     aaaTGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTATGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGgcgc                                                 >  3:2114121/1‑149 (MQ=255)
                                                                                                         gTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACTTTCTCCCTTTTTGGTGACAAATTTCCCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGAtt                                             <  2:3467/149‑1 (MQ=255)
                                                                                                          tAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGTTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTa                                            <  2:2400291/149‑1 (MQ=255)
                                                                                                           ggctttATTGTTCCCTTCCTTGCTACAGAGTTCGCAAGATATCCCGCTATGACATTCTCCCTTTTTGGTAAAAAATTTACCCGCCACACCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTAt                                           <  2:1666474/145‑1 (MQ=255)
                                                                                                            gCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATg                                          >  1:1537860/1‑149 (MQ=255)
                                                                                                              gTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCt                                        >  3:963212/1‑149 (MQ=255)
                                                                                                              gTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCt                                        >  2:591122/1‑149 (MQ=255)
                                                                                                               ttATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTc                                       >  4:1883777/1‑149 (MQ=255)
                                                                                                               ttATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTc                                       >  2:459562/1‑149 (MQ=255)
                                                                                                                  ttGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGaa                                                                         >  2:246745/1‑112 (MQ=255)
                                                                                                                    gTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTcggcgg                                  >  2:1149962/1‑149 (MQ=255)
                                                                                                                     ttACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTcggcggc                                 >  1:2278005/1‑149 (MQ=255)
                                                                                                                     ttACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTcggcggc                                 >  4:2224318/1‑149 (MQ=255)
                                                                                                                     ttACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTcggcggc                                 >  1:351231/1‑149 (MQ=255)
                                                                                                                     ttACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTcggcggc                                 >  4:1579774/1‑149 (MQ=255)
                                                                                                                      tACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTcggcggcg                                >  2:2232139/1‑149 (MQ=255)
                                                                                                                       aCCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTcggcggcgg                               >  4:245928/1‑149 (MQ=255)
                                                                                                                        ccttccttGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGt                              >  2:2034918/1‑149 (MQ=255)
                                                                                                                         cttccttGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGTCATTACGCTGTTGATGGAAGATCTGAACAACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTa                             >  2:1964128/1‑149 (MQ=255)
                                                                                                                            ccttCCTACAGAGTTCGAAAGATACCCCGCTATGACCTTCTCCTTTTTTGGTGACAAATTAACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATc                          <  4:2152221/149‑1 (MQ=255)
                                                                                                                             cttGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTcggcggcggagatcg                         >  3:1765617/1‑143 (MQ=255)
                                                                                                                              ttGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCg                        >  2:1426797/1‑149 (MQ=255)
                                                                                                                                   aCAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGaa                                                                                 >  3:2021652/1‑87 (MQ=255)
                                                                                                                                      gagTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATcc                         >  4:1172865/1‑140 (MQ=255)
                                                                                                                                        gTTCGACAGATATCCCGTTATGACTTTCTCCCTTTTTGTTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATcc              <  3:336349/149‑1 (MQ=255)
                                                                                                                                         ttCGACAGATATCCCGCTATGACCTTCTCACTTTTTGGTGACAAATTTACCCGCCACACCGGCATTACGCTGTTGATGGAAGATCTGAAAGACGGTTTACGCACGCATCGCGCGATTATGCTCGGCGGCGGTCATCCGGCGCAGATccc             >  4:1328788/1‑149 (MQ=255)
                                                                                                                                         ttCGACAAATACCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATccc             <  1:430042/149‑1 (MQ=255)
                                                                                                                                           cGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCgg           >  1:191479/1‑149 (MQ=255)
                                                                                                                                            gACAGATATCCTGCTATGACATTCTCCCTTTTTGGTGACAAATTTACACGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGa          >  1:728402/1‑149 (MQ=255)
                                                                                                                                            gACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGGACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGa          >  3:2315122/1‑149 (MQ=255)
                                                                                                                                            gACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGa          >  3:1861098/1‑149 (MQ=255)
                                                                                                                                            gACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGa          >  3:1563248/1‑149 (MQ=255)
                                                                                                                                            gACAGATATCCCGCTATGACATTCCCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGcgccggattaccgccgccgagcataatcgcgccaggcgtgcgtaaaccgtcgttcagatcttcca          <  3:1309460‑M1/149‑66 (MQ=255)
                                                                                                                                             aCAGATATCCCGCTATGACATTCCCCCTTTTTGGTGAAAAATTTCCCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGaa         <  1:207239/149‑1 (MQ=255)
                                                                                                                                              cAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGaaa        <  2:766254/149‑1 (MQ=255)
                                                                                                                                              cAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGaaa        >  1:1164217/1‑149 (MQ=255)
                                                                                                                                               aGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAAt       >  1:561524/1‑149 (MQ=255)
                                                                                                                                                 atatCCCGATATAACATTCTCCTTTTTGGGTAACAATTTTCCCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGc     <  1:256272/149‑1 (MQ=255)
                                                                                                                                                   atCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAg   >  4:1250021/1‑149 (MQ=255)
                                                                                                                                                   atCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAg   <  3:9003/149‑1 (MQ=255)
                                                                                                                                                   atCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAg   >  2:1633356/1‑149 (MQ=255)
                                                                                                                                                   atCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCT                                                                            <  1:501659‑M1/76‑1 (MQ=255)
                                                                                                                                                    tCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAgg  >  4:5385/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGAAGATCTGAACGACGGTTTACGCACGCCTGGCGCGATTATGCTCGGCGGCGGTAATCCGGCGCAGATCCCGGAAATGCAGG  >  NC_000913/3739550‑3739846

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: