Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,883,678 G→A P18S (CCA→TCA)  dnaA ← chromosomal replication initiator protein DnaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,883,6780GA97.1% 199.0 / ‑2.1 68P18S (CCA→TCA) dnaAchromosomal replication initiator protein DnaA
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (27/39);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

GCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGGTAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCT  >  NC_000913/3883434‑3883922
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
gCAGGGGCCGCGACTTGGCGCGCCACTGCCGCTTGTGGCGTTTGCTTACCCGGTTTGGGGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGGTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACt                                                                                                                                                                                                                                                 <  2:400218/250‑1 (MQ=255)
  aGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCAt                                                                                                                                                                                                                                              >  1:169168/1‑251 (MQ=255)
       ccGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCAATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAACACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATCACTCATCCTGc                                                                                                                                                                                                                                         >  2:171130/1‑251 (MQ=255)
              gTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGACGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAATACGGTTTGGCGCGTCCAGGCCCAGCGTGTTATCGCTCAGTTCCGCCTGCACTGGGCGTATCCACATACTGCATTCTGTGGCTGATAACTCATCCTGCAATCggg                                                                                                                                                                                                                                  >  1:284127/1‑251 (MQ=255)
              gTTGCTCGTCACTGCCGCTTGAGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTTATCACTCATCCTGCAATCggg                                                                                                                                                                                                                                  >  2:336857/1‑251 (MQ=255)
                 gCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCa                                                                                                                                                                                                                                >  1:62393/1‑250 (MQ=255)
                    cGTCACTGCCGCGTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCAGCCTGCAATCGGGCAAGa                                                                                                                                                                                                                             <  1:358027/250‑1 (MQ=255)
                           gCCGCTTGCGGCGTTTGAGTCACCGGTTTGGTGCCGACTTCACAAAGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATAGAGGACAAAACGGTTTGGCGCGTACTGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCACTGGGCGTATCCACATACTGACTTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAAACActgctg                                                                                                                                                                                                                     >  2:208600/1‑251 (MQ=255)
                           gCCGCTGGTGGCGGGTGCGTCCCCGGTTGTGGGCCTACCTCAAAACGCAGCTGTGGGGCCTCCGCTCCGCAGAAACTGGATAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAAGCGAGGACAAAACGGTTTGGCGCGTACAGGGCCATCGTGTTATCGCTCAGTTCCGTCCGCAATGGGCGTATCCACCTACTGAATTCTGTGGCTGATAACTCAGCCTGCAATCGGGCAAGCCActgc                                                                                                                                                                                                                       <  2:91493/249‑1 (MQ=255)
                                     gCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGctaa                                                                                                                                                                                                           <  2:169168/251‑4 (MQ=255)
                                                aCCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACAcggc                                                                                                                                                                                                >  1:459296/1‑251 (MQ=255)
                                                    gTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCACTGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGAc                                                                                                                                                                                            >  1:284333/1‑251 (MQ=255)
                                                         ggtgctgactcaaaaaCGCATCTTTGTGGCAACCGCTCCGCAGAATCTGGTTATCAGTCCAGTGATATTATTAAGGTATGTGTCCCGTACCCAATCGAGTATATAACGGTTTGGCGCGTACAGGGCCAGCTTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCCCATACTGAATTCTGTGGCTGATAACTGATCCTGCCAGCGGGCAAGACTCTGCTGCCAAAGCGAAAGTGACACGGCCGAcgcc                                                                                                                                                                                         <  2:318385/238‑4 (MQ=255)
                                                         ggTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCa                                                                                                                                                                                        <  1:159108/251‑1 (MQ=255)
                                                              ccacattaaaacggtggtgtgGGGCTTCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGa                                                                                                                                                                                   <  2:340441/234‑1 (MQ=255)
                                                              cGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGa                                                                                                                                                                                   <  1:62669/251‑1 (MQ=255)
                                                              cGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTAAAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCTAAGCGAAAGTGACACGGCGGACTCCACTCg                                                                                                                                                                                    >  1:208241/1‑250 (MQ=255)
                                                                   tCAAAACGCAGTTGTGGGTCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACa                                                                                                                                                                                <  2:138881/249‑1 (MQ=255)
                                                                   tCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACaa                                                                                                                                                                               >  1:201728/1‑250 (MQ=255)
                                                                    cAAAACGCAGGTGTGGGGCTTCCGCACCGCAGAAACTAGTTAGCAGTCCAGTGATATTATTAAGGTGCTTGTCCCGTACCCAAACGAGTACAAAACGGTTTGGCGCGTACAGGGTCAGCGTGTTATCGCTGATGTCCGCCTGCAATGGGTGTAGCCACATACTGAATTCTGTGTCTGATAACTCATCCTGCTATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACaaaa                                                                                                                                                                             <  2:158787/251‑1 (MQ=255)
                                                                     aaaaCGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAg                                                                                                                                                                            <  1:240727/251‑1 (MQ=255)
                                                                      aaaCGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATCTTCTAAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGATGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAACGAAAAAGACACGGCGGACACCAATCGCAAAAAAGt                                                                                                                                                                           >  1:439340/1‑251 (MQ=255)
                                                                       aacgcgtcttttgggCATCCGCTCCGGAGAAACTGGTTAGCAGTCTATTGATTGTATTAGGGTGCTGGTCACGTGCCCATTCTAGGACAAAACGGTTTAGCGCGTGCAGGGCCAGCGTGTTATCGCGCAGTTCCGCCTGCGATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCAGCCTGCAATCGGGCAAGCCACTGCTGCCCAAGCGAAAGTGACTCGGCGGACTCCTCTCGAACCAAAGtt                                                                                                                                                                          <  2:301926/239‑2 (MQ=255)
                                                                                  ggggCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATAGTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGTCCAGCGTGTTATCGCGCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGATTCCACTCGAACAAAAGTCGATAATGACt                                                                                                                                                                <  1:172578/250‑1 (MQ=255)
                                                                                   gggCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTACGGTACTTGTCCCGTCCCCAATCGAGGACAAAACGGTTTGGCGCGTACCGGGCCAGCGTGTTATCGCTAAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAAACCATCTTGACATCGGGCAAGACACTGCTGCCCAAGCGAAAGAGACACGGCGGCCTCCAATCGAAccaacgccggtaatgccta                                                                                                                                                               >  2:41791/1‑232 (MQ=255)
                                                                                   gggCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGAAAATGACTa                                                                                                                                                               >  1:55069/1‑250 (MQ=255)
                                                                                          cgctccgcAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGa                                                                                                                                                        >  1:80946/1‑250 (MQ=255)
                                                                                              ccgcATAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGGCACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACa                                                                                                                                                    <  2:427198/250‑1 (MQ=255)
                                                                                              ccgcAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAAc                                                                                                                                                     <  2:37579/249‑1 (MQ=255)
                                                                                                  aGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTTCCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCAt                                                                                                                                               <  2:337705/251‑1 (MQ=255)
                                                                                                      aCTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGAtt                                                                                                                                           <  1:336857/251‑1 (MQ=255)
                                                                                                      aCTGGTTAACAGACCATTGATATTATTAAGGTACTTGGCCCGTTCCCAATCGAGGACAAAACGGTTGGGCGCGGACAGTGCCAGCGTGTTATCGCTCATTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCTCTCGAACAAAAGTCTATAATGACTAAGGCTGAAACATTCATGAtt                                                                                                                                           <  2:5543/251‑1 (MQ=255)
                                                                                                                 gTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCCCTCGAACAAAAGTCGATAATGACTAAGGCTGAAACAtt                                                                                                                                                  >  2:202840/1‑233 (MQ=255)
                                                                                                                 gTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACAtt                                                                                                                                                  <  1:202840/233‑1 (MQ=255)
                                                                                                                   ccATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGacgtacgt                                                                                                                               <  1:56084/250‑1 (MQ=255)
                                                                                                                       tGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGAACGTCGaa                                                                                                                           >  1:321973/1‑250 (MQ=255)
                                                                                                                              attaAGGTACTTGTCCCGTTCCCAATCGAGGACAAAACGGTTTGGCGCGTACATGGCCAGCGTGTTATCGATCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGAccc                                                                                                                    <  1:121810/250‑1 (MQ=255)
                                                                                                                                  aGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGGTAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGAAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCttt                                                                                                              >  2:107627/1‑250 (MQ=255)
                                                                                                                                         tGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCCCTCGAACAAAAGGCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTg                                                                                                        >  2:382031/1‑251 (MQ=255)
                                                                                                                                               gTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAc                                                                                                  <  1:78182/251‑1 (MQ=255)
                                                                                                                                               gTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAc                                                                                                  <  2:85853/251‑1 (MQ=255)
                                                                                                                                                   ccAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAAc                                                                                              <  1:414960/251‑1 (MQ=255)
                                                                                                                                                    cAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAAcc                                                                                             <  2:459296/251‑1 (MQ=255)
                                                                                                                                                        cGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGAGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGGACAAAAGTCGCTAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAGGACCCTGCTTGAGGGTGACGCACGAACCGCTg                                                                                         >  2:175301/1‑251 (MQ=255)
                                                                                                                                                           ggACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCt                                                                                      <  1:244992/251‑1 (MQ=255)
                                                                                                                                                                   cGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCATGCTTGAGGGTGACGCCCGAACCGCTGTCTGCGGTtat                                                                              >  2:201575/1‑251 (MQ=255)
                                                                                                                                                                     gTTTGGCGGGTGCAGGGCCTGCGTGTGATCGCTAATTCCCGCCTGCAAGGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTGATCCTGCAATCGGGCAAGACACTGCGGCCAAAGCGATAGTGACACGGCGTACTCCACTCGATCAAAAGTCGATAATGACTAAGTCTGAAACATTCATGATCGTTGCCGGACGTCGAAAAGACCCTGCTTGAGGGTTACGCACGAACCGCTGTTTGCGGTTatat                                                                            <  2:48366/251‑1 (MQ=255)
                                                                                                                                                                     gTTTGGCGCTTACAGGGCCAGCGTGTGATCGCTCAGTTCCGCCTTCAATGGGCGTATCCACATACGTAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAATACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCAGTTatat                                                                            <  2:394561/251‑1 (MQ=255)
                                                                                                                                                                        tGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGc                                                                          <  2:84912/250‑1 (MQ=255)
                                                                                                                                                                        tGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGc                                                                          <  1:458483/250‑1 (MQ=255)
                                                                                                                                                                           cgcgTACAGAGCAAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGACTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAACGAAAGTGAGACGGCGGACTCCCCTCGAACACGAGTCGATAATGCCTAAAGCTGAACCATTCATGATTGTTGACGTACGTCGAAACGACCCTGCTTGGGAGTGACGCACGAACCGCTGTTTGCGGTGATATGCCCgc                                                                      >  2:388213/1‑250 (MQ=255)
                                                                                                                                                                                   gggCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGTGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATc                                                              <  1:290936/251‑1 (MQ=255)
                                                                                                                                                                                     gCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAACGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATAATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCt                                                            >  1:18144/1‑251 (MQ=255)
                                                                                                                                                                                          cGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaa                                                       <  1:408953/251‑1 (MQ=255)
                                                                                                                                                                                            tgtTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGGTAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                     >  1:93055/1‑251 (MQ=255)
                                                                                                                                                                                             gtTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGACTTCTGTGGCTGATAACTCATCCTGCAACCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGGGGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCTCTGTCTGCGGTTATATGCCCGATCAAGATCCTGTAAAAc                                                     <  1:167567/250‑1 (MQ=255)
                                                                                                                                                                                                          tCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGGACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGGTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCgg                                       >  1:230208/1‑251 (MQ=255)
                                                                                                                                                                                                                  cAATGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCAGCCTGCAATCGGGCAAGACACTTCTGACAAAGCGAAAGTGACACGGCGGACTTCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGGTTGAGGGTGACGCACGAACTGCTGTCTGCGGTTATAGGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGc                               <  2:426787/251‑1 (MQ=255)
                                                                                                                                                                                                                  cAATGGGCGTATCCACATACTGAAGTCGGGGGCTGATAACTCATCCTGCAAGCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGc                               <  2:468084/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTGGGCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGAGCCGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCCTAGCCt                             >  2:337852/1‑251 (MQ=255)
                                                                                                                                                                                                                     tgagCGTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGCCTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGGGACGCACGAACAGCTGTCTGCGGTTCTATGCACGGTCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCt                             >  1:435559/4‑250 (MQ=255)
                                                                                                                                                                                                                          gTATCCACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTg                       >  1:443170/1‑251 (MQ=255)
                                                                                                                                                                                                                              ccACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCgg                                                                                  <  2:319443/188‑1 (MQ=255)
                                                                                                                                                                                                                              ccACATACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCgg                                                                                  >  1:319443/1‑188 (MQ=255)
                                                                                                                                                                                                                               cacaTACTGAATTCTGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCaa                  <  1:463645/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aCTGAATTCTGGGGCTGATAACTCAGCCTGCAATCGGGCAAGACACTGCGGCCAAAGCGAAAGTGACACGGCGGACTCCAACCGAACAAAAGTCGATAATGACTAAGGCTTAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAt             <  2:183593/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gtgGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGttt    <  1:171130/250‑1 (MQ=255)
                                                                                                                                                                                                                                               tgGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtct  >  1:96502/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
GCAGGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTCTGTGGCTGGTAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCT  >  NC_000913/3883434‑3883922

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