Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,629,498 T→C 5.4% Q14R (CAA→CGA)  rbbA ← ribosome‑associated ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,629,4980TC 5.4% 291.3 / 4.8 93Q14R (CAA→CGA) rbbAribosome‑associated ATPase
Reads supporting (aligned to +/- strand):  ref base T (39/49);  new base C (3/2);  total (42/51)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGAG  >  NC_000913/3629350‑3629645
                                                                                                                                                    |                                                                                                                                                   
gACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTGACGGCCGCCAGTTg                                                                                                                                                    <  2:236513/150‑1 (MQ=255)
gACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTg                                                                                                                                                    <  2:482054/150‑1 (MQ=255)
 aCGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTgc                                                                                                                                                   >  2:770137/1‑150 (MQ=255)
  cgcgGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTgcg                                                                                                                                                  <  2:1208785/150‑1 (MQ=255)
      ggCACCGGAAATCAACTACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGAc                                                                                                                                              <  2:1064396/150‑1 (MQ=255)
      ggCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGAc                                                                                                                                              <  2:198236/150‑1 (MQ=255)
      ggCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGAc                                                                                                                                              >  1:458159/1‑150 (MQ=255)
       gCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACa                                                                                                                                             <  2:122973/150‑1 (MQ=255)
          ccGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTgcgc                                                                                                                                                 <  2:1320822/143‑1 (MQ=255)
          ccGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTgcgc                                                                                                                                                 >  1:1320822/1‑143 (MQ=255)
          ccGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGc                                                                                                                                          <  2:174586/150‑1 (MQ=255)
          ccGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGc                                                                                                                                          >  1:600928/1‑150 (MQ=255)
           cGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAACCCGACCATACGGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcg                                                                                                                                         >  2:684838/1‑150 (MQ=255)
                aTCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacg                                                                                                                                    >  2:889919/1‑150 (MQ=255)
                 tCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCTGTCGCGCAACTGGcgggcgtg                                                                                                                                   <  1:986909/150‑5 (MQ=255)
                 tCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                                                   <  1:691399/150‑1 (MQ=255)
                  cAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacggg                                                                                                                                  >  2:531352/1‑150 (MQ=255)
                   aaCGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacggga                                                                                                                                 >  1:904585/1‑150 (MQ=255)
                     cGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAAc                                                                                                                               >  2:209843/1‑150 (MQ=255)
                          acaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGtt                                                                                                                          >  2:911407/1‑150 (MQ=255)
                          acaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGtt                                                                                                                          >  1:1128867/1‑150 (MQ=255)
                          acaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGtt                                                                                                                          >  1:1241650/1‑150 (MQ=255)
                           caaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTc                                                                                                                         <  2:931651/150‑1 (MQ=255)
                           caaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTc                                                                                                                         <  1:89373/150‑1 (MQ=255)
                             aGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCa                                                                                                                       <  1:531352/150‑1 (MQ=255)
                               cTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGa                                                                                                                     >  1:421956/1‑150 (MQ=255)
                                 cGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATg                                                                                                                   >  1:1026098/1‑150 (MQ=255)
                                   aCTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCg                                                                                                                 <  2:84705/150‑1 (MQ=255)
                                   aCTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCg                                                                                                                 <  2:1387557/150‑1 (MQ=255)
                                    cTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGt                                                                                                                >  2:794745/1‑150 (MQ=255)
                                    cTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGt                                                                                                                <  2:904585/150‑1 (MQ=255)
                                      tCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCa                                                                                                              <  1:1016181/150‑1 (MQ=255)
                                          gACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGc                                                                                                          <  1:884788/150‑1 (MQ=255)
                                          gACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGc                                                                                                          <  1:943867/150‑1 (MQ=255)
                                            cGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCgg                                                                                                        <  1:976239/150‑1 (MQ=255)
                                                gTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAAc                                                                                                    >  2:449151/1‑150 (MQ=255)
                                                gTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAAc                                                                                                    >  2:590323/1‑150 (MQ=255)
                                                 tCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAAcc                                                                                                   <  2:944425/150‑1 (MQ=255)
                                                       ccAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTcacca                                                                                             >  1:529355/1‑150 (MQ=255)
                                                         aaTCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCa                                                                                                                       >  1:1198199/1‑122 (MQ=255)
                                                         aaTCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCa                                                                                                                       <  2:1198199/122‑1 (MQ=255)
                                                            cAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAgg                                                                                        <  1:1241034/150‑1 (MQ=255)
                                                            cAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAgg                                                                                        <  1:911407/150‑1 (MQ=255)
                                                             aGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgt                                                                                       <  1:495502/150‑1 (MQ=255)
                                                              gCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtc                                                                                      >  2:1358164/1‑150 (MQ=255)
                                                                ccGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                                                                    >  1:575200/1‑150 (MQ=255)
                                                                     cATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGc                                                                               <  1:209843/150‑1 (MQ=255)
                                                                     cATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGc                                                                               >  1:987400/1‑150 (MQ=255)
                                                                          aGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACgg                                                                          <  2:271176/150‑1 (MQ=255)
                                                                                ccGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttc                                                                    >  1:1137207/1‑150 (MQ=255)
                                                                                  ggAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttctt                                                                  >  2:1278170/1‑150 (MQ=255)
                                                                                   gAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttcttc                                                                 <  2:1349322/150‑1 (MQ=255)
                                                                                   gAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttcttc                                                                 <  1:934159/150‑1 (MQ=255)
                                                                                     atatCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcat                                                               <  2:1008866/150‑1 (MQ=255)
                                                                                       atCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcattc                                                             <  2:1137207/150‑1 (MQ=255)
                                                                                        tCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcattca                                                            >  1:134202/1‑150 (MQ=255)
                                                                                         cGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCAc                                                                            <  1:1000317/133‑1 (MQ=255)
                                                                                         cGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCAc                                                                            >  2:1000317/1‑133 (MQ=255)
                                                                                            gagTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCg                                                        <  2:1241650/150‑1 (MQ=255)
                                                                                            gagTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCg                                                        <  1:1027221/150‑1 (MQ=255)
                                                                                            gagTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCg                                                        >  2:498387/1‑150 (MQ=255)
                                                                                              gTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCa                                                      <  2:1533453/150‑1 (MQ=255)
                                                                                                gATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAcc                                                    >  2:39317/1‑150 (MQ=255)
                                                                                                 atatTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAccc                                                   >  1:72272/1‑150 (MQ=255)
                                                                                                      gTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCt                                              >  2:328907/1‑150 (MQ=255)
                                                                                                       ttCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTa                                             <  1:859559/150‑1 (MQ=255)
                                                                                                       ttCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTa                                             <  1:1358164/150‑1 (MQ=255)
                                                                                                       ttCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTa                                             <  1:449151/150‑1 (MQ=255)
                                                                                                        tCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTAc                                            >  1:1094836/1‑150 (MQ=255)
                                                                                                         cAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCCGGCCACGGAAGTTCTTCATTCACCCGCACGCCCGCGAcc                                           >  2:985931/1‑149 (MQ=255)
                                                                                                               aaCGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGc                                     <  2:575200/150‑1 (MQ=255)
                                                                                                                aCGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCa                                    >  2:1395587/1‑150 (MQ=255)
                                                                                                                   gTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAAc                                 <  1:1397591/150‑1 (MQ=255)
                                                                                                                   gTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAAc                                 <  2:239942/150‑1 (MQ=255)
                                                                                                                       ttCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGt                             >  1:702849/1‑150 (MQ=255)
                                                                                                                          cATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt                          <  1:770137/150‑1 (MQ=255)
                                                                                                                            tAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGa                        <  2:600928/150‑1 (MQ=255)
                                                                                                                               tGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACat                     >  1:766899/1‑150 (MQ=255)
                                                                                                                                   ggCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTc                 <  1:684838/150‑1 (MQ=255)
                                                                                                                                    gCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTcc                <  2:529355/150‑1 (MQ=255)
                                                                                                                                    gCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTcc                <  1:889919/150‑1 (MQ=255)
                                                                                                                                     cTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCa               <  1:842931/150‑1 (MQ=255)
                                                                                                                                        aCGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAt            <  2:1433484/150‑1 (MQ=255)
                                                                                                                                         gtcccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCCCGAGGTCGTCAGGCCACGGAAGTTCTTCATTCCCCCGCACCCACGCTACGCCCGGCAACCCGGTTTTGACATATTCCAGAtg           >  1:828511/3‑150 (MQ=255)
                                                                                                                                         gtcccgccTGTCGCGCAACTGGCGGGCGTGAGCCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgc          <  1:52089/148‑1 (MQ=38)
                                                                                                                                          tcccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGAGTCATTGCGGCAACCTCACCGCGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCAGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgc          >  2:1063208/2‑150 (MQ=255)
                                                                                                                                           cccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgct         >  2:1187018/1‑150 (MQ=255)
                                                                                                                                           cccgccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgct         >  1:488857/1‑150 (MQ=255)
                                                                                                                                           cccgccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgct         >  1:1267777/1‑150 (MQ=255)
                                                                                                                                           cccgccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgct         >  2:1313809/1‑150 (MQ=255)
                                                                                                                                               ccAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGctgctg     <  2:1128867/150‑1 (MQ=255)
                                                                                                                                                cAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTgg    <  2:72272/150‑1 (MQ=255)
                                                                                                                                                  gTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGCAg  >  2:1417802/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGAG  >  NC_000913/3629350‑3629645

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: