Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,981,814 T→A 5.1% T192S (ACG→TCG)  otsB ← trehalose‑6‑phosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,981,8140TA 5.1% 546.1 / 9.9 168T192S (ACG→TCG) otsBtrehalose‑6‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (86/73);  new base A (5/4);  total (91/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

TTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATG  >  NC_000913/1981665‑1981963
                                                                                                                                                     |                                                                                                                                                     
ttCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGt                                                                                                                                                       <  1:1429280/150‑1 (MQ=255)
     gCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGcc                                                                                                                                                  >  2:852853/1‑150 (MQ=255)
     cccAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTGAAATCATCGCCCAGAAATACGGGCGTTCGcc                                                                                                                                                  >  1:2296924/2‑150 (MQ=255)
       cAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCg                                                                                                                                                <  1:693093/150‑1 (MQ=255)
           tCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATaa                                                                                                                                            >  1:476789/1‑150 (MQ=255)
           tCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGAACGCCCgtatt                                                                                                                                            <  1:839247/150‑5 (MQ=255)
            ccAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGAACGCCCgta                                                                                                                                              <  2:655652/147‑3 (MQ=255)
             cAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAg                                                                                                                                          >  2:157395/1‑150 (MQ=255)
              aGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAgg                                                                                                                                         <  1:1333498/150‑1 (MQ=255)
              aGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGAACGCCCgtatttct                                                                                                                                         <  1:2174894/150‑8 (MQ=255)
              aGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGAACGCCCgtatttct                                                                                                                                         <  2:1527987/150‑8 (MQ=255)
               gACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAggg                                                                                                                                        <  2:444962/150‑1 (MQ=255)
                aCATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGa                                                                                                                                       <  1:115907/150‑1 (MQ=255)
                 cATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCNGATAAAGGGAg                                                                                                                                      <  1:2407362/150‑1 (MQ=255)
                 cATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAg                                                                                                                                      >  1:1671187/1‑150 (MQ=255)
                  aTCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGc                                                                                                                                     <  1:1822000/150‑1 (MQ=255)
                   tCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCt                                                                                                                                    >  1:991074/1‑150 (MQ=255)
                     cGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTc                                                                                                                                  >  2:1422618/1‑150 (MQ=255)
                        cacaCCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTg                                                                                                                               <  2:1373893/150‑1 (MQ=255)
                         acacCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGc                                                                                                                              <  2:1383134/150‑1 (MQ=255)
                         acacCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGc                                                                                                                              <  2:476789/150‑1 (MQ=255)
                                ccAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAg                                                                                                                       <  2:620087/150‑1 (MQ=255)
                                 cAGTAGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGc                                                                                                                      <  1:1087825/150‑1 (MQ=255)
                                   gTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTg                                                                                                                    <  2:1358057/150‑1 (MQ=255)
                                     cGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCa                                                                                                                  >  1:1988487/1‑150 (MQ=255)
                                          tgatgCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGc                                                                                                             >  2:2102738/1‑150 (MQ=255)
                                          tgatgCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGc                                                                                                             >  2:308579/1‑150 (MQ=255)
                                           gatgCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGcc                                                                                                            <  2:1269967/150‑1 (MQ=255)
                                           gatgCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGcc                                                                                                            <  2:1853534/150‑1 (MQ=255)
                                               ccTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCAc                                                                                                        >  2:1130425/1‑150 (MQ=255)
                                               ccTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCAc                                                                                                        >  2:929060/1‑150 (MQ=255)
                                                     ttGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTAc                                                                                                  <  2:2198443/150‑1 (MQ=255)
                                                      tGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACt                                                                                                 >  2:1169597/1‑150 (MQ=255)
                                                         agcTGTGCAAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCTCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGt                                                                                              <  2:1705929/148‑1 (MQ=255)
                                                           cTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTAc                                                                                            <  2:2392113/150‑1 (MQ=255)
                                                           cTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTAc                                                                                            <  1:1316223/150‑1 (MQ=255)
                                                            tgtgCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTAcc                                                                                           <  1:887281/150‑1 (MQ=255)
                                                             gtgCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACct                                                                                          >  1:1206706/1‑150 (MQ=255)
                                                               gCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCTGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACctcg                                                                                        >  2:1967772/1‑149 (MQ=255)
                                                               gCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACctct                                                                                        <  2:2296924/150‑1 (MQ=255)
                                                                ccAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTAcc                                                                                           >  1:1465086/1‑146 (MQ=255)
                                                                ccAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTAcc                                                                                           <  2:1465086/146‑1 (MQ=255)
                                                                ccAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCtctc                                                                                       >  2:496538/1‑150 (MQ=255)
                                                                ccAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCtctc                                                                                       >  2:2136125/1‑150 (MQ=255)
                                                                  aaTTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCgg                                                                                     <  2:207536/150‑1 (MQ=255)
                                                                  aaTTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCgg                                                                                     <  2:1916695/150‑1 (MQ=255)
                                                                    tttttACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGtt                                                                                   >  1:1695170/1‑150 (MQ=255)
                                                                    tttttACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGtt                                                                                   >  2:1623148/1‑150 (MQ=255)
                                                                      tttACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTg                                                                                 >  1:2190835/1‑150 (MQ=255)
                                                                       ttACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGa                                                                                <  1:971971/150‑1 (MQ=255)
                                                                       ttACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGa                                                                                <  1:308579/150‑1 (MQ=255)
                                                                          cTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAtct                                                                             >  1:1539191/1‑150 (MQ=255)
                                                                          cTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAtct                                                                             >  2:2258757/1‑150 (MQ=255)
                                                                          cTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGATACCTCTCGGTTTGAtct                                                                             >  1:646702/1‑150 (MQ=255)
                                                                            gACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCg                                                                           >  2:879319/1‑150 (MQ=255)
                                                                              cATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGac                                                                         >  1:1357525/1‑150 (MQ=255)
                                                                               aTTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTGACCTTTACTGGTACCTCTCGGTTTGATCTCGaca                                                                        >  1:1529971/1‑150 (MQ=255)
                                                                               aTTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGaca                                                                        >  2:1072706/1‑150 (MQ=255)
                                                                                ttCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGacaa                                                                       >  1:541420/1‑150 (MQ=255)
                                                                                  cttcgcAGTCGTTAACCGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGacacca                                                                     >  2:199411/6‑150 (MQ=255)
                                                                                  cttcgcAGTCGTTAACCGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGacaaca                                                                     >  2:725736/6‑150 (MQ=255)
                                                                                  cttcgcAGTCGTTAACCGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGacaaca                                                                     >  1:748958/6‑150 (MQ=255)
                                                                                   ttcgcAGTCGTTAACCGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATTCGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAacac                                                                    >  2:1379916/5‑150 (MQ=255)
                                                                                   cGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAacac                                                                    <  2:1671187/150‑1 (MQ=255)
                                                                                      gcAGTCGTTAACCGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACt                                                                 >  1:661381/2‑150 (MQ=255)
                                                                                      ccAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACt                                                                 <  1:2094337/150‑1 (MQ=255)
                                                                                         gTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTc                                                              >  1:389632/1‑150 (MQ=255)
                                                                                            ggTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCct                                                           >  1:1143569/1‑150 (MQ=255)
                                                                                               tAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCctgct                                                        <  1:157395/150‑1 (MQ=255)
                                                                                               tAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCctgct                                                        <  1:1893157/150‑1 (MQ=255)
                                                                                                  cGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTa                                                     <  2:118020/150‑1 (MQ=255)
                                                                                                     cTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACg                                                  >  1:265458/1‑150 (MQ=255)
                                                                                                      tGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGc                                                 >  2:1991090/1‑150 (MQ=255)
                                                                                                       gCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGcc                                                >  2:114883/1‑150 (MQ=255)
                                                                                                           aGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCAttt                                            >  1:2060388/1‑150 (MQ=255)
                                                                                                            gCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTtg                                           <  1:1877172/150‑1 (MQ=255)
                                                                                                            gCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTtg                                           <  1:2059033/150‑1 (MQ=255)
                                                                                                               aGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTgg                                        <  2:2101337/150‑1 (MQ=255)
                                                                                                                   tCAGCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGCGCTTCCTGCATAAAAGCTGCACTTGCCTCACCTTTACTGGTACCTCTCGGGTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTCGCCAg                                    >  1:1418357/1‑150 (MQ=255)
                                                                                                                        ggTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTg                               >  1:1424182/1‑150 (MQ=255)
                                                                                                                        ggTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTg                               <  1:557115/150‑1 (MQ=255)
                                                                                                                        ggTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTg                               >  2:77234/1‑150 (MQ=255)
                                                                                                                         gTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGa                              >  1:819020/1‑150 (MQ=255)
                                                                                                                            aaaTCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCACTTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTa                           <  1:2312561/150‑1 (MQ=255)
                                                                                                                              atcatcGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAAt                         <  1:325790/150‑1 (MQ=255)
                                                                                                                                 atcGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACg                      >  2:908048/1‑150 (MQ=255)
                                                                                                                                    gCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTgc                  >  2:1069815/1‑150 (MQ=255)
                                                                                                                                     cccAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTgc                  <  1:1169597/150‑1 (MQ=255)
                                                                                                                                       cAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTgcgc                >  2:1849330/1‑150 (MQ=255)
                                                                                                                                       cAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTgcgc                >  2:1516496/1‑150 (MQ=255)
                                                                                                                                       cAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTgcgc                <  2:1357525/150‑1 (MQ=255)
                                                                                                                                        cctttatcGGGCGAACGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCt               >  2:2176807/8‑150 (MQ=255)
                                                                                                                                        cctttatcGGGCGAACGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCt               >  1:1494276/8‑150 (MQ=255)
                                                                                                                                        aGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCt               <  2:380151/150‑1 (MQ=255)
                                                                                                                                          aaaTACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTaa             >  1:1096813/1‑150 (MQ=255)
                                                                                                                                           aaTACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCAGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAAt            >  1:1997598/1‑150 (MQ=255)
                                                                                                                                            aTACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATg           <  1:1798373/150‑1 (MQ=255)
                                                                                                                                                gggCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCAt       >  2:2404264/1‑150 (MQ=255)
                                                                                                                                                gggCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCAt       <  1:2394695/150‑1 (MQ=255)
                                                                                                                                                     ttCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATg  <  2:1753065/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAATG  >  NC_000913/1981665‑1981963

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