Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,563,295 Δ1 bp 6.0% intergenic (‑219/+39) ddpX ← / ← dosP D‑alanyl‑D‑alanine dipeptidase/oxygen‑sensing c‑di‑GMP phosphodiesterase DosP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,563,2950G. 6.0% 219.5 / 7.5 67intergenic (‑219/+39)ddpX/dosPD‑alanyl‑D‑alanine dipeptidase/oxygen‑sensing c‑di‑GMP phosphodiesterase DosP
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base . (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

AAGCGCCAGAATGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCAT  >  NC_000913/1563146‑1563444
                                                                                                                                                      |                                                                                                                                                     
aaGCGCCAGAATGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACg                                                                                                                                                       <  2:658719/150‑1 (MQ=255)
 aGCGCCAGAATGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGt                                                                                                                                                      >  2:769762/1‑150 (MQ=255)
          aTGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAACGTGACACTGTCACTTAAAtactg                                                                                                                                             <  2:8979/150‑5 (MQ=255)
             cTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAAcc                                                                                                                                          <  1:1235248/150‑1 (MQ=255)
                gCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCcgc                                                                                                                                       >  1:771158/1‑150 (MQ=255)
                gCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCcgc                                                                                                                                       >  1:208841/1‑150 (MQ=255)
                 cTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCcgcg                                                                                                                                      <  2:1377876/150‑1 (MQ=255)
                   tAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAg                                                                                                                                    >  1:1054691/1‑150 (MQ=255)
                   tAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAg                                                                                                                                    >  2:924724/1‑150 (MQ=255)
                      aGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGc                                                                                                                                 <  2:600068/150‑1 (MQ=255)
                          atatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGc                                                                                                                             >  2:494872/1‑150 (MQ=255)
                          atatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGc                                                                                                                             >  1:578836/1‑150 (MQ=255)
                          atatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGc                                                                                                                             <  2:729103/150‑1 (MQ=255)
                          atatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGc                                                                                                                             <  1:431689/150‑1 (MQ=255)
                           tatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCg                                                                                                                            >  1:715562/1‑150 (MQ=255)
                           tatTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCg                                                                                                                            <  2:34067/150‑1 (MQ=255)
                                tttAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAgg                                                                                                                       <  1:460388/150‑1 (MQ=255)
                                    aCTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAAtt                                                                                                                   <  2:927622/150‑1 (MQ=255)
                                          gATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAg                                                                                                             <  2:1054691/150‑1 (MQ=255)
                                          gATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAg                                                                                                             <  1:920930/150‑1 (MQ=255)
                                          gATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAg                                                                                                             >  1:924514/1‑150 (MQ=255)
                                                 aTTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCa                                                                                                      <  1:627327/150‑1 (MQ=255)
                                                     gtgCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCgg                                                                                                  >  1:277634/1‑150 (MQ=255)
                                                         aGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAac                                                                                              <  2:512479/150‑1 (MQ=255)
                                                          gCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCTCGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAaca                                                                                             <  2:250420/150‑1 (MQ=255)
                                                          gCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAaca                                                                                             >  1:1417004/1‑150 (MQ=255)
                                                              aGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACgct                                                                                         >  2:1268185/1‑150 (MQ=255)
                                                                       ttGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCa                                                                                <  2:1199498/150‑1 (MQ=255)
                                                                          gTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGcc                                                                             >  2:1115966/1‑150 (MQ=255)
                                                                                gagGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAAt                                                                       <  1:676024/150‑1 (MQ=255)
                                                                                       aGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcg                                                                >  2:733600/1‑150 (MQ=255)
                                                                                        gTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcgg                                                               <  2:1266088/150‑1 (MQ=255)
                                                                                          tAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcggcg                                                             <  2:1417004/150‑1 (MQ=255)
                                                                                          tAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcggcg                                                             >  2:1037349/1‑150 (MQ=255)
                                                                                           aaggaaggAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcggcgg                                                            >  2:672413/1‑150 (MQ=255)
                                                                                             ggaaggaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGt                                                          <  1:769762/150‑1 (MQ=255)
                                                                                              gaaggaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcgg                                                               >  2:574414/1‑144 (MQ=255)
                                                                                              gaaggaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcgg                                                               <  1:574414/144‑1 (MQ=255)
                                                                                               aaggaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAg                                                        <  2:765768/150‑1 (MQ=255)
                                                                                                  gaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcggcgg                                                            <  2:1021980/143‑1 (MQ=255)
                                                                                                  gaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcggcgg                                                            >  1:1021980/1‑143 (MQ=255)
                                                                                                        cccGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGa                                               >  2:1030689/1‑150 (MQ=255)
                                                                                                         ccGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGaa                                              <  2:578836/150‑1 (MQ=255)
                                                                                                         ccGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGaa                                              >  2:1407353/1‑150 (MQ=255)
                                                                                                            gAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAg                                           <  1:1037349/150‑1 (MQ=255)
                                                                                                               tCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGaaa                                        <  1:57734/150‑1 (MQ=255)
                                                                                                                    gCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATcc                                   >  2:1210464/1‑150 (MQ=255)
                                                                                                                      tGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCt                                 >  1:737439/1‑150 (MQ=255)
                                                                                                                       taaaTCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGAATTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTg                                <  2:1308437/149‑1 (MQ=255)
                                                                                                                       gAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTg                                <  2:924514/150‑1 (MQ=255)
                                                                                                                       gAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTg                                <  1:924724/150‑1 (MQ=255)
                                                                                                                         aaTCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGaa                              >  2:659752/1‑150 (MQ=255)
                                                                                                                            cacaGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATaa                           >  2:511330/1‑150 (MQ=255)
                                                                                                                             acaGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATacc                          >  2:956996/1‑148 (MQ=255)
                                                                                                                             acaGAAATTAAGTTACAGTAACACTTTAAATGAAAAAACAACACTGAGGTCGACATCAAATTATAAAAATTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAAc                          >  2:651376/1‑150 (MQ=255)
                                                                                                                                  aTTTAACGTGACACTGTCACTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGaaa                                        <  2:1292793/131‑1 (MQ=255)
                                                                                                                                  aTTTAACGTGACACTGTCACTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGaaa                                        >  1:1292793/1‑131 (MQ=255)
                                                                                                                                  aTTTAACGTGACACTGTCACTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGAc                     >  1:1071262/1‑150 (MQ=255)
                                                                                                                                  aTTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGAc                     >  2:387722/1‑150 (MQ=255)
                                                                                                                                     tAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGt                  >  2:244115/1‑150 (MQ=255)
                                                                                                                                      aaGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTg                 >  2:1396659/1‑150 (MQ=255)
                                                                                                                                       aGTGACAGTGTCACGTTAAATGAAAACCCGCGAATGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTgg                >  2:192147/1‑150 (MQ=255)
                                                                                                                                         gTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGa               >  2:1009240/1‑150 (MQ=255)
                                                                                                                                         gTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGa               >  2:319800/1‑150 (MQ=255)
                                                                                                                                             cAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCtt           <  1:1115966/150‑1 (MQ=255)
                                                                                                                                                   cACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAg     >  2:351187/1‑150 (MQ=255)
                                                                                                                                                     cGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCa   <  1:494872/150‑1 (MQ=255)
                                                                                                                                                     cGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCa   >  1:1063123/1‑150 (MQ=255)
                                                                                                                                                      gTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCAt  <  2:715562/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AAGCGCCAGAATGCTGGCTTAAAGTTATATTTTTTAACTTTTGATCAACATTTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCAT  >  NC_000913/1563146‑1563444

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: