Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,167,883 T→C 5.1% K122E (AAG→GAG)  mqsA ← antitoxin for MqsR toxin; transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,167,8830TC 5.1% 368.2 / 5.9 118K122E (AAG→GAG) mqsAantitoxin for MqsR toxin; transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (59/53);  new base C (2/4);  total (61/57)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

AGAAATTATAAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACTATAAATT  >  NC_000913/3167734‑3168032
                                                                                                                                                     |                                                                                                                                                     
agaAATTATAAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCtt                                                                                                                                                       >  1:1674863/1‑150 (MQ=255)
agaAATTATAAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCtt                                                                                                                                                       <  2:400364/150‑1 (MQ=255)
 taaaTTATAAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTa                                                                                                                                                      <  2:1078009/149‑1 (MQ=255)
    aTTATAAGCCCATCAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAACGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATcc                                                                                                                                                   <  2:1200099/150‑1 (MQ=255)
    aTTATAAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCCAACTTCTCGCGAGTGATTTGGCGCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATcc                                                                                                                                                   <  2:1327329/150‑1 (MQ=255)
    aTTATAAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATcc                                                                                                                                                   <  2:1179281/150‑1 (MQ=255)
     ttATAAGCCCAACAACAATACCCCTGTGGCATTGTTCGCTCAAACTTATCTCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCa                                                                                                                                                  <  2:1262313/150‑1 (MQ=255)
         aaGCCCAACTACACTACGCCTGTGGAATTGTTCGCTAAAACTTATCGCGATTGATTTGGCTAACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTGTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATAAGCATCCAGAac                                                                                                                                              <  2:319451/150‑1 (MQ=255)
         aaGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAac                                                                                                                                              <  2:886507/150‑1 (MQ=255)
         aaGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATAAGCATCCAGAac                                                                                                                                              <  2:35236/150‑1 (MQ=255)
          aGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTATCGGAATTCATTCAATAGTTCTGGATGCTTATCCAGAaca                                                                                                                                             >  1:288539/1‑150 (MQ=255)
          aGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAaca                                                                                                                                             >  2:419375/1‑150 (MQ=255)
          aGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATAAGCATCCAGAact                                                                                                                                             <  2:871598/150‑2 (MQ=255)
             ccAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGt                                                                                                                                          <  1:1717593/150‑1 (MQ=255)
              caacaacaaTACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAATAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTa                                                                                                                                         <  1:116496/150‑1 (MQ=255)
                acaacaaTACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTaaa                                                                                                                                       <  1:659543/150‑1 (MQ=255)
                 caacaaTACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTaaaa                                                                                                                                      >  1:13372/1‑150 (MQ=255)
                   acaaTACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGt                                                                                                                                    >  2:1401554/1‑150 (MQ=255)
                    caaTACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGtt                                                                                                                                   >  1:1104084/1‑150 (MQ=255)
                       tACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGa                                                                                                                                <  1:1839495/150‑1 (MQ=255)
                         cGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGAtt                                                                                                                              >  1:373198/1‑150 (MQ=255)
                         cGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGAtt                                                                                                                              >  2:1402686/1‑150 (MQ=255)
                           ccTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCTCACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTATTTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAAGACGTAAAATTTTGATtgt                                                                                                                            >  2:47058/1‑150 (MQ=255)
                             tgtgGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTgg                                                                                                                          <  2:1364954/150‑1 (MQ=255)
                             tgtgGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTgg                                                                                                                          <  2:1816342/150‑1 (MQ=255)
                               tgGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGaa                                                                                                                        <  2:834936/150‑1 (MQ=255)
                                 gCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCCCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACAAGTAAAAGTTTGATTGTGGAAgg                                                                                                                      >  1:196325/1‑150 (MQ=255)
                                    ttGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATg                                                                                                                   <  2:1742092/150‑1 (MQ=255)
                                      gTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAg                                                                                                                 >  2:1886799/1‑150 (MQ=255)
                                         cGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGtt                                                                                                              <  1:448224/150‑1 (MQ=255)
                                          gCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTg                                                                                                             >  2:1444848/1‑150 (MQ=255)
                                           cTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTgg                                                                                                            >  1:895470/1‑150 (MQ=255)
                                             cAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGc                                                                                                          >  1:1538230/1‑150 (MQ=255)
                                                 cTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTg                                                                                                      <  2:557733/150‑1 (MQ=255)
                                                 cTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTg                                                                                                      >  1:647227/1‑150 (MQ=255)
                                                      cgcgAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCttt                                                                                                 <  2:1893913/150‑1 (MQ=255)
                                                       gcgAGTGATTTGGCTCAGACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCtttt                                                                                                <  1:163658/150‑1 (MQ=255)
                                                       gcgAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCtttt                                                                                                >  1:47162/1‑150 (MQ=255)
                                                        cgAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCttttt                                                                                               <  1:563710/150‑1 (MQ=255)
                                                          aGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCg                                                                                             <  2:1503793/150‑1 (MQ=255)
                                                            tGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTa                                                                                           <  1:744957/150‑1 (MQ=255)
                                                                   gCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGaa                                                                                    >  2:1403906/1‑150 (MQ=255)
                                                                    cTCACACTCCGGTAATGAATAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCATTAGTTCTGGTTGGTTATCCTGTACAGGTAAATGTTTGATTGTGGTAGGATGAGGTTGGGCTTTGCCTTTTTCGTAACGCGaaa                                                                                   <  2:1691900/150‑1 (MQ=255)
                                                                    cTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTAGATTGTGGAAGTATGAGGTTGGGCATTGCCTTTTTCGTAACGCGata                                                                                   >  1:1235303/1‑148 (MQ=255)
                                                                     tCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGaaaa                                                                                  <  1:826563/150‑1 (MQ=255)
                                                                     ccacaCTACCGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGGTGGGCATTGCCTTTTTCGTAACGCGaaaa                                                                                  <  1:246764/149‑1 (MQ=255)
                                                                      cacaCTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAAc                                                                                 >  1:1645810/1‑150 (MQ=255)
                                                                        cacTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGc                                                                               >  1:1230570/1‑150 (MQ=255)
                                                                         acTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCa                                                                              <  1:1369449/150‑1 (MQ=255)
                                                                         acTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCa                                                                              <  1:1151161/150‑1 (MQ=255)
                                                                          cTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGATCACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCAt                                                                             >  1:63904/1‑150 (MQ=255)
                                                                          cTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCAt                                                                             <  1:690856/150‑1 (MQ=255)
                                                                          cTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCAt                                                                             >  1:1261410/1‑150 (MQ=255)
                                                                          cTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCAt                                                                             >  2:1069677/1‑150 (MQ=255)
                                                                              ggTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTac                                                                         <  1:230388/150‑1 (MQ=255)
                                                                              ggTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTac                                                                         <  2:745704/150‑1 (MQ=255)
                                                                                tAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATTAGGTTGGGCATTGACTTTTTCGTAACGCGAAAACGAATTTacac                                                                       >  1:512072/1‑150 (MQ=255)
                                                                                 aaaGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCAGGCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAACCGCATTTacaac                                                                      >  1:1863636/1‑148 (MQ=255)
                                                                                  aaGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCt                                                                     <  2:274122/150‑1 (MQ=255)
                                                                                  aaGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGAAGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGATAACGCATTTACACCt                                                                     >  1:1038488/1‑150 (MQ=255)
                                                                                   aGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTc                                                                    >  1:1268914/1‑150 (MQ=255)
                                                                                    gAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTcc                                                                   >  1:1603335/1‑150 (MQ=255)
                                                                                      aaaaaGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGTATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACCCCTcccc                                                                 >  1:917447/1‑150 (MQ=255)
                                                                                        aaaGCGAGTTGCCCCCGCTATTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCa                                                               >  1:1820883/1‑150 (MQ=255)
                                                                                           gCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCaaaa                                                            >  2:1895350/1‑150 (MQ=255)
                                                                                                 tGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCg                                                      <  1:1684158/150‑1 (MQ=255)
                                                                                                 tGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCg                                                      >  1:1821671/1‑150 (MQ=255)
                                                                                                 tGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCg                                                      <  2:1533727/150‑1 (MQ=255)
                                                                                                 tGCCCACGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAATGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCTTTTACACCTCCCCCAAAAATTTCg                                                      >  1:256429/1‑150 (MQ=255)
                                                                                                  gCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGc                                                     >  2:633311/1‑150 (MQ=255)
                                                                                                    ccccGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCtt                                                   <  2:199865/150‑1 (MQ=255)
                                                                                                       cGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGcc                                                <  1:1303570/150‑1 (MQ=255)
                                                                                                       cGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGcc                                                <  2:793968/150‑1 (MQ=255)
                                                                                                          ttttCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGATTGATTGTGGATGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGACAAGGCATTTACACCTCCCCCAAAAATTTCGCTTGCatct                                             >  2:491399/1‑150 (MQ=255)
                                                                                                              ccATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGTTAGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCATCttttt                                         >  1:899003/1‑150 (MQ=255)
                                                                                                              ccATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCttttt                                         >  1:734292/1‑150 (MQ=255)
                                                                                                                 ttaattaaCGGATTTGATTCAATAGTTATGGATGCTTATCCAGAACACGTAATAGTTTGATTGTGGAAGGATTAGGTTGGTCATTGCCTTTTTGGTCACGCGAAAACGCATTTACACCTCCCCCAAATATTTCGCTTGACTCTTTTTggg                                      <  1:1722957/150‑1 (MQ=255)
                                                                                                                  taattaaCGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGt                                     >  2:118641/1‑150 (MQ=255)
                                                                                                                         cGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAagagag                              >  1:300601/1‑150 (MQ=255)
                                                                                                                         cGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAagagag                              <  2:1471530/150‑1 (MQ=255)
                                                                                                                          ggATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACTTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCACAAAAATTTCGCTTGCCTCTTTTTGGGTAagagaga                             >  1:1767090/1‑150 (MQ=255)
                                                                                                                          ggATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAagagaga                             >  1:1021445/1‑150 (MQ=255)
                                                                                                                          ggATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAagagaga                             >  1:321347/1‑150 (MQ=255)
                                                                                                                            aTTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGc                           <  2:1261410/150‑1 (MQ=255)
                                                                                                                                cattcaATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGGtttt                       >  2:1685404/1‑150 (MQ=255)
                                                                                                                                cattcaATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCtttt                       >  2:1463755/1‑150 (MQ=255)
                                                                                                                                 attcaATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCttttt                      >  2:650760/1‑150 (MQ=255)
                                                                                                                                    caATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCg                   >  1:1543548/1‑150 (MQ=255)
                                                                                                                                     acgtgTTCTGGATAAGCATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCACAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGa                  >  1:374846/5‑150 (MQ=255)
                                                                                                                                      aTAGTTCTTGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGaa                 <  1:1444848/150‑1 (MQ=255)
                                                                                                                                        aGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGCTTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCTAAAACGGATTTACACCTCCCCCTAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAAcc               <  2:288539/150‑1 (MQ=255)
                                                                                                                                        aGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAAcc               >  2:1879082/1‑150 (MQ=255)
                                                                                                                                        aGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAAcc               <  2:1104084/150‑1 (MQ=255)
                                                                                                                                              ggATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACt         >  2:71898/1‑150 (MQ=255)
                                                                                                                                              ggATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACt         >  2:1057102/1‑150 (MQ=255)
                                                                                                                                              ggATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACt         <  2:1532298/150‑1 (MQ=255)
                                                                                                                                              ggATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCACCCAAAAATTTCGCTTGCCTCTATTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACt         >  1:1332001/1‑150 (MQ=255)
                                                                                                                                                   cTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGACAGAGCTTTTTTCGAACCTTCACTATaaa    >  2:286985/1‑150 (MQ=255)
                                                                                                                                                    ttATCCAGAACAAGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGAAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACTATAAAt   <  1:450564/150‑1 (MQ=255)
                                                                                                                                                     tATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGATGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCACCCAAAAATTTCGCTTGCCTCTTTTAGCGTTAGAGAGAGCTTTTTTCGAACCTTCACTATAAAtt  >  1:749516/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGAAATTATAAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACTATAAATT  >  NC_000913/3167734‑3168032

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