Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,167,881 C→A 6.6% K122N (AAG→AAT mqsA ← antitoxin of the MqsRA toxin‑antitoxin system/DNA‑binding transcriptional repressor MqsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,167,8810CA 6.6% 367.1 / 8.4 111K122N (AAG→AATmqsAantitoxin of the MqsRA toxin‑antitoxin system/DNA‑binding transcriptional repressor MqsA
Reads supporting (aligned to +/- strand):  ref base C (56/47);  new base A (2/6);  total (58/53)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

GAGAAATTATAAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACTATAAA  >  NC_000913/3167733‑3168030
                                                                                                                                                    |                                                                                                                                                     
gagaAATTATAAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACTGATTTCATTCAATAGTTCTGGATGCt                                                                                                                                                      <  2:1547320/150‑1 (MQ=255)
gagaAATTATAAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCt                                                                                                                                                      <  1:67820/150‑1 (MQ=255)
  gaAATTATAAGCCCAACAACAATACGCCTGTGGCATTGTTGGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTa                                                                                                                                                    <  1:104623/150‑1 (MQ=255)
  gaAATTATAAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTa                                                                                                                                                    <  1:664498/150‑1 (MQ=255)
   aaaTTATAAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTGCATTAATTAACGGATTTCATTAAATAGTTATTGATGCTTATc                                                                                                                                                  >  2:174850/1‑150 (MQ=255)
         taAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCGGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGaa                                                                                                                                             >  2:1307670/1‑150 (MQ=255)
          aaGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAac                                                                                                                                            <  1:364518/150‑1 (MQ=255)
          aaGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAac                                                                                                                                            >  2:146506/1‑150 (MQ=255)
          aaGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATAAGCATCCAGAac                                                                                                                                            <  1:452310/150‑1 (MQ=255)
          aaGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATAAGCATCCAGAac                                                                                                                                            <  1:1439541/150‑1 (MQ=255)
           aGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATAAGCATCCAGAact                                                                                                                                           <  2:151062/150‑2 (MQ=255)
             cccAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATAAGCATCCAGAactat                                                                                                                                         >  2:1541873/1‑147 (MQ=255)
              ccAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATAAGCATCCAGAactatt                                                                                                                                        <  1:1541873/150‑5 (MQ=255)
               caacaacaaTACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTa                                                                                                                                       <  1:1471353/150‑1 (MQ=255)
               caacaacaaTACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATAAGCATCCAGAactattg                                                                                                                                       <  1:1056943/150‑6 (MQ=255)
               caacaacaaTACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGAAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATAAGCATCCAGAactattg                                                                                                                                       <  2:1198509/150‑6 (MQ=255)
                  caacaaTACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTaaaa                                                                                                                                    >  1:261222/1‑150 (MQ=255)
                  caacaaTACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTaaaa                                                                                                                                    <  1:1507508/150‑1 (MQ=255)
                    acaaTACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGt                                                                                                                                  <  1:594989/150‑1 (MQ=255)
                    acaaTACGCCTGTAGCATTGTTAGCTCAAAATTATAGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGt                                                                                                                                  >  2:1285500/1‑150 (MQ=255)
                     aaaTACGCCTGTGGAATTGTTCGCTCAAACTTATCGCCAGTGATTTGGCCCACACTCCGGTAAACAAAAAACCGAGTTGCCCCCGCTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGttt                                                                                                                                <  1:174850/149‑1 (MQ=255)
                       aTACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTg                                                                                                                               >  2:1447906/1‑150 (MQ=255)
                        tACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGa                                                                                                                              <  2:45462/150‑1 (MQ=255)
                          cGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGAtt                                                                                                                            >  1:124815/1‑150 (MQ=255)
                             cTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATAAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATtgtg                                                                                                                         <  2:199487/150‑1 (MQ=255)
                               gtgGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGa                                                                                                                       <  1:1159410/150‑1 (MQ=255)
                               gtgGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGa                                                                                                                       >  1:1294541/1‑150 (MQ=255)
                                  gCATTTTTCGCTCAAACTTATCCCGATTGATTTGGCTCACAATCCGGTAAATAAAAAAGCGATTTTCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAgg                                                                                                                    <  1:1077262/150‑1 (MQ=255)
                                  gCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAgg                                                                                                                    >  1:1422836/1‑150 (MQ=255)
                                     ttGTTCGCTCAAACTTATCGCTAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATg                                                                                                                 <  2:288182/150‑1 (MQ=255)
                                      tGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGa                                                                                                                <  2:64444/150‑1 (MQ=255)
                                       gTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAg                                                                                                               >  1:516180/1‑150 (MQ=255)
                                            cTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTgg                                                                                                          >  2:1458936/1‑150 (MQ=255)
                                            cTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTgg                                                                                                          <  1:1058422/150‑1 (MQ=255)
                                              cAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGc                                                                                                        <  2:516180/150‑1 (MQ=255)
                                                       cgcgAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCttt                                                                                               >  2:76252/1‑150 (MQ=255)
                                                        gcgAGTTATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCtttt                                                                                              <  2:575852/150‑1 (MQ=255)
                                                        gcgAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCtttt                                                                                              >  2:1458712/1‑150 (MQ=255)
                                                        gcgAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCtttt                                                                                              <  2:899156/150‑1 (MQ=255)
                                                          gAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTc                                                                                            >  2:142664/1‑150 (MQ=255)
                                                               aTTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAAc                                                                                       >  1:65474/1‑150 (MQ=255)
                                                                tttGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAAcg                                                                                      <  2:770192/150‑1 (MQ=255)
                                                                 ttGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAAcgc                                                                                     <  1:133037/150‑1 (MQ=255)
                                                                  tGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAAcgcg                                                                                    >  2:7730/1‑150 (MQ=255)
                                                                     cTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGaaa                                                                                 >  1:1126040/1‑150 (MQ=255)
                                                                     cTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGaaa                                                                                 <  2:1627347/150‑1 (MQ=255)
                                                                      tCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGaaaa                                                                                >  1:1256188/1‑150 (MQ=255)
                                                                          acTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCa                                                                            <  1:385948/150‑1 (MQ=255)
                                                                          acTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCa                                                                            <  1:1447906/150‑1 (MQ=255)
                                                                           cTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCAt                                                                           >  1:521864/1‑150 (MQ=255)
                                                                           cTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCAt                                                                           >  2:1501632/1‑150 (MQ=255)
                                                                              cGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTa                                                                        >  1:877802/1‑150 (MQ=255)
                                                                                  aaaGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACAcc                                                                    >  2:1521398/1‑150 (MQ=255)
                                                                                   aaGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCt                                                                   >  2:301614/1‑150 (MQ=255)
                                                                                    aGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGATTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACATc                                                                  >  2:1203071/1‑150 (MQ=255)
                                                                                       aaaaaGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTcccc                                                               <  2:727657/150‑1 (MQ=255)
                                                                                         aaaGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCa                                                             >  2:1005979/1‑150 (MQ=255)
                                                                                            gCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCaaaa                                                          <  1:1458712/150‑1 (MQ=255)
                                                                                                  tGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCg                                                    >  2:21144/1‑150 (MQ=255)
                                                                                                  tGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCg                                                    <  1:622727/150‑1 (MQ=255)
                                                                                                   gCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGc                                                   <  2:1126040/150‑1 (MQ=255)
                                                                                                   gCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGc                                                   >  2:1306254/1‑150 (MQ=255)
                                                                                                    cccccGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCt                                                  >  1:1588935/1‑150 (MQ=255)
                                                                                                      ctcGCTTTTCCATTAATTAACGGCTTTCATTCAAAAGTTCTGGCTGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTg                                                <  2:682775/148‑1 (MQ=255)
                                                                                                      cccGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCAAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTg                                                >  1:544357/1‑150 (MQ=255)
                                                                                                       ccGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGc                                               >  2:1159780/1‑150 (MQ=255)
                                                                                                              tCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCtttt                                        >  2:1022594/1‑150 (MQ=255)
                                                                                                                cATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTg                                      >  2:861637/1‑150 (MQ=255)
                                                                                                                cATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTg                                      >  1:1585568/1‑150 (MQ=255)
                                                                                                                cATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTg                                      <  2:2167/150‑1 (MQ=255)
                                                                                                                  ttaattaaCGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTggg                                    <  1:1239501/150‑1 (MQ=255)
                                                                                                                    aattaaCGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTa                                  <  2:1422836/150‑1 (MQ=255)
                                                                                                                    aattaaCGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTa                                  <  1:1623559/150‑1 (MQ=255)
                                                                                                                       taaCGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAaga                               <  1:632509/150‑1 (MQ=255)
                                                                                                                        aaCGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAagag                              <  1:464780/150‑1 (MQ=255)
                                                                                                                        aaCGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAagag                              <  1:975382/150‑1 (MQ=255)
                                                                                                                         aCGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAagaga                             >  2:394063/1‑150 (MQ=255)
                                                                                                                           ggATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAagagaga                           >  2:1251010/1‑150 (MQ=255)
                                                                                                                                tcattcaATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCttt                      >  1:88733/1‑150 (MQ=255)
                                                                                                                                tcattcaATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCttt                      <  1:1458936/150‑1 (MQ=255)
                                                                                                                                 cattcaATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCtttt                     >  2:1491125/1‑150 (MQ=255)
                                                                                                                                   ttcaATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCtttttt                   <  1:1134848/150‑1 (MQ=255)
                                                                                                                                    tctaTAGTTCTTGATGCTTATCGAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGGATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTc                  <  1:90208/147‑1 (MQ=255)
                                                                                                                                     caATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCg                 >  1:685387/1‑150 (MQ=255)
                                                                                                                                     caATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCg                 >  2:658395/1‑150 (MQ=255)
                                                                                                                                      aaTAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGa                >  1:1439897/1‑150 (MQ=255)
                                                                                                                                       aTAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGaa               >  2:488529/1‑150 (MQ=255)
                                                                                                                                          gTTCTGGATGCTTATCCAGACCACGAAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCAATTTTCGTAAAGCGAAAACGCATTTACACCTCCCCCAAAAATATCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTATTTCGAACCt            >  2:1279993/1‑150 (MQ=255)
                                                                                                                                          gTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCt            >  2:1095715/1‑150 (MQ=255)
                                                                                                                                            tCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTc          >  2:1169353/1‑150 (MQ=255)
                                                                                                                                             cTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCa         >  2:1193967/1‑150 (MQ=255)
                                                                                                                                              tGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCAc        >  1:1333515/1‑150 (MQ=255)
                                                                                                                                               ggATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACt       >  2:1008607/1‑150 (MQ=255)
                                                                                                                                               ggATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACt       >  1:724032/1‑150 (MQ=255)
                                                                                                                                                gATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACta      <  1:1159780/150‑1 (MQ=255)
                                                                                                                                                  tGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCg                 >  2:1230172/1‑137 (MQ=255)
                                                                                                                                                  tGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCg                 <  1:1230172/137‑1 (MQ=255)
                                                                                                                                                  tGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACtata    <  2:1585568/150‑1 (MQ=255)
                                                                                                                                                    cTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACTATaaa  >  2:1239304/1‑150 (MQ=255)
                                                                                                                                                    cTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACTATaaa  >  2:1613937/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAGAAATTATAAGCCCAACAACAATACGCCTGTGGCATTGTTCGCTCAAACTTATCGCGAGTGATTTGGCTCACACTCCGGTAAAGAAAAAAGCGAGTTGCCCCCGCTTTTCCATTAATTAACGGATTTCATTCAATAGTTCTGGATGCTTATCCAGAACACGTAAAAGTTTGATTGTGGAAGGATGAGGTTGGGCATTGCCTTTTTCGTAACGCGAAAACGCATTTACACCTCCCCCAAAAATTTCGCTTGCCTCTTTTTGGGTAAGAGAGAGCTTTTTTCGAACCTTCACTATAAA  >  NC_000913/3167733‑3168030

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: