Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,663,482 A→C 7.1% K14T (AAG→ACG)  suhB → inositol monophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,663,4820AC 7.1% 378.2 / 14.6 119K14T (AAG→ACG) suhBinositol monophosphatase
Reads supporting (aligned to +/- strand):  ref base A (45/65);  new base C (4/5);  total (49/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

CATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCCACAGCACA  >  NC_000913/2663334‑2663628
                                                                                                                                                    |                                                                                                                                                  
cATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCCGCGCGCAAg                                                                                                                                                   >  1:1162880/1‑150 (MQ=255)
cATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAg                                                                                                                                                   >  2:1654382/1‑150 (MQ=255)
cATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAg                                                                                                                                                   >  2:1651099/1‑150 (MQ=255)
cATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAg                                                                                                                                                   >  2:801653/1‑150 (MQ=255)
cATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCACGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAAAATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAg                                                                                                                                                   >  2:2377104/1‑150 (MQ=255)
 aTCTAAAGATTCTTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAgg                                                                                                                                                  <  1:693586/150‑1 (MQ=255)
 aTCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAgg                                                                                                                                                  <  1:693699/150‑1 (MQ=255)
  tCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGc                                                                                                                                                 <  2:1915679/150‑1 (MQ=255)
    tAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                                                                               <  1:1601834/150‑1 (MQ=255)
     acaGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTGTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCGGCGCGCACGGCgat                                                                                                                                              <  2:2290900/148‑3 (MQ=255)
     aaaGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgat                                                                                                                                              <  2:1327222/150‑3 (MQ=255)
     aaaGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgat                                                                                                                                              <  2:2001989/150‑3 (MQ=255)
     aaaGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgat                                                                                                                                              <  2:1637194/150‑3 (MQ=255)
        gATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGTGAGTGACCGATGCATCCGATGCTGAACATCGCCGTGGGGGGGGCGCGCAAGGCGGGtaa                                                                                                                                           >  1:169645/1‑150 (MQ=255)
        gATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaa                                                                                                                                           <  1:1789440/150‑1 (MQ=255)
         attattCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCCAGGCGGGtact                                                                                                                                          >  1:1323063/1‑148 (MQ=255)
         attattCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaat                                                                                                                                          >  1:301365/1‑150 (MQ=255)
            attCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTGTACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaattta                                                                                                                                       <  1:2224923/150‑1 (MQ=255)
            attCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaattta                                                                                                                                       <  2:1776389/150‑1 (MQ=255)
            attCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaattta                                                                                                                                       <  2:1776388/150‑1 (MQ=255)
             ttCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCAGTTCTTTAACATCCAGTGAGAGATACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaa                                                                                                                                      <  1:1484913/150‑1 (MQ=255)
              tCACGCATCTTATCAGAAAACGCAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTGTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaat                                                                                                                                     <  1:1871839/150‑1 (MQ=255)
               cACGCATCTGATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTTCGCCGTTTTCCCTTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGAGCAGCGCGCAAGGCGGGtaatttaatt                                                                                                                                    <  2:2218701/150‑1 (MQ=255)
                aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTg                                                                                                                                   <  1:665971/150‑1 (MQ=255)
                aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTg                                                                                                                                   <  1:665959/150‑1 (MQ=255)
                aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTg                                                                                                                                   <  1:1153569/150‑1 (MQ=255)
                    aTCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaa                                                                                                                               <  1:1172339/150‑1 (MQ=255)
                    aTCTGATCAGAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCTCCGTTTTCCCGTTCCGTCACAGCCAGGGAGAGAGACCGATGCATCCGCTGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaa                                                                                                                               <  1:1172310/150‑1 (MQ=255)
                     tCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaaa                                                                                                                              <  1:2100721/150‑1 (MQ=255)
                     tCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCCCGCACGGCGGGTAATTTAATTTCCaaa                                                                                                                              >  1:829400/1‑150 (MQ=255)
                       ttATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaaaaa                                                                                                                            <  2:2345058/150‑1 (MQ=255)
                       ttATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaaaaa                                                                                                                            <  1:2278344/150‑1 (MQ=255)
                       ttATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaaaaa                                                                                                                            <  1:2278428/150‑1 (MQ=255)
                       ttATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaaaaa                                                                                                                            >  2:322014/1‑150 (MQ=255)
                        tATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAAc                                                                                                                           >  1:1790092/1‑150 (MQ=255)
                          tCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTa                                                                                                                         <  2:2268158/150‑1 (MQ=255)
                           cATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                                        <  2:2236146/150‑1 (MQ=255)
                           cATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                                        <  1:973722/150‑1 (MQ=255)
                            aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCCCCGCGCCAGGCGGGGACTGTAATTGCCAAAAACTATg                                                                                                                       >  1:2159851/1‑150 (MQ=255)
                            aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATg                                                                                                                       >  2:1025594/1‑150 (MQ=255)
                                aaCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAc                                                                                                                   <  1:38094/150‑1 (MQ=255)
                                 aCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAcc                                                                                                                  >  2:755587/1‑150 (MQ=255)
                                         aGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGGGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCCGCGCGCAAGGCGGGGAATTTAATTTCCAAAAAAAATGAAAACCCCGACGc                                                                                                          >  1:1771360/1‑150 (MQ=255)
                                             gCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTa                                                                                                      >  2:1862577/1‑150 (MQ=255)
                                             gCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTa                                                                                                      <  2:1082008/150‑1 (MQ=255)
                                             gCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCAGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTa                                                                                                      <  2:1008535/150‑1 (MQ=255)
                                             gCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCAGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTa                                                                                                      <  2:1763375/150‑1 (MQ=255)
                                                gATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGaa                                                                                                   <  2:347645/150‑1 (MQ=255)
                                                gATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGaa                                                                                                   <  2:1935276/150‑1 (MQ=255)
                                                      gctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACAACCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCCGCGCCCAAGGCGGGGAATTTAAATGCCAAAAAAAATGAAACCCCCGACCCGGTAGAAGCGAGc                                                                                             >  1:883551/1‑150 (MQ=255)
                                                            ataCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGaca                                                                                       >  2:2014467/1‑148 (MQ=255)
                                                            ataCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGaaa                                                                                       >  2:2014648/1‑150 (MQ=255)
                                                              aCTCTGCGCCGTTTTCCCGTTCTTTAACATCCCGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCCGCGCGCCCCCCGGGGACTTTAATTGCCCAAAACAATGAAACCCCCGACCCGGGAGAAGCGCGccccaaaag                                                                                     >  1:504479/1‑143 (MQ=255)
                                                                tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                                                                                   <  2:64229/150‑1 (MQ=255)
                                                                tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                                                                                   <  2:64227/150‑1 (MQ=255)
                                                                     gacGTTTACCCGTTCTTTAACATCAAGTGAGAGAGACCGATGCATCCGATCCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGACAAAAACTATGAAACCACGGAAGCTGTAGAAGCGAGCCAGAAAGGCAGTAAc                                                                              >  2:861101/3‑150 (MQ=255)
                                                                      ccGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCGAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACg                                                                             >  2:35544/1‑150 (MQ=255)
                                                                      ccGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAAATTAATTTCCAAACACTATGAAAACCCGGACCCTGTAGAAGCGAGCCAGAAAGGCACTAACg                                                                             >  1:2196511/1‑150 (MQ=255)
                                                                         ttttCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAtt                                                                          <  1:87212/150‑1 (MQ=255)
                                                                             cccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGt                                                                      >  1:2067533/1‑150 (MQ=255)
                                                                                  tCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCa                                                                 <  1:1789437/150‑1 (MQ=255)
                                                                                  tCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCa                                                                 <  1:844351/150‑1 (MQ=255)
                                                                                  tCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCa                                                                 <  1:844377/150‑1 (MQ=255)
                                                                                    ctgaaCATCCAGTGAGAGAGACCGATGCAGCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTCATTGCCAAAACCTATGACACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAAc                                                               <  2:1355039/147‑1 (MQ=255)
                                                                                              aGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAg                                                     <  1:683259/150‑1 (MQ=255)
                                                                                              aGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAg                                                     <  1:800721/150‑1 (MQ=255)
                                                                                              aGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAg                                                     <  1:136392/150‑1 (MQ=255)
                                                                                                tGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCt                                                   <  1:904060/150‑1 (MQ=255)
                                                                                                tGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCt                                                   <  1:21057/150‑1 (MQ=255)
                                                                                                tGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCt                                                   >  2:438079/1‑150 (MQ=255)
                                                                                                      agaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGaa                                             <  2:1242005/150‑1 (MQ=255)
                                                                                                      agaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGaa                                             <  1:3993/150‑1 (MQ=255)
                                                                                                      agaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGaa                                             <  1:1737789/150‑1 (MQ=255)
                                                                                                      agaCCGAGGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGaa                                             <  2:1885878/150‑1 (MQ=255)
                                                                                                          cGATGCATCCGATGCTGAACATCGCCGTGCGCGCCGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCGGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAATCTGCCGAAGCgg                                         >  1:2186230/1‑150 (MQ=255)
                                                                                                          cGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCgg                                         >  2:1543646/1‑150 (MQ=255)
                                                                                                           gATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGt                                        <  2:885508/150‑1 (MQ=255)
                                                                                                           gATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGt                                        <  2:885503/150‑1 (MQ=255)
                                                                                                                  ccGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATc                                 >  2:791512/1‑150 (MQ=255)
                                                                                                                    gATGCTGAACATCGCCGTGCGCGCCGCGCGCAAGGCGGGTAAATTAATTTCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCTTGACCAACGTACATAAATCAGCCGAAACGGTGATTATCGa                               >  1:1507197/1‑150 (MQ=255)
                                                                                                                    gATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAACAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAACGGCAGTAACGATTTCGTGACCACCGTAGATAAAGCTGCCGAAGCGGTGATTATCGa                               <  2:1684282/150‑1 (MQ=255)
                                                                                                                    gATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGa                               <  1:1092572/150‑1 (MQ=255)
                                                                                                                      tGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGaca                             <  2:1008918/150‑1 (MQ=255)
                                                                                                                       gCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGacac                            <  2:1773000/150‑1 (MQ=255)
                                                                                                                       gCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGacac                            <  2:1633340/150‑1 (MQ=255)
                                                                                                                       gCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGacac                            >  2:2075873/1‑150 (MQ=255)
                                                                                                                         tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGa                          >  2:1937883/1‑150 (MQ=255)
                                                                                                                         tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGACAAAAAATATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAAGATTTAGTGACCAACCTAGATAAACAAGCCGAAGCGGTGATTATAGACACGa                          >  1:1565183/1‑150 (MQ=255)
                                                                                                                          gAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGAt                         >  2:775656/1‑150 (MQ=255)
                                                                                                                              aTCGCCGTGCGCGCAGCGCGCCAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGAAGCTGTAGAAGCGAGCCAGAAAGGCACTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGt                     >  1:929814/1‑150 (MQ=255)
                                                                                                                               ccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGaa                                             <  1:1993977/124‑1 (MQ=255)
                                                                                                                               ccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGaa                                             >  2:1993977/2‑125 (MQ=255)
                                                                                                                                 gCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAATCTGCCGAAGCGGTGATTATCGACACGCTTCGTaaa                  >  1:1483443/1‑150 (MQ=255)
                                                                                                                                   cGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGACAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATc                >  2:333254/1‑150 (MQ=255)
                                                                                                                                   cGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATc                >  2:642575/1‑150 (MQ=255)
                                                                                                                                          gcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCca         >  1:153135/1‑150 (MQ=255)
                                                                                                                                            agcgcgcaAGGCGTGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCcaca       <  1:1491190/150‑1 (MQ=255)
                                                                                                                                             gcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTCTGACACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCCACAg      <  2:1771360/150‑1 (MQ=255)
                                                                                                                                             gcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCCACAg      <  1:322341/150‑1 (MQ=255)
                                                                                                                                                 gcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCCACAGcaca  >  1:1138931/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACCCACAGCACA  >  NC_000913/2663334‑2663628

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑