Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,130,008 T→C 6.7% I58I (ATT→ATC metL → Bifunctional aspartokinase/homoserine dehydrogenase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,130,0080TC 6.7% 443.2 / 16.5 149I58I (ATT→ATCmetLBifunctional aspartokinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base T (74/65);  new base C (5/5);  total (79/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

Only 100 of 152 total aligned reads displayed.

GCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGC  >  NC_000913/4129859‑4130156
                                                                                                                                                     |                                                                                                                                                    
gCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCATTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGAACTCTCAGCCTGACGATATGATTGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGAtt                                                                                                                                                      <  2:935896/150‑1 (MQ=255)
gCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGAtt                                                                                                                                                      >  2:803470/1‑150 (MQ=255)
  gAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTaa                                                                                                                                                    >  2:1633923/1‑150 (MQ=255)
   aaaGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAAc                                                                                                                                                   >  2:54907/1‑150 (MQ=255)
      ggTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTgg                                                                                                                                                <  2:615891/150‑1 (MQ=255)
       gtcgtcAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGt                                                                                                                                               <  1:1854397/150‑1 (MQ=255)
        gcgtcAGCTGCATAACTTTGGTGGCAGTAGTCAGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTAAGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGtt                                                                                                                                              <  2:1378374/149‑1 (MQ=255)
         cgtcAGCTGCATTAATTTTGTTGCAGTAGTATTGCTGATGTGAAGTGTTATTTGAGTGTCGCGGGCATTATGGCGGATTACTCTCAGACTGACGATATGATGGTGGTATCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGtta                                                                                                                                             <  1:148260/150‑2 (MQ=255)
         cgtcAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTg                                                                                                                                             <  1:191591/150‑1 (MQ=255)
          gtcAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGa                                                                                                                                            >  1:191503/1‑150 (MQ=255)
          ggcAGCTGCATAAATTTGGTGGCAGGAGGCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGa                                                                                                                                            <  2:1460277/148‑1 (MQ=255)
           gcAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGttagt                                                                                                                                           <  2:1890613/149‑4 (MQ=255)
             aGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGGCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATTATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAAc                                                                                                                                         <  2:1148019/150‑1 (MQ=255)
             aGCTGCATAAATTTGGTAGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGttagtgg                                                                                                                                         <  1:939264/150‑6 (MQ=255)
             aGCCGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGttagtgg                                                                                                                                         <  2:15204/150‑6 (MQ=255)
              ggggCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGCCGATATGATGGTCGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAATCt                                                                                                                                        >  2:960262/4‑150 (MQ=255)
                tGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTaa                                                                                                                                      >  1:2189217/1‑150 (MQ=255)
                 gCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAg                                                                                                                                     >  1:1261142/1‑150 (MQ=255)
                  cATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGc                                                                                                                                    >  1:1087572/1‑150 (MQ=255)
                   aTAAATTTGTTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGcc                                                                                                                                   <  1:1472368/150‑1 (MQ=255)
                   aTAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGcc                                                                                                                                   <  1:8658/150‑1 (MQ=255)
                     aaaTTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAg                                                                                                                                 <  1:1408696/150‑1 (MQ=255)
                           ggtggCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGAt                                                                                                                           >  1:748327/1‑150 (MQ=255)
                             tggCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCg                                                                                                                         <  1:2021409/150‑1 (MQ=255)
                             gggCAGTAGGCTGGCTGATGTGAAGTGTTATTTGCGGGGCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCg                                                                                                                         <  2:173320/149‑1 (MQ=255)
                              ggCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGt                                                                                                                        >  2:251751/1‑150 (MQ=255)
                                 agtagtCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGtctc                                                                                                                     >  1:1577103/1‑150 (MQ=255)
                                      tCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTgc                                                                                                                <  1:1909462/150‑1 (MQ=255)
                                           cTGATGTGTATTGTTATTTGCGTGTCGCGGTCATTTTGGCGGAGTACTATCAGCCTGACGTTATGATGGTGGTTTCGGCCGCCGGTAGCACCACTTTCCAGTTGATTTACTGGTTGAAACTAATGCAGACCGATCTTCTCTCTGCGCATc                                                                                                           <  2:1251822/150‑1 (MQ=255)
                                            tGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCa                                                                                                          <  1:1940505/150‑1 (MQ=255)
                                                gtgAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGtt                                                                                                      <  1:1795469/150‑1 (MQ=255)
                                                 tgAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTc                                                                                                     <  2:2189217/150‑1 (MQ=255)
                                                  gAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTca                                                                                                    <  2:2133182/150‑1 (MQ=255)
                                                  gAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTca                                                                                                    >  2:1127841/1‑150 (MQ=255)
                                                  gAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTca                                                                                                    >  2:1375985/1‑150 (MQ=255)
                                                  gAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGCTTTCCGCCGCCGGTAGCACCACTAAACAGTTGATTATCTGGTTGAAACTAAGCCAGACCGCTCGTCTCTCTGCGCATCAGGTTca                                                                                                    >  2:222614/1‑150 (MQ=255)
                                                   aaGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTcaa                                                                                                   <  1:2183815/150‑1 (MQ=255)
                                                       gtTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACaaa                                                                                               >  1:1098951/1‑150 (MQ=255)
                                                         tATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACg                                                                                             >  2:50519/1‑150 (MQ=255)
                                                          aTTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGc                                                                                            >  2:36813/1‑150 (MQ=255)
                                                            ttGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTg                                                                                          <  2:936728/150‑1 (MQ=255)
                                                            ttGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTg                                                                                          <  1:1050605/150‑1 (MQ=255)
                                                            ttGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTg                                                                                          >  1:95750/1‑150 (MQ=255)
                                                             tGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACCAACGCTGc                                                                                         >  1:304562/1‑150 (MQ=255)
                                                              gCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCg                                                                                        <  1:504004/150‑1 (MQ=255)
                                                               cGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGGTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGAGCATCAGGTTCAACAAACGCTGCGt                                                                                       >  1:1499590/1‑150 (MQ=255)
                                                                 tgtCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCg                                                                                     <  2:2112713/150‑1 (MQ=255)
                                                                  gtCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGc                                                                                    <  1:1019679/150‑1 (MQ=255)
                                                                   tCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCt                                                                                   <  1:1860290/150‑1 (MQ=255)
                                                                    cgcgGGCATTATGGCGGAGTACTCTCAGCCTGACGATTTGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAACCTAAGCCAGCCCGATCGTCGCTCTGCGCATCAGGTACACCAAACGCTGCGTCGCTa                                                                                  >  2:129231/1‑150 (MQ=255)
                                                                    cgcgGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTa                                                                                  >  1:1494326/1‑150 (MQ=255)
                                                                        ggCATTATTGCGGAGTACGCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACAAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAg                                                                              <  1:775169/150‑1 (MQ=255)
                                                                        ggCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTACCCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTCCGTCGCTATCAg                                                                              >  1:975940/1‑150 (MQ=255)
                                                                         gCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGt                                                                             >  1:1027536/1‑150 (MQ=255)
                                                                         gCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGt                                                                             >  1:1649342/1‑150 (MQ=255)
                                                                          cATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTg                                                                            >  1:295336/1‑150 (MQ=255)
                                                                           aTTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGc                                                                           <  1:251751/150‑1 (MQ=255)
                                                                               tGGCGGAGTACTCTCAGCCTGACGATATGTTGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATc                                                                       <  2:1601952/150‑1 (MQ=255)
                                                                                 gCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTg                                                                     >  2:1305639/1‑150 (MQ=255)
                                                                                  cGGAGTACTCTCAGCCTGACGATATGATGGTGTTTTCCGCCGCCGGTAGCACCACTAACCAGTTGAGTAACTGGGTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGa                                                                    <  1:111704/150‑1 (MQ=255)
                                                                                   ggAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGAt                                                                   <  2:752154/150‑1 (MQ=255)
                                                                                     aGTACTCTCAGCCTGACGATATGATGGTGGTAGCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGGCGCTATCAGTGCGATCTGATTa                                                                 <  2:402075/150‑1 (MQ=255)
                                                                                       tACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGc                                                               <  2:1010908/150‑1 (MQ=255)
                                                                                        aCTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCg                                                              >  1:1131289/1‑150 (MQ=255)
                                                                                         ctctCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACGACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGTCCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCgg                                                             <  1:572953/150‑1 (MQ=255)
                                                                                         ctctCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCgg                                                             >  1:429687/1‑150 (MQ=255)
                                                                                           ctcAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAACCGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTc                                                           >  1:1050233/1‑150 (MQ=255)
                                                                                           ctcAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTc                                                           >  1:2136534/1‑150 (MQ=255)
                                                                                               gCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCt                                                       <  1:860674/150‑1 (MQ=255)
                                                                                                  tGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTAcc                                                    >  1:1691767/1‑150 (MQ=255)
                                                                                                  tGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTAcc                                                    >  1:685006/1‑150 (MQ=255)
                                                                                                   gACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTTAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTAccc                                                   <  1:1557387/150‑1 (MQ=255)
                                                                                                     cGATATGATGGTGGTTTCAGCCGGCGGAAGCACCTCTAAACAGTTGATTAACTGGTTGAAACTAAGCCATACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGc                                                 <  1:1375018/150‑1 (MQ=255)
                                                                                                          tgatgGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATGGTCTCTATGCGCATCAGGTTCAACAAACGCTGCGTCGCTTTCAGGGCGAGCTGATTAGAGGTCTGCTTCACGCTgaag                                            >  1:361577/1‑150 (MQ=255)
                                                                                                           gatgGTGGTTTCCGCCGCCGGTAGCACCAATAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTgaaga                                           <  1:167504/150‑1 (MQ=255)
                                                                                                               gtgggTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGcc                                       <  2:800849/150‑1 (MQ=255)
                                                                                                                 ggtTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGa                                     <  2:646592/150‑1 (MQ=255)
                                                                                                                   tttCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATa                                   >  2:1566433/1‑150 (MQ=255)
                                                                                                                       cgccgccgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGCCCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCt                               >  1:78321/1‑150 (MQ=255)
                                                                                                                        gccgccgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTc                              >  1:12020/1‑150 (MQ=255)
                                                                                                                          cgccgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCAt                            >  2:1835600/1‑150 (MQ=255)
                                                                                                                           gccgGTAGCACAACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCACTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGACGCCGATAGCCTCAtt                           >  2:1925124/1‑150 (MQ=255)
                                                                                                                            ccgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTa                          >  1:1949332/1‑150 (MQ=255)
                                                                                                                             cgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGATCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCATCATTAg                         >  1:364913/1‑150 (MQ=255)
                                                                                                                              ggTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGCTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAgc                        <  1:1258134/150‑1 (MQ=255)
                                                                                                                                    accacTAACCAGTTGATTCACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCtttt                  <  2:1717532/150‑1 (MQ=255)
                                                                                                                                     gtttcaaCCAGTTAATCAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTg                 >  1:116868/6‑150 (MQ=255)
                                                                                                                                     ccacTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGCTCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTg                 >  1:281274/1‑150 (MQ=255)
                                                                                                                                      cacTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGt                >  2:274406/1‑150 (MQ=255)
                                                                                                                                       acTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTc               <  1:660354/150‑1 (MQ=255)
                                                                                                                                       acTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTc               >  2:722711/1‑150 (MQ=255)
                                                                                                                                         caaCCTGTTAATCAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGCAGCGGATGGCCTCTTTAGCGCTTTTGGCAg             >  2:117857/2‑150 (MQ=255)
                                                                                                                                          aaCCAGTTAATCAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGACTCATTAGCGCTTTTGTCAGc            >  1:1965286/1‑150 (MQ=255)
                                                                                                                                           aCCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTGCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCg           <  2:1050233/150‑1 (MQ=255)
                                                                                                                                             cAGTTGATTAACTGGTTGAAACTAAGCCAGACCGCTCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTGTTGTCAGCGAc         >  2:1070678/1‑150 (MQ=255)
                                                                                                                                             cAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGAc         >  1:450921/1‑150 (MQ=255)
                                                                                                                                               gTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCt       >  1:1232071/1‑150 (MQ=255)
                                                                                                                                                  gATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGa    >  2:1793077/1‑150 (MQ=255)
                                                                                                                                                   aTTATCTGGTTGAAACTAAGCCATACCTATCGTCTCTCTGCGCATCAGGGTCAACAAACGCTTCGTCGCTATCATTGCGATATGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAg   <  1:1114080/150‑1 (MQ=255)
                                                                                                                                                    ttAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAgc  >  1:2129936/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGC  >  NC_000913/4129859‑4130156

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: