Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,883,834 T→A 5.0% intergenic (‑105/‑502) 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,8340TA 5.0% 364.5 / 6.3 119intergenic (‑105/‑502)dnaA/rpmHchromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator/50S ribosomal subunit protein L34
Reads supporting (aligned to +/- strand):  ref base T (43/69);  new base A (4/2);  total (47/72)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 119 total aligned reads displayed.

CTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCC  >  NC_000913/3883688‑3883981
                                                                                                                                                  |                                                                                                                                                   
cTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTc                                                                                                                                                  >  1:684344/1‑150 (MQ=255)
cTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTc                                                                                                                                                  >  1:1280464/1‑150 (MQ=255)
   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)
                  ctgctgCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGAt                                                                                                                                >  2:230762/1‑150 (MQ=255)
                  ctgctgCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGAt                                                                                                                                <  2:502284/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)
                       gcCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGa                                                                                                                           <  2:1007281/150‑1 (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)
                          aaaGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATcc                                                                                                                        <  1:1216765/150‑1 (MQ=255)
                           aatcGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAACCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCt                                                                                                                       <  2:340275/147‑1 (MQ=255)
                           aaGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCt                                                                                                                       <  2:809598/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)
                                aaaGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaa                                                                                                                  >  1:353784/1‑150 (MQ=255)
                                 aaGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaaa                                                                                                                 <  2:585967/150‑1 (MQ=255)
                                 aaGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaaa                                                                                                                 <  1:1038533/150‑1 (MQ=255)
                                  cttGACAGGGAGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                                                                                <  1:1785643/148‑1 (MQ=255)
                                  aGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTGGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                                                                                <  1:1475468/150‑1 (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)
                                                 tCCACTCGAACAAAAGTCGAAAATGACTAAGGCTGAAACATTCATGTTTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATAAGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGa                                                                                                 <  2:477143/150‑1 (MQ=255)
                                                  ccACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGAt                                                                                                <  1:1513751/150‑1 (MQ=255)
                                                    aCTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCa                                                                                              <  1:859133/150‑1 (MQ=255)
                                                     cTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCAt                                                                                             <  2:1280464/150‑1 (MQ=255)
                                                      tCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATa                                                                                            <  1:1631409/150‑1 (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)
                                                                     tAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCaa                                                                             >  2:1722935/1‑150 (MQ=255)
                                                                            tAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGAGCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtca                                                                      <  2:676188/150‑2 (MQ=255)
                                                                             aaGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtctt                                                                     <  1:1181069/150‑1 (MQ=255)
                                                                              aGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGTCCCTGGTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTTCAAAACGATCGGGACCGCGGATCATAGCCTACACTGCGCAAGAGAtcttc                                                                    <  2:1009125/150‑1 (MQ=255)
                                                                               ggCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttct                                                                   <  1:1803660/150‑1 (MQ=255)
                                                                                     aaCATTCATGATTGTTGANGTACGTCGAAAAGACCCTGNTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtct                                                             >  2:1061834/1‑150 (MQ=255)
                                                                                      aCATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtctc                                                            <  1:172412/150‑1 (MQ=255)
                                                                                        aTTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTcac                                                          >  1:21891/1‑150 (MQ=255)
                                                                                         ttCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCATGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCATGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTcaca                                                         <  2:130710/150‑1 (MQ=255)
                                                                                           cATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGa                                                       <  2:815057/150‑1 (MQ=255)
                                                                                            aTGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGAt                                                      >  1:1708180/1‑150 (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)
                                                                                                       cgtacgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGAttt                                           >  1:15372/1‑150 (MQ=255)
                                                                                                        gtacgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTa                                          <  1:592912/150‑1 (MQ=255)
                                                                                                            gtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCca                                      >  2:1800476/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)
                                                                                                                aaaaGACCCTGCTTGAGGGTGACGCACGAACCTCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAgg                                  <  1:784008/150‑1 (MQ=255)
                                                                                                                 aaaGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGa                                 >  1:1657403/1‑150 (MQ=255)
                                                                                                                     aCCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACttt                             <  2:655665/150‑1 (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:1440889/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)
                                                                                                                                    gACGCACGAACCGCTGTCTGCGGTTATATGCCCCATCAAGATCCTGCAAAACGATCGGGACCGCGGGTCATAGCCTAAACTGCGCAAGAGATCTTATGTTTCTCACAGATTATCCCGATTTATCCACAGGACTTTCCAGACCTCGCTAAg              >  2:1162985/1‑150 (MQ=255)
                                                                                                                                     ggcatataaCCGCAGACAGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTACCGATTTATCCACAGGACTTTCCAGAACTCGCTAAgt             >  1:1145561/8‑150 (MQ=255)
                                                                                                                                      cGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAgtg            <  2:780997/150‑1 (MQ=255)
                                                                                                                                       catataaCCGCAGACAGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAgtgt           >  1:578427/6‑150 (MQ=255)
                                                                                                                                       catataaCCGCAGACAGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAgtgt           >  1:1423144/6‑150 (MQ=255)
                                                                                                                                       catataaCCGCAGACAGCGGTTATATGCCAGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAGCTTCTGTTTCTCACAGATTTTCCAGATTTATCCAGAGGACTTTCCAGCACTCGCTAAgtgt           >  1:262586/6‑150 (MQ=255)
                                                                                                                                        cACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCGTAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTa          >  1:946022/1‑150 (MQ=255)
                                                                                                                                        cACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTa          <  1:1820062/150‑1 (MQ=255)
                                                                                                                                        cACGAACCGCTGGCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTa          <  1:353868/150‑1 (MQ=255)
                                                                                                                                              ccGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGAt    >  1:395862/1‑150 (MQ=255)
                                                                                                                                               cGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATc   <  2:382425/150‑1 (MQ=255)
                                                                                                                                                gCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGAGCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATcc  <  2:337660/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCC  >  NC_000913/3883688‑3883981

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: