Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,079,615 C→T 5.9% intergenic (‑285/+29) cmtB ← / ← tktA putative mannitol‑specific enzyme IIA component of PTS/transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,079,6150CT 5.9% 416.1 / 9.6 139intergenic (‑285/+29)cmtB/tktAputative mannitol‑specific enzyme IIA component of PTS/transketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (73/56);  new base T (4/4);  total (79/60)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

Only 100 of 143 total aligned reads displayed.

GTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGG  >  NC_000913/3079466‑3079764
                                                                                                                                                     |                                                                                                                                                     
gTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGcc                                                                                                                                                       >  1:1555243/1‑150 (MQ=255)
  tataGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGa                                                                                                                                                     >  2:198074/1‑150 (MQ=255)
     aGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAATGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGc                                                                                                                                                  >  2:1515297/1‑150 (MQ=255)
      gATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAATCg                                                                                                                                                 >  1:97543/1‑150 (MQ=255)
        tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGAc                                                                                                                                               >  1:872660/1‑150 (MQ=255)
          taACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCt                                                                                                                                             <  2:1727463/150‑1 (MQ=255)
           aacaacAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCtt                                                                                                                                            <  2:1081486/150‑1 (MQ=255)
             caacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCtttt                                                                                                                                          <  1:1312981/150‑1 (MQ=255)
             caacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCtttt                                                                                                                                          <  1:45806/150‑1 (MQ=255)
              aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCttttt                                                                                                                                         <  1:793047/150‑1 (MQ=255)
               acaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCtttttt                                                                                                                                        <  1:762851/150‑2 (MQ=255)
               acaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCtttttt                                                                                                                                        <  2:1049017/150‑2 (MQ=255)
                caaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAACAGGTCGCCGAAGCGACCTTTTTAc                                                                                                                                       <  1:1008782/150‑1 (MQ=255)
                 aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTAcc                                                                                                                                      >  1:936453/1‑150 (MQ=255)
                   aTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCg                                                                                                                                    <  2:1347763/150‑1 (MQ=255)
                   aTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCg                                                                                                                                    <  2:1025688/150‑1 (MQ=255)
                    tGTCAATACGCATATCGTGATTCATATTACATAATTGATATTTAATTCAAATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGa                                                                                                                                   <  2:803090/150‑1 (MQ=255)
                       cAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAAt                                                                                                                                <  2:1682740/150‑1 (MQ=255)
                       cAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAAt                                                                                                                                >  2:946749/1‑150 (MQ=255)
                       cAATACACATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAAt                                                                                                                                <  1:1792621/150‑1 (MQ=255)
                          tACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCt                                                                                                                             <  2:1189176/150‑1 (MQ=255)
                            cGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTaa                                                                                                                           >  1:631649/1‑150 (MQ=255)
                            cGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTaa                                                                                                                           >  1:1054174/1‑150 (MQ=255)
                            cGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGAGACCTTTTTACCCGAAATGCTaa                                                                                                                           >  1:625181/1‑150 (MQ=255)
                                 atCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAca                                                                                                                      <  2:503165/150‑1 (MQ=255)
                                  tCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcag                                                                                                                     >  1:905765/1‑150 (MQ=255)
                                   cGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcagc                                                                                                                    <  1:1304328/150‑1 (MQ=255)
                                     tGATGCATATTACATAATTGATATTTATTTCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcagcag                                                                                                                  <  1:1177769/150‑1 (MQ=255)
                                        tGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGAGACCTTTTTACACGAAATGCTAATTACAGCAGtta                                                                                                               >  1:1682764/1‑149 (MQ=255)
                                               tACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCt                                                                                                        <  2:1325626/150‑1 (MQ=255)
                                                aCATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGAAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATTCTAATAACAGCAGTACTTTTGCtt                                                                                                       >  1:192896/1‑150 (MQ=255)
                                                 cATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCttt                                                                                                      <  1:948249/150‑1 (MQ=255)
                                                 cATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCttt                                                                                                      >  2:109660/1‑150 (MQ=255)
                                                   tAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCg                                                                                                    <  2:617429/150‑1 (MQ=255)
                                                   tAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCg                                                                                                    <  1:1198664/150‑1 (MQ=255)
                                                     aTTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCTca                                                                                                  >  1:1118912/1‑150 (MQ=255)
                                                      ttGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTATTTACAGCAGTTCTTTTGCTTTCGcaa                                                                                                 >  1:148749/1‑150 (MQ=255)
                                                           atttattGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACg                                                                                            >  1:562487/1‑150 (MQ=255)
                                                              tattGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTa                                                                                         >  2:1271940/1‑150 (MQ=255)
                                                               attGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAACGTCGCCGAAGCGACCTTTTTACACGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACCTTAt                                                                                        >  1:1570/1‑150 (MQ=255)
                                                                  gCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCaa                                                                                     >  1:263219/1‑150 (MQ=255)
                                                                  gCATATAAATATTTGTGATATACAACACGCCTTATCTATTGCTTTTCCTCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCTACAACGTTATCaa                                                                                     <  2:249172/150‑1 (MQ=255)
                                                                       taAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTg                                                                                <  1:1547373/150‑1 (MQ=255)
                                                                       taAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTg                                                                                >  2:1584774/1‑150 (MQ=255)
                                                                       taAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAAAGGACATATCAAGGTAATAAAAAAGGTCGCGGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCTACAGTg                                                                                >  1:859445/1‑150 (MQ=255)
                                                                          atatTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAg                                                                             <  1:823950/150‑1 (MQ=255)
                                                                           tatTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTCGCTTTCGCAACAACGTTATCAACAGTGAAGc                                                                            <  2:375797/150‑1 (MQ=255)
                                                                            atTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGcc                                                                           <  1:1751325/150‑1 (MQ=255)
                                                                             tttGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCAGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCg                                                                          >  2:1249966/1‑150 (MQ=255)
                                                                              ttGTGATCTACAACACGCCTTATCTATTGCTTGTCCACAACCGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGCCCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGa                                                                         >  2:300030/1‑150 (MQ=255)
                                                                                gtgATCTACAACACGCCTTATCTATTGCTTGTCCGCATACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTGCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAAc                                                                       <  2:1214009/150‑1 (MQ=255)
                                                                                gtgATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAAc                                                                       <  2:97543/150‑1 (MQ=255)
                                                                                  gATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAActc                                                                     >  1:723290/1‑150 (MQ=255)
                                                                                  gATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAActc                                                                     >  1:568904/1‑150 (MQ=255)
                                                                                   aTCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAActct                                                                    >  1:1255364/1‑150 (MQ=255)
                                                                                    tCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCtt                                                                   <  2:1572670/150‑1 (MQ=255)
                                                                                    tCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCtt                                                                   <  1:1216582/150‑1 (MQ=255)
                                                                                     cTACAACAAGCCTTATCTATTGCTTGTCCGCACACGGACATATCAAGGTACTAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTGTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTc                                                                  <  1:303218/150‑1 (MQ=255)
                                                                                       acaacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCaa                                                                <  1:337851/150‑1 (MQ=255)
                                                                                       acaacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCaa                                                                <  1:1129505/150‑1 (MQ=255)
                                                                                        caacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCaaa                                                               >  2:1174020/1‑150 (MQ=255)
                                                                                             cGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAAcagca                                                          >  2:1810740/1‑150 (MQ=255)
                                                                                              gCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAAcagcag                                                         >  1:1074414/1‑150 (MQ=255)
                                                                                               ccTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACagcagc                                                        <  1:1073716/150‑1 (MQ=255)
                                                                                                 ttATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGctc                                                      >  1:285990/1‑150 (MQ=255)
                                                                                                       aTTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGCCCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCgg                                                >  2:955453/1‑150 (MQ=255)
                                                                                                       aTTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCgg                                                >  1:450702/1‑150 (MQ=255)
                                                                                                          gCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGGCGAAGCTACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTGACGCAACAACGTTATCAACAGAGAAGCCGAACTCGTCAAACAGCAGCTCTGCCGGAGc                                             >  1:1389187/1‑150 (MQ=255)
                                                                                                           cTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAATCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCa                                            >  2:312965/1‑150 (MQ=255)
                                                                                                           cTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCa                                            <  1:359526/150‑1 (MQ=255)
                                                                                                           cTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCa                                            >  1:1478113/1‑150 (MQ=255)
                                                                                                              gTCCGCATACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGAt                                         >  1:1262772/1‑150 (MQ=255)
                                                                                                              gTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAATCTACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGAt                                         <  1:1379648/150‑1 (MQ=255)
                                                                                                              gTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGAt                                         >  1:1466430/1‑150 (MQ=255)
                                                                                                              gTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGAt                                         <  1:1524585/150‑1 (MQ=255)
                                                                                                              gTCAGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTGATTACCCGAAATGCTAATTTCCTCAGTTCTTTTGCTTTCTCAAACACGTTATCACCAGTGAAGCCGACCGCTACACACCTCAGCTCTTCCGGAGCAGAt                                         >  2:141521/1‑150 (MQ=255)
                                                                                                                 cGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCa                                      >  2:140253/1‑150 (MQ=255)
                                                                                                                 cGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCa                                      >  1:944154/1‑150 (MQ=255)
                                                                                                                   cAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCTCAACAACGTTATCAACAGTGAAGCCGAACTCATCAAACAGCAGCTCTGCCGGAGCAGATTCAcc                                    >  2:974163/1‑150 (MQ=255)
                                                                                                                       cGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAg                                >  2:548651/1‑150 (MQ=255)
                                                                                                                        ggACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCTGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAgg                               >  2:34432/1‑150 (MQ=255)
                                                                                                                           cATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggtgg                            <  2:905765/150‑1 (MQ=255)
                                                                                                                           cATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggtgg                            >  1:1730100/1‑150 (MQ=255)
                                                                                                                            atatCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggtggt                           <  2:18026/150‑1 (MQ=255)
                                                                                                                             taacAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATAACAGCAGTTCTTTTCCTTTCGCAACACCGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggtggac                          >  1:638659/4‑148 (MQ=255)
                                                                                                                               tCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCAt                        >  1:1701946/1‑150 (MQ=255)
                                                                                                                                       ataaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACcgacga                <  2:625181/150‑1 (MQ=255)
                                                                                                                                       ataaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACcgacga                <  2:631649/150‑1 (MQ=255)
                                                                                                                                       ataaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACAGAAGGTGGTCATACcgaaga                >  1:324426/1‑150 (MQ=255)
                                                                                                                                        taaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGAt               >  2:1111879/1‑150 (MQ=255)
                                                                                                                                        gtaaaaaGGTCGCTTCGGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGAt               >  2:1304676/3‑150 (MQ=255)
                                                                                                                                         aaaaaaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGGTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATa              <  1:304756/150‑1 (MQ=255)
                                                                                                                                          aaaaaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAg             >  2:1364574/1‑150 (MQ=255)
                                                                                                                                          aaaaaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAAAAACGTTATCAACAGTGAAGACGAACTATTCAAACAGCAGCTCTGCCGGAGCAGATTCACAGAAGGTGTTCATACCGACGATAg             >  1:189291/1‑150 (MQ=255)
                                                                                                                                                gTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGt       >  2:165620/1‑150 (MQ=255)
                                                                                                                                                  cGCCGAAGCGACCTTTTTAACCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTGCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTc     >  1:1003129/1‑150 (MQ=255)
                                                                                                                                                   gCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGCTAGCACCGTTCa    >  2:421017/1‑150 (MQ=255)
                                                                                                                                                   gCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCa    >  1:60825/1‑150 (MQ=255)
                                                                                                                                                   gCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCa    >  1:913663/1‑150 (MQ=255)
                                                                                                                                                     cGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAgg  <  1:1515297/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGG  >  NC_000913/3079466‑3079764

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: