Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,174,035 T→G E454D (GAA→GAC parE ← DNA topoisomerase IV, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,174,0350TG97.1% 230.6 / NA 70E454D (GAA→GACparEDNA topoisomerase IV, subunit B
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (37/31);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGATCAATACGGTAGAGCGGTGGCAGTGCGACGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTGGTCAGCTTTTTACGCACCACTTTTTTGGCCGCACGCATACGGCGCT  >  NC_000913/3173804‑3174267
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
gagaTCAATACGGTAGGGCGGTGGCAGTGCAACGTAAACGTGACCGGGTTTCACCAACGCGCGGAAATGTTTGACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCAcgtcgt                                                                                                                                                                                                                                      <  2:527863/236‑1 (MQ=255)
         gcGGTAGAGCGGTGGCAGTGCGACGGAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAATGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAgaga                                                                                                                                                                                                                             <  2:660765/235‑1 (MQ=255)
         aCGGTAGAGCGGTGGCAGTGCGACGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAgaga                                                                                                                                                                                                                             <  1:265102/236‑1 (MQ=255)
         aCGGTAGAGCGGTGCCAGTGCGACGTAAACGTGACCGTGTTTCCCCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAgaga                                                                                                                                                                                                                             <  2:271591/236‑1 (MQ=255)
             tAGAGCGGTGGCAGTGCGACGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGGACTTGCTGCGAAGCCAGCACGTCGTCGGaaatgacgtc                                                                                                                                                                                                                         >  2:420422/1‑228 (MQ=255)
             tAGAGCGGTGGCAGTGCGACGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGACCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGAGTCCTACCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTc                                                                                                                                                                                                                         >  2:387190/1‑236 (MQ=255)
                ggCGGTTGCAGTTCAACGTAAACGTGACCGTGTTTCACCAACGGGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTccc                                                                                                                                                                                                                       <  2:250992/234‑1 (MQ=255)
                agCGGTGGCAGTGCGACGTAACCGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTccc                                                                                                                                                                                                                       <  1:339550/235‑1 (MQ=255)
                agCGGTGGCAGTGCGACGTAAACGTGACCGTGTTTCACCAACGCTCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTccc                                                                                                                                                                                                                       <  2:651609/235‑1 (MQ=255)
                agCGGTGGCAGTGCGACGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTccc                                                                                                                                                                                                                       <  1:127799/235‑1 (MQ=255)
                agCGGTGGCAGTGCGACGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTccc                                                                                                                                                                                                                       <  1:1763/235‑1 (MQ=255)
                agCGGTGGCAGTGCGACGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTccc                                                                                                                                                                                                                       <  1:260568/235‑1 (MQ=255)
                     tGGCAGTGCGACGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGgtgt                                                                                                                                                                                                                 >  1:160967/1‑236 (MQ=255)
                       gCAGTGCGACGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGtt                                                                                                                                                                                                                <  2:746615/235‑1 (MQ=255)
                       gCAGTGCGACGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTa                                                                                                                                                                                                               >  1:27443/1‑236 (MQ=255)
                               aCGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTccc                                                                                                                                                                                                                       >  2:728115/1‑220 (MQ=255)
                               aCGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTccc                                                                                                                                                                                                                       <  1:728115/220‑1 (MQ=255)
                               aCGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTaa                                                                                                                                                                                                              <  1:387190/229‑1 (MQ=255)
                               aCGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCt                                                                                                                                                                                                        <  2:634565/235‑1 (MQ=255)
                                    aaCGTGACCGTGTGTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCtt                                                                                                                                                                                                  <  2:397983/236‑1 (MQ=255)
                                          aCCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTg                                                                                                                                                                                            >  1:18574/1‑236 (MQ=255)
                                               ttttacccAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATg                                                                                                                                                                                       <  2:523858/230‑1 (MQ=255)
                                               gtTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCATTGGCAt                                                                                                                                                                                        >  1:523858/1‑235 (MQ=255)
                                                         cgcgcgGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATAATCGCCTGa                                                                                                                                                                              >  1:484982/1‑235 (MQ=255)
                                                                      tttACGAACAAAGCGAAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATcgcg                                                                                                                                                                <  2:160967/236‑1 (MQ=255)
                                                                          cGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGGGTTAAGGATCTTACCTTTCAGTGGCATGATCGACTGATATTCTCGATcgcgcg                                                                                                                                                              >  2:222345/1‑234 (MQ=255)
                                                                                        agcagcGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATGGCCTGATATTCTCGATCGCGCGCCTGCTTGGCAGATcc                                                                                                                                              >  1:537486/1‑236 (MQ=255)
                                                                                                  cAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGGGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGa                                                                                                                                     >  1:580614/1‑235 (MQ=255)
                                                                                                       tgCAGACCATCAGAGTCCGCATCCGCGATGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCAcc                                                                                                                               <  2:724155/236‑1 (MQ=255)
                                                                                                       tgCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCCGATCCGCCTGCGCAGTAAc                                                                                                                                >  1:279043/1‑235 (MQ=255)
                                                                                                       ggCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACCAAGCTGGCTCAGTTCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCAc                                                                                                                                <  2:630694/234‑1 (MQ=255)
                                                                                                                      tCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAg                                                                                                                >  2:667482/1‑236 (MQ=255)
                                                                                                                        cGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAAc                                                                                                                  >  1:250568/1‑232 (MQ=255)
                                                                                                                            tCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGt                                                                                                          <  2:27443/236‑1 (MQ=255)
                                                                                                                            tCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGt                                                                                                          <  1:674671/236‑1 (MQ=255)
                                                                                                                              cgcgAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCATCCCCACGGCACAGCTCGGTAc                                                                                                        >  2:674671/1‑236 (MQ=255)
                                                                                                                                     aTACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTCCGGATCCACCTTCCCCAAGGAACAGCCCGGTACGGTTaa                                                                                                  >  1:182591/1‑235 (MQ=255)
                                                                                                                                       aCAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGGACGGTTaa                                                                                                  >  1:531822/1‑233 (MQ=255)
                                                                                                                                             tttGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACACGGAACAGCTCGGTATGGTTAAGGGCCTgc                                                                                          >  2:203597/1‑235 (MQ=255)
                                                                                                                                             tttGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTAc                                                                                                        <  2:690277/221‑1 (MQ=255)
                                                                                                                                             tttGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTAc                                                                                                        >  1:690277/1‑221 (MQ=255)
                                                                                                                                                   aTAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCtt                                                                                                                             <  1:435744/194‑1 (MQ=255)
                                                                                                                                                     aaCGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTgcg                                                                                         <  2:644066/228‑1 (MQ=255)
                                                                                                                                                     aaCGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTgcg                                                                                         >  1:644066/1‑228 (MQ=255)
                                                                                                                                                          aGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTTGCACATCCGCCTGCGGAGTGACCTTCCCCAAg                                                                                                                      >  2:435744/1‑194 (MQ=255)
                                                                                                                                                                     aTCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCgg                                                                                                           <  1:117254/194‑1 (MQ=255)
                                                                                                                                                                     aTCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCgg                                                                                                           >  2:117254/1‑194 (MQ=255)
                                                                                                                                                                     aTCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAg                                                                 >  2:223743/1‑236 (MQ=255)
                                                                                                                                                                         tcgCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTGCGGTTAAGGTCCTgc                                                                                          >  1:290841/1‑207 (MQ=255)
                                                                                                                                                                         tcgCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTGCGGTTAAGGTCCTgc                                                                                          <  2:290841/207‑1 (MQ=255)
                                                                                                                                                                         tcgCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCCGCGCGGTACAATCAGCCAGTTTGCCAGGCa                                                              >  1:724894/1‑235 (MQ=255)
                                                                                                                                                                                         tACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGATGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTAACCTTCCACAAGGAACAGCCCGGTACGGTTAAGGGCCTGGCCGGAACAATCAGCAAGATTGCCAGGCAACGCCGGGCCGTTgg                                               >  2:705322/1‑234 (MQ=255)
                                                                                                                                                                                         tACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCGGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTgg                                               >  1:578483/1‑234 (MQ=255)
                                                                                                                                                                                         tACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCGTGATCGCCTGATATTCCCGATCCCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGGACGGGTAAGGGCCTGCCCGGACCAATCAGCCAGGTTGCAAGGCAAGTCCCGGCCGCTCGTc                                             >  2:159749/1‑236 (MQ=255)
                                                                                                                                                                                         tACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTgg                                               >  1:486183/1‑234 (MQ=255)
                                                                                                                                                                                                   tACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTGGTCAGCTTTTTa                                    >  1:81746/1‑235 (MQ=255)
                                                                                                                                                                                                               gTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTa                                                                                                   <  2:663787/160‑1 (MQ=255)
                                                                                                                                                                                                               gTGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTa                                                                                                   >  1:663787/1‑160 (MQ=255)
                                                                                                                                                                                                                tGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTGGTCa                                            >  1:268245/1‑214 (MQ=255)
                                                                                                                                                                                                                tGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTGGTCa                                            <  2:268245/214‑1 (MQ=255)
                                                                                                                                                                                                                tGCACTTCCTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGACGTACAATCAGCCAGTTTGCCAGGCAACTCCGGGCCGCTGGTCAGCCTTTTACGGACCACt                           >  2:412416/1‑231 (MQ=255)
                                                                                                                                                                                                                   ccTTCCTGCGAAGCCAGCACTGCATCGGAAGAGACTTCCCTGTTGTTAAGGACCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTGGTCAGCTTTTTACGCACCACTTTTTTGGc                   <  2:182591/235‑1 (MQ=255)
                                                                                                                                                                                                                    ctcacTGCAAAGCCAGCACGTCGTCGGAAGAGACTGCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGACCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTGGTCAGCTTTTTACGCACCACTTTTTTGGc                   <  2:486183/231‑1 (MQ=255)
                                                                                                                                                                                                                      tcatgagAAGCCAGCACGTCGGCGGAAGAGGCTTCCCAGGGTTTAAGGATCTTCCCTTTCAGTTGCATGCACGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTGGTCAGCTTTTTACGCACCACTTTTTTGGCcg                 <  2:531822/229‑1 (MQ=255)
                                                                                                                                                                                                                       ccTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTGGt                                              >  2:676286/1‑205 (MQ=255)
                                                                                                                                                                                                                       ccTGCGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTGGt                                              <  1:676286/205‑1 (MQ=255)
                                                                                                                                                                                                                         atcgaatgcAGCACGGCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGGGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTGGTCAGCTTTTTACGCACCACTTTTTTGGCcg                 <  2:279043/224‑1 (MQ=255)
                                                                                                                                                                                                                           cGAAGCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCACATTGCCAGGCAACGCCGGGCCGCTGGTCAGCTTTCTACTCACCACCTTTTTCAccgcccgc            >  1:368889/1‑235 (MQ=255)
                                                                                                                                                                                                                               gCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTGGTCAGCTTTTTACGCACCACTTTTTTGGCCGCACGCATAc        >  2:580773/1‑235 (MQ=255)
                                                                                                                                                                                                                               gCCAGCACGTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGCACAATCAGCCAGCTTGCCAGGCAACGCCGGGCCGCTGGTCAGCCTTTTACGCACCACTTTTTTGTCCGCAAGCAttc        >  2:545926/1‑233 (MQ=255)
                                                                                                                                                                                                                                    cACTTCGTCGGAAGAGACTTCACTGGTGTTAAGGATCTTACCTTTCAGTGGCATGAACGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTGGTCAGCTTTTTACGCACCACTTTTTTGGCCGCACGCATACggcgct  <  2:391929/236‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
GAGATCAATACGGTAGAGCGGTGGCAGTGCGACGTAAACGTGACCGTGTTTCACCAACGCGCGGAAATGTTTTACGAACAAAGCGCAGAGCAGCGTGGCAATGTGCAGACCATCAGAGTCCGCATCCGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTCCCAGGTGTTAAGGATCTTACCTTTCAGTGGCATGATCGCCTGATATTCGCGATCGCGCGCCTGCTTGGCAGATCCGCCTGCGGAGTCACCTTCCACAAGGAACAGCTCGGTACGGTTAAGGTCCTGCGCGGTACAATCAGCCAGTTTGCCAGGCAACGCCGGGCCGCTGGTCAGCTTTTTACGCACCACTTTTTTGGCCGCACGCATACGGCGCT  >  NC_000913/3173804‑3174267

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