Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,883,838 T→A 5.0% intergenic (‑109/‑498) dnaA ← / → rpmH chromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator/50S ribosomal subunit protein L34

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,883,8380TA 5.0% 359.6 / 6.8 119intergenic (‑109/‑498)dnaA/rpmHchromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator/50S ribosomal subunit protein L34
Reads supporting (aligned to +/- strand):  ref base T (44/69);  new base A (4/2);  total (48/71)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

Only 100 of 120 total aligned reads displayed.

CAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGC  >  NC_000913/3883691‑3883987
                                                                                                                                                   |                                                                                                                                                     
cAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGc                                                                                                                                                     >  2:353868/1‑150 (MQ=255)
 aaTCGGGCAAGACACTGCTGCCAAAGCTAAAGTGACACGGCGGACTCCACTCGCACAAAAGTCGATAATGACTAAGGCTGAACAATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCg                                                                                                                                                    <  2:960364/150‑1 (MQ=255)
  aTCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCgg                                                                                                                                                   >  1:788824/1‑150 (MQ=255)
  aTCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCAGGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGCCGCACGAACCGCTGTCTGCgg                                                                                                                                                   >  2:915985/1‑150 (MQ=255)
  aTCGGGCAAGACACTGCAGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAACAGGCGATAATGACTAAGGCTGAAACAGTCATGATTGTTGACGTCCGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCgg                                                                                                                                                   <  1:267305/150‑1 (MQ=255)
   tCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGt                                                                                                                                                  <  2:1059553/150‑1 (MQ=255)
   tCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGt                                                                                                                                                  <  1:995703/150‑1 (MQ=255)
     gggCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCCTGATTGTTGACGTACGTCGAAAAGCCCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTta                                                                                                                                                <  1:404447/150‑1 (MQ=255)
       gCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTtata                                                                                                                                              <  1:973269/150‑1 (MQ=255)
        cAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGGCTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCAGACAGCGGttcgtg                                                                                                                                             <  1:730367/150‑5 (MQ=255)
        cAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCAGACAGCGGttcgtg                                                                                                                                             <  2:1268721/150‑5 (MQ=255)
             cacTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCg                                                                                                                                        >  1:655665/1‑150 (MQ=255)
             cacTGCGGCCAAAGCGAAAGTGACCCGGCGGACTCCACTCTAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCg                                                                                                                                        <  2:1244008/150‑1 (MQ=255)
                 gctgcCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCa                                                                                                                                    <  2:475340/150‑1 (MQ=255)
                 gctgcCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATGGTTGAGGTACGTCGAAAAGACCGTGCTTGAGGGTGCCGCACGACCCGCTGTCTGCGGTTATATGCCCGATCa                                                                                                                                    >  2:275639/1‑150 (MQ=255)
                  ctgcCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGGACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCaa                                                                                                                                   <  2:919783/150‑1 (MQ=255)
                    gcCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGCCGCACGAACCGCTGTCTGCGGTTATATGCCCGATCATGa                                                                                                                                 >  1:1704013/1‑150 (MQ=255)
                    gcCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGa                                                                                                                                 >  2:200886/1‑150 (MQ=255)
                     ccAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGAt                                                                                                                                >  2:642781/1‑150 (MQ=255)
                     ccAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGAt                                                                                                                                <  2:414110/150‑1 (MQ=255)
                      cAACGCGAAAGGGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATc                                                                                                                               <  1:336625/150‑1 (MQ=255)
                      cAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGAGGTTATATGCCCGATCAAGATc                                                                                                                               >  1:979171/1‑150 (MQ=255)
                       aaaGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATcc                                                                                                                              <  1:819144/150‑1 (MQ=255)
                        aatcGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAACCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCt                                                                                                                             <  2:340275/147‑1 (MQ=255)
                        aaGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCt                                                                                                                             <  2:809598/150‑1 (MQ=255)
                        aaGCGAAAGTGACACGGCGGACTCCACTCGAAAAAAAGTCGATATTGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCt                                                                                                                             <  2:1255292/150‑1 (MQ=255)
                         aGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTg                                                                                                                            <  1:1402189/150‑1 (MQ=255)
                          gCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGc                                                                                                                           <  2:522961/150‑1 (MQ=255)
                             aaaGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaa                                                                                                                        >  1:1359244/1‑150 (MQ=255)
                              aaGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaaa                                                                                                                       <  1:1038533/150‑1 (MQ=255)
                              aaGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaaa                                                                                                                       <  2:585967/150‑1 (MQ=255)
                               cttGACAGGGAGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                                                                                      <  1:1785643/148‑1 (MQ=255)
                               aGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                                                                                      <  1:810864/150‑1 (MQ=255)
                               aGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTGGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                                                                                      <  1:1475468/150‑1 (MQ=255)
                                gTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACg                                                                                                                     >  1:1258178/1‑150 (MQ=255)
                                  gACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGAt                                                                                                                   >  2:1043597/1‑150 (MQ=255)
                                    cacGGCGGACTCCAGTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCg                                                                                                                 <  2:909769/150‑1 (MQ=255)
                                      cggcggACTCCACTCGAACAACAGTCGATAATGACTAAGGCTGATACATTCATGATTGTTGACGTACGTCGAAATGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCggg                                                                                                               <  1:20079/150‑1 (MQ=255)
                                      cggcggACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTTCGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCggg                                                                                                               <  1:292748/150‑1 (MQ=255)
                                          ggACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACACTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACcg                                                                                                           >  1:1215/1‑150 (MQ=255)
                                           gACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACcgc                                                                                                          <  2:759245/150‑1 (MQ=255)
                                            aCTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACcgcg                                                                                                         <  1:1591741/150‑1 (MQ=255)
                                              tCCACTCGAACAAAAGTCGAAAATGACTAAGGCTGAAACATTCATGTTTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATAAGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGa                                                                                                       <  2:477143/150‑1 (MQ=255)
                                               ccACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGAt                                                                                                      <  1:1513751/150‑1 (MQ=255)
                                                  cTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCAt                                                                                                   <  2:1280464/150‑1 (MQ=255)
                                                   tCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATa                                                                                                  <  1:1631409/150‑1 (MQ=255)
                                                     gAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAAGAACCTCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGc                                                                                                >  2:226554/1‑150 (MQ=255)
                                                      aaCAAAAGTCGATAATGACTTAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGcc                                                                                               <  2:230743/150‑1 (MQ=255)
                                                      aaCAAAAGTCGAAAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGcc                                                                                               <  1:1820181/150‑1 (MQ=255)
                                                       aCAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCt                                                                                              >  1:732132/1‑150 (MQ=255)
                                                          aaaGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTaaa                                                                                           >  1:1683753/1‑150 (MQ=255)
                                                              tCGATAATGACTACTGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTgc                                                                                       <  1:178433/150‑1 (MQ=255)
                                                                       ccTAAGGCTGAAACATTGATGATTGTTGAAGTACGTCGAAAAGACCCTGGTTGAGGGTGTCGCACGAACCGCTGTCTGCGGTTTTATGCCCGTTCAAGTTCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAt                                                                              <  1:392164/149‑1 (MQ=255)
                                                                         tAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGAGCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtca                                                                            <  2:676188/150‑2 (MQ=255)
                                                                          aaGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtctt                                                                           <  1:1181069/150‑1 (MQ=255)
                                                                           aGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGTCCCTGGTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTTCAAAACGATCGGGACCGCGGATCATAGCCTACACTGCGCAAGAGAtcttc                                                                          <  2:1009125/150‑1 (MQ=255)
                                                                           aGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttc                                                                          >  1:1753027/1‑150 (MQ=255)
                                                                            ggCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttct                                                                         <  1:1803660/150‑1 (MQ=255)
                                                                             gCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTg                                                                        <  1:275505/150‑1 (MQ=255)
                                                                             gCTGAAACATTCATGATTGTTGACGTACGTAGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTg                                                                        <  1:1848257/150‑1 (MQ=255)
                                                                                  aaCATTCATGATTGTTGANGTACGTCGAAAAGACCCTGNTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtct                                                                   >  2:1061834/1‑150 (MQ=255)
                                                                                   aCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtctc                                                                  <  2:1036168/150‑1 (MQ=255)
                                                                                   aCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtctc                                                                  <  2:353784/150‑1 (MQ=255)
                                                                                   aCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtctc                                                                  <  1:172412/150‑1 (MQ=255)
                                                                                     aTTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTcac                                                                >  1:21891/1‑150 (MQ=255)
                                                                                     aTTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTcac                                                                <  1:96535/150‑1 (MQ=255)
                                                                                        cATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGa                                                             <  2:815057/150‑1 (MQ=255)
                                                                                             ttgttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTc                                                        >  1:1270416/1‑150 (MQ=255)
                                                                                               gttgACGTACGTCGAACAGACCCTGCTTGAGGGTGACGCACGAACCGCGGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTccc                                                      <  1:1853585/150‑1 (MQ=255)
                                                                                               gttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTccc                                                      <  1:1286941/150‑1 (MQ=255)
                                                                                                tagacggacgTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCTGTTATATGCCCGATCTAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCg                                                     <  1:295159/148‑1 (MQ=255)
                                                                                                     gtacgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTa                                                <  1:592912/150‑1 (MQ=255)
                                                                                                         gtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCca                                            >  2:1800476/1‑150 (MQ=255)
                                                                                                         gtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCca                                            >  1:655804/1‑150 (MQ=255)
                                                                                                           cGAAAAGACCCGGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGGTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGAGTTATCcacc                                          >  2:1388245/1‑149 (MQ=255)
                                                                                                            gAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAg                                         >  2:1598626/1‑150 (MQ=255)
                                                                                                            gAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAg                                         >  2:540788/1‑150 (MQ=255)
                                                                                                            gAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAg                                         >  1:1309179/1‑150 (MQ=255)
                                                                                                              aaaGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGa                                       >  1:1657403/1‑150 (MQ=255)
                                                                                                                   cccTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTc                                  <  1:427127/150‑1 (MQ=255)
                                                                                                                    ccTGCTTGAGGGGGACGCACGAACCGCTGTCTGCGGTTATATGCACGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTcc                                 <  2:1189508/150‑1 (MQ=255)
                                                                                                                      tGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAg                               <  2:211523/150‑1 (MQ=255)
                                                                                                                       gCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATACACAGGACTTTCCAGa                              >  2:406714/1‑150 (MQ=255)
                                                                                                                         ttGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAAc                            <  2:538243/150‑1 (MQ=255)
                                                                                                                          tGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACt                           >  2:50141/1‑150 (MQ=255)
                                                                                                                          tGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACt                           >  2:1440889/1‑150 (MQ=255)
                                                                                                                            aGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCg                         >  2:1571551/1‑150 (MQ=255)
                                                                                                                             gggTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGCCCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGc                        >  1:1045292/1‑150 (MQ=255)
                                                                                                                             gggTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGc                        <  1:820992/150‑1 (MQ=255)
                                                                                                                                 gACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAg                    >  1:856246/1‑150 (MQ=255)
                                                                                                                                  ggcatataaCCGCAGACAGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTACCGATTTATCCACAGGACTTTCCAGAACTCGCTAAgt                   >  1:1145561/8‑150 (MQ=255)
                                                                                                                                   cGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAgtg                  <  2:780997/150‑1 (MQ=255)
                                                                                                                                    catataaCCGCAGACAGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAgtgt                 >  1:1423144/6‑150 (MQ=255)
                                                                                                                                     cACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCGTAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTa                >  1:946022/1‑150 (MQ=255)
                                                                                                                                     cACGAACCGCTGGCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTa                <  1:353868/150‑1 (MQ=255)
                                                                                                                                           ccGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGAt          >  1:395862/1‑150 (MQ=255)
                                                                                                                                             gCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGAGCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATcc        <  2:337660/150‑1 (MQ=255)
                                                                                                                                                gTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGc     >  2:968690/1‑150 (MQ=255)
                                                                                                                                                  cTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCcg   >  1:1528271/1‑150 (MQ=255)
                                                                                                                                                   tGCGGTTACATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCcgc  >  2:1684650/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGC  >  NC_000913/3883691‑3883987

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