Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,079,617 A→C 5.7% intergenic (‑287/+27) 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,6170AC 5.7% 416.8 / 13.0 140intergenic (‑287/+27)cmtB/tktAputative mannitol‑specific enzyme IIA component of PTS/transketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base A (75/57);  new base C (4/4);  total (79/61)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

Only 100 of 141 total aligned reads displayed.

TATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGG  >  NC_000913/3079468‑3079764
                                                                                                                                                     |                                                                                                                                                   
tataGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGa                                                                                                                                                     >  2:198074/1‑150 (MQ=255)
  taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAg                                                                                                                                                   <  1:816684/150‑1 (MQ=255)
  taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAg                                                                                                                                                   <  1:522311/150‑1 (MQ=255)
      tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGAc                                                                                                                                               >  1:872660/1‑150 (MQ=255)
      tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGAc                                                                                                                                               <  1:1469301/150‑1 (MQ=255)
         aacaacAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCtt                                                                                                                                            <  2:1081486/150‑1 (MQ=255)
          acaacaAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCttt                                                                                                                                           <  2:1345166/150‑1 (MQ=255)
           caacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCtttt                                                                                                                                          <  1:45806/150‑1 (MQ=255)
            aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCttttt                                                                                                                                         <  1:793047/150‑1 (MQ=255)
              caaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAACAGGTCGCCGAAGCGACCTTTTTAc                                                                                                                                       <  1:1008782/150‑1 (MQ=255)
                aaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTAccc                                                                                                                                     <  1:754358/150‑1 (MQ=255)
                 aTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCg                                                                                                                                    <  2:1025688/150‑1 (MQ=255)
                 aTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCg                                                                                                                                    <  2:1347763/150‑1 (MQ=255)
                  tGTCAATACGCATATCGTGATTCATATTACATAATTGATATTTAATTCAAATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGa                                                                                                                                   <  2:803090/150‑1 (MQ=255)
                     cAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAAt                                                                                                                                >  2:946749/1‑150 (MQ=255)
                     cAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAAt                                                                                                                                <  2:1682740/150‑1 (MQ=255)
                     cAATACACATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAAt                                                                                                                                <  1:1792621/150‑1 (MQ=255)
                       aTACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGc                                                                                                                              >  2:1027611/1‑150 (MQ=255)
                         aCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTa                                                                                                                            <  1:1535647/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)
                             atatCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTa                                                                                                                        >  2:1538509/1‑150 (MQ=255)
                               atCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAca                                                                                                                      <  2:503165/150‑1 (MQ=255)
                                 cGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcagc                                                                                                                    <  1:1304328/150‑1 (MQ=255)
                                   tGATGCATATTACATAATTGATATTTATTTCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcagcag                                                                                                                  <  1:1177769/150‑1 (MQ=255)
                                      tGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTc                                                                                                               <  1:932267/150‑1 (MQ=255)
                                      tGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGAGACCTTTTTACACGAAATGCTAATTACAGCAGtta                                                                                                               >  1:1682764/1‑149 (MQ=255)
                                              aCATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGAAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATTCTAATAACAGCAGTACTTTTGCtt                                                                                                       >  1:192896/1‑150 (MQ=255)
                                               cATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCttt                                                                                                      <  1:948249/150‑1 (MQ=255)
                                                 tAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCg                                                                                                    <  2:617429/150‑1 (MQ=255)
                                                 tAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCg                                                                                                    >  2:1214789/1‑150 (MQ=255)
                                                            tattGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTa                                                                                         >  2:1271940/1‑150 (MQ=255)
                                                            tattGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTACGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTa                                                                                         >  1:1249028/1‑150 (MQ=255)
                                                                gCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCaa                                                                                     >  1:263219/1‑150 (MQ=255)
                                                                 cATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAAc                                                                                    >  2:297921/1‑150 (MQ=255)
                                                                    ataAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGt                                                                                 >  1:946307/1‑150 (MQ=255)
                                                                     taAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTg                                                                                <  1:1547373/150‑1 (MQ=255)
                                                                     taAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAAAGGACATATCAAGGTAATAAAAAAGGTCGCGGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCTACAGTg                                                                                >  1:859445/1‑150 (MQ=255)
                                                                      aaaTATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGa                                                                               >  2:1411559/1‑150 (MQ=255)
                                                                         tatTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTCGCTTTCGCAACAACGTTATCAACAGTGAAGc                                                                            <  2:375797/150‑1 (MQ=255)
                                                                          atTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGcc                                                                           <  1:1751325/150‑1 (MQ=255)
                                                                            ttGTGATCTACAACACGCCTTATCTATTGCTTGTCCACAACCGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGCCCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGa                                                                         >  2:300030/1‑150 (MQ=255)
                                                                              gtgATCTACAACACGCCTTATCTATTGCTTGTCCGCATACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTGCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAAc                                                                       <  2:1214009/150‑1 (MQ=255)
                                                                              gtgATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAAc                                                                       <  2:97543/150‑1 (MQ=255)
                                                                               tgATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAAct                                                                      >  2:779693/1‑150 (MQ=255)
                                                                                gATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAActc                                                                     >  1:723290/1‑150 (MQ=255)
                                                                                gATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAActc                                                                     >  1:568904/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)
                                                                                    tACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCa                                                                 <  1:163331/150‑1 (MQ=255)
                                                                                     acaacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAATGGGCGCCGAAGCGAGCTTTTTACCCGAAATGCTAATTACAGCAGGTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCaa                                                                <  1:880009/150‑1 (MQ=255)
                                                                                     acaacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCaa                                                                <  1:1129505/150‑1 (MQ=255)
                                                                                     acaacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCaa                                                                <  1:337851/150‑1 (MQ=255)
                                                                                      caacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCaaa                                                               >  2:1174020/1‑150 (MQ=255)
                                                                                          acGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAAcagc                                                           <  2:196863/150‑1 (MQ=255)
                                                                                          acGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAAcagc                                                           >  1:1502808/1‑150 (MQ=255)
                                                                                           cGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGAGAAGCCGAACTCTTCAAAcagca                                                          >  1:318713/1‑150 (MQ=255)
                                                                                            gCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAAcagcag                                                         >  1:1074414/1‑150 (MQ=255)
                                                                                             ccTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACagcagc                                                        <  1:1073716/150‑1 (MQ=255)
                                                                                              cTTATCTATTGCTTGTCCGCAACCGGACATATCAAGGTAATAAAAACGGTCGCCGAAGCGAGCTTTTTACCCGAAATGCTACTTACAGCAGTTCTTTTGCTTTCGGAACAACGTTTTCAACAGTGACCCCGAACTCTTCAAACagcagca                                                       >  1:1150917/1‑149 (MQ=255)
                                                                                               ttATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGctc                                                      >  1:285990/1‑150 (MQ=255)
                                                                                                     aTTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGCCCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCgg                                                >  2:955453/1‑150 (MQ=255)
                                                                                                     aTTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCgg                                                >  1:450702/1‑150 (MQ=255)
                                                                                                         cTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAATCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCa                                            >  2:312965/1‑150 (MQ=255)
                                                                                                         cTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCa                                            >  1:1478113/1‑150 (MQ=255)
                                                                                                         cTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCa                                            <  1:359526/150‑1 (MQ=255)
                                                                                                           tGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGa                                          <  2:421263/150‑1 (MQ=255)
                                                                                                            gTCCGCATACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGAt                                         >  1:1262772/1‑150 (MQ=255)
                                                                                                            gTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAATCTACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGAt                                         <  1:1379648/150‑1 (MQ=255)
                                                                                                            gTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGAt                                         >  1:1466430/1‑150 (MQ=255)
                                                                                                            gTCAGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTGATTACCCGAAATGCTAATTTCCTCAGTTCTTTTGCTTTCTCAAACACGTTATCACCAGTGAAGCCGACCGCTACACACCTCAGCTCTTCCGGAGCAGAt                                         >  2:141521/1‑150 (MQ=255)
                                                                                                              ccGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTc                                       >  2:1145772/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)
                                                                                                                  aaaCGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCg                                   >  2:1622151/1‑150 (MQ=255)
                                                                                                                  aaaCGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCg                                   <  1:1279322/150‑1 (MQ=255)
                                                                                                                     cGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAg                                >  2:548651/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)
                                                                                                                          atatCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAAAGCAGTTCTTTTGATTTCGAAACAACGTTATAAAATGTGAAGAAGAACTCTTCAAACAGAAGATCAGCCGGAGCAGATACCACCAAggtggt                           >  2:34841/1‑150 (MQ=255)
                                                                                                                           taacAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATAACAGCAGTTCTTTTCCTTTCGCAACACCGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggtggac                          >  1:638659/4‑148 (MQ=255)
                                                                                                                             tCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCAt                        >  1:1701946/1‑150 (MQ=255)
                                                                                                                                     ggtaaaaaGGTCGCTTCGGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACcgacga                >  2:1273497/4‑150 (MQ=255)
                                                                                                                                     ataaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACcgacga                <  2:625181/150‑1 (MQ=255)
                                                                                                                                      taaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGAt               >  2:1111879/1‑150 (MQ=255)
                                                                                                                                      gtaaaaaGGTCGCTTCGGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCACCAACGTTTTCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATTCCGACGAt               >  1:1373579/3‑150 (MQ=255)
                                                                                                                                      gtaaaaaGGTCGCTTCGGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGAt               >  2:1304676/3‑150 (MQ=255)
                                                                                                                                       taaaaaGGTCGCTTCGGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATa              >  2:1089253/2‑150 (MQ=255)
                                                                                                                                       aaaaaaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGGTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATa              <  1:304756/150‑1 (MQ=255)
                                                                                                                                        aaaaaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAAAAACGTTATCAACAGTGAAGACGAACTATTCAAACAGCAGCTCTGCCGGAGCAGATTCACAGAAGGTGTTCATACCGACGATAg             >  1:189291/1‑150 (MQ=255)
                                                                                                                                              gTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGt       >  2:165620/1‑150 (MQ=255)
                                                                                                                                                cGCCGAAGCGACCTTTTTAACCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTGCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTc     >  1:1003129/1‑150 (MQ=255)
                                                                                                                                                 gCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCa    >  1:913663/1‑150 (MQ=255)
                                                                                                                                                 gCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCa    >  1:578442/1‑150 (MQ=255)
                                                                                                                                                   cGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAgg  >  2:1797862/1‑150 (MQ=255)
                                                                                                                                                   cGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAgg  <  1:1515297/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGG  >  NC_000913/3079468‑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: