Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,882,515 A→C 5.8% R405R (CGT→CGG)  dnaA ← chromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,882,5150AC 5.8% 241.9 / 3.1 68R405R (CGT→CGG) dnaAchromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base C (2/2);  total (37/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 = 1.00e+00

GATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGAAG  >  NC_000913/3882368‑3882659
                                                                                                                                                   |                                                                                                                                                
gATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATATGGCGCGGAc                                                                                                                                                 >  3:5605/1‑149 (MQ=255)
    tCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACgcgc                                                                                                                                             >  1:533875/1‑149 (MQ=255)
     cGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGTCCGCGcc                                                                                                                                            <  2:16604/149‑1 (MQ=255)
      gTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCa                                                                                                                                           >  2:365192/1‑149 (MQ=255)
           tctctTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAg                                                                                                                                      >  2:656626/1‑149 (MQ=255)
                 cACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGAt                                                                                                                                >  2:433190/1‑149 (MQ=255)
                  aCGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATc                                                                                                                               <  2:218039/149‑1 (MQ=255)
                  aCGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATc                                                                                                                               >  3:541210/1‑149 (MQ=255)
                    gCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGa                                                                                                                             >  4:360377/1‑149 (MQ=255)
                       aCTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGc                                                                                                                          >  4:447444/1‑149 (MQ=255)
                         tGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCtt                                                                                                                        <  2:533875/149‑1 (MQ=255)
                           cTCGATCTTACGGCAGGCGGGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTgg                                                                                                                      >  2:113168/1‑149 (MQ=255)
                            tCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGa                                                                                                                     >  2:654208/1‑149 (MQ=255)
                            tCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGa                                                                                                                     >  2:395170/1‑149 (MQ=255)
                            tCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGa                                                                                                                     >  3:372887/1‑149 (MQ=255)
                              gATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGaaa                                                                                                                   <  2:490947/149‑1 (MQ=255)
                              gATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGaaa                                                                                                                   <  1:378296/149‑1 (MQ=255)
                              gATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGaaa                                                                                                                   <  4:372887/149‑1 (MQ=255)
                                  ttACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGgag                                                                                                               >  2:553442/1‑149 (MQ=255)
                                  ttACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGACATCTGGCGCGGACGCGCCACCGAGCGAGATCGACGCTTGCAAAGgag                                                                                                               >  1:533225/1‑149 (MQ=255)
                                    aCGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGAt                                                                                                             >  3:314917/1‑149 (MQ=255)
                                        caggcaTGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCcgc                                                                                                         >  4:271662/1‑149 (MQ=255)
                                           gcaTGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGAc                                                                                                      >  3:497370/1‑149 (MQ=255)
                                                 aGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCAGCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGAt                                                                                                <  4:5605/149‑1 (MQ=255)
                                                  gCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATc                                                                                               <  1:395170/149‑1 (MQ=255)
                                                  gCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATc                                                                                               <  1:656626/149‑1 (MQ=255)
                                                    aCCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCtt                                                                                             >  2:579289/1‑149 (MQ=255)
                                                     tggtcgtGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTg                                                                                            >  1:339158/3‑149 (MQ=255)
                                                      cgtcgtGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgt                                                                                           <  4:497370/149‑1 (MQ=255)
                                                               gTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTCTGATCTTGTAGTACTCCg                                                                                  >  2:72400/1‑149 (MQ=255)
                                                                tCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGc                                                                                 >  4:42944/1‑149 (MQ=255)
                                                                  aCGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCa                                                                               >  1:659263/1‑149 (MQ=255)
                                                                     gCCACGAAACGCATCGCCAATCTCCGGCAGACTCTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCg                                                                            >  1:196835/1‑149 (MQ=255)
                                                                       caccaAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtc                                                                          >  2:482603/1‑149 (MQ=255)
                                                                         ccaAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtctt                                                                        <  4:386185/149‑1 (MQ=255)
                                                                          caAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAACGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtcttc                                                                       <  1:365192/149‑1 (MQ=255)
                                                                           aaaCGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtcttct                                                                      <  1:678465/149‑1 (MQ=255)
                                                                              cGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGaa                                                                   <  1:579289/149‑1 (MQ=255)
                                                                                 aTCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAatat                                                                >  1:373388/1‑149 (MQ=255)
                                                                                  tCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATAtt                                                               <  3:217780/149‑1 (MQ=255)
                                                                                    gCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGt                                                             >  3:24347/1‑149 (MQ=255)
                                                                                      cAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgt                                                                                           >  2:366814/1‑117 (MQ=255)
                                                                                      cAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgt                                                                                           <  1:366814/117‑1 (MQ=255)
                                                                                         tctcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATg                                                        >  4:108888/1‑149 (MQ=255)
                                                                                         tctcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATg                                                        >  2:310406/1‑149 (MQ=255)
                                                                                           tcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGt                                                      <  2:196835/149‑1 (MQ=255)
                                                                                             cGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGa                                                    >  2:350414/1‑149 (MQ=255)
                                                                                                     tgtgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGAt                                                         >  3:244834/1‑136 (MQ=255)
                                                                                                     tgtgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGAt                                                         <  4:244834/136‑1 (MQ=255)
                                                                                                       tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTc                                          <  1:268688/149‑1 (MQ=255)
                                                                                                       tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCAAACGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAAGATTGTCGATGGTGACCAGTTTTTc                                          <  1:113168/149‑1 (MQ=255)
                                                                                                           tAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGc                                      <  3:360377/149‑1 (MQ=255)
                                                                                                            aGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGCTCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCa                                     <  4:541210/149‑1 (MQ=255)
                                                                                                              tCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAAt                                   >  2:412294/1‑149 (MQ=255)
                                                                                                                  ctctTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCa                               >  3:498393/1‑149 (MQ=255)
                                                                                                                     tttCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGc                                 >  2:523574/1‑144 (MQ=255)
                                                                                                                     tttCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGc                                 <  1:523574/144‑1 (MQ=255)
                                                                                                                          ccAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTcg                       <  2:551669/149‑1 (MQ=255)
                                                                                                                                ccatcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAgc                 <  1:72400/149‑1 (MQ=255)
                                                                                                                                ccatcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAgc                 <  3:108888/149‑1 (MQ=255)
                                                                                                                                     gccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTc            <  4:314917/149‑1 (MQ=255)
                                                                                                                                     gccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTc            <  2:373388/149‑1 (MQ=255)
                                                                                                                                          gtGGCGCGTCCGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTcacgca       >  4:59377/2‑149 (MQ=255)
                                                                                                                                           tGGCGCGTCCGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTc                                                            <  4:117293/95‑1 (MQ=255)
                                                                                                                                           tGGCGCGTCCGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTc                                                            >  3:117293/1‑95 (MQ=255)
                                                                                                                                             gcgcGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTcc                                                                                   >  4:331995/1‑70 (MQ=255)
                                                                                                                                             gcgcGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTcc                                                                                   <  3:331995/70‑1 (MQ=255)
                                                                                                                                               gcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGAAg  <  3:271662/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGAAG  >  NC_000913/3882368‑3882659

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑