Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,662,163 A→C N584H (AAC→CAC)  mdtF → anaerobic multidrug efflux transporter, ArcA‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,662,1630AC98.4% 204.8 / ‑4.1 64N584H (AAC→CAC) mdtFanaerobic multidrug efflux transporter, ArcA‑regulated
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (30/33);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

GCCCTGTTCGCACGCTTCAACACGCTGTTTGAAAAATCAACTCAACACTATACCGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTAACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCGGGCAATGATTGCGTTAAGCAGTATCAATAAAGC  >  NC_000913/3661938‑3662389
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  
ccacTGTTCCGCACGGTTCAACCAGTAGATGGAAAACTCAACGCAACGCTGTACAGATCGCACCCGCTCGTTGTTGGGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGa                                                                                                                                                                                                                           <  2:436343/233‑1 (MQ=255)
    gtttCGCACGTTTACACCCGCTTATTGAAAATTCAACCCAACGCTAGCCCGATTGCACCCGCTCGTTGTTGCCTTGTACCGGTCGCTACATGGTTGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGtg                                                                                                                                                                                                                      <  2:668707/234‑1 (MQ=255)
                 tCAACACGCTGTTTGAAAAATCAACTCAACACTATACCGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTg                                                                                                                                                                                                          <  1:354832/236‑1 (MQ=255)
                   aaCGCGCTGTTTGAAACATCAACTCAACCCTATACCGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGAc                                                                                                                                                                                                        <  1:113651/236‑1 (MQ=255)
                   aaCACGCTGTTTGAAAAATCAACTCAACACTATACCGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGAc                                                                                                                                                                                                        <  1:210762/236‑1 (MQ=255)
                                    aTCAACTCAACACTATACCGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAg                                                                                                                                                                                                                        <  2:111200/203‑1 (MQ=255)
                                    aTCAACTCAACACTATACCGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAg                                                                                                                                                                                                                        >  1:111200/1‑203 (MQ=255)
                                       aactcaacACTATACCGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAag                                                                                                                                                                                    <  1:551108/236‑1 (MQ=255)
                                        actcaacACTATACCGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGa                                                                                                                                                                                                                           >  1:141939/1‑196 (MQ=255)
                                        actcaacACTATACCGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGa                                                                                                                                                                                                                           <  2:141939/196‑1 (MQ=255)
                                                   taCCGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATg                                                                                                                                                                        >  1:345133/1‑236 (MQ=255)
                                                   taCCGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACCAAAGCGAAAGCTaaag                                                                                                                                                                        >  1:337682/1‑234 (MQ=255)
                                                      cGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGaaa                                                                                                                                                                                                                         <  1:595485/184‑1 (MQ=255)
                                                      cGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGaaa                                                                                                                                                                                                                         >  2:595485/1‑184 (MQ=255)
                                                        aTAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGCGACGGATTATTATCTGACTAAAGAGAAAGATAATGTcaat                                                                                                                                                                   >  2:576453/1‑233 (MQ=255)
                                                        aTAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGGGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGGTaaatgccaa                                                                                                                                                                   >  2:357367/1‑229 (MQ=255)
                                                           gCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGAt                                                                                                                                                                            >  2:579481/1‑224 (MQ=255)
                                                           gCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGAt                                                                                                                                                                            <  1:579481/224‑1 (MQ=255)
                                                                                           tACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCgg                                                                                                                                >  1:429188/1‑236 (MQ=255)
                                                                                              aTGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTc                                                                                                                              <  2:429188/235‑1 (MQ=255)
                                                                                                     tCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCaa                                                                                                                      >  1:428193/1‑236 (MQ=255)
                                                                                                         ccTGCTGATTTGCGCCGGGATGGCGGTGCTGTGCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACGCCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAaaca                                                                                                                  <  1:64127/236‑1 (MQ=255)
                                                                                                              tgATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGCTTATTATCTGACTACAGAGACAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAaacaa                                                                                                                 >  1:610260/1‑232 (MQ=255)
                                                                                                                   tGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTg                                                                                                                                                                                                          >  2:330308/1‑138 (MQ=255)
                                                                                                                   tGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTg                                                                                                                                                                                                          <  1:330308/138‑1 (MQ=255)
                                                                                                                   tGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGGTTACCGTTGGCGGCCTTGGCTTCATCGGTCATGGGCATAACAACGGccctgc                                                                                                        >  2:45402/1‑232 (MQ=255)
                                                                                                                          ggATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGt                                                                                                      >  1:629983/1‑231 (MQ=255)
                                                                                                                          ggATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTAt                                                                                                  <  2:345133/235‑1 (MQ=255)
                                                                                                                                             gcgcACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGGCTg                                                                              >  1:366201/1‑236 (MQ=255)
                                                                                                                                             gcgcACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGGTTATCAGTCTCAAGCCGTGGTCTg                                                                              >  1:652246/1‑236 (MQ=255)
                                                                                                                                                             ttcttACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAgg                                                               >  2:212106/1‑235 (MQ=255)
                                                                                                                                                                  aCCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTAt                                                                                                  >  2:389809/1‑195 (MQ=255)
                                                                                                                                                                  aCCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTAt                                                                                                  <  1:389809/195‑1 (MQ=255)
                                                                                                                                                                       aagaGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACgtg                                                                        >  1:621574/1‑216 (MQ=255)
                                                                                                                                                                       aagaGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACgtg                                                                        <  2:621574/216‑1 (MQ=255)
                                                                                                                                                                               cgggtggTATTGAGGACCCACGCCCAGTTCCCTGCCGGTGCCACCATGGTTCAACCCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGa                                             <  2:337682/230‑1 (MQ=255)
                                                                                                                                                                                 gggggTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCa                                                                                                                                    >  2:166484/1‑146 (MQ=255)
                                                                                                                                                                                 gggggTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCa                                                                                                                                    <  1:166484/146‑1 (MQ=255)
                                                                                                                                                                                 gggggTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATc                                           >  2:311517/1‑235 (MQ=255)
                                                                                                                                                                                      tATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAgggg                                                                                                                         >  1:548059/1‑152 (MQ=255)
                                                                                                                                                                                      tATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAgggg                                                                                                                         <  2:548059/152‑1 (MQ=255)
                                                                                                                                                                                      tATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGcc                                                                                      >  2:611402/1‑187 (MQ=255)
                                                                                                                                                                                      tATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGcc                                                                                      <  1:611402/187‑1 (MQ=255)
                                                                                                                                                                                       aTTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGGTTATCAGTCTAAGCCGTGGTCTGAACGTTTCGGTGAGGAAAACTCGGTTACCTCGGTAAtttaact                                   >  2:387288/1‑231 (MQ=255)
                                                                                                                                                                                         ttATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCg                                   <  2:629983/235‑1 (MQ=255)
                                                                                                                                                                                              accaccGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGtctc                                                                                           <  2:6311/174‑1 (MQ=255)
                                                                                                                                                                                              accaccGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGtctc                                                                                           >  1:6311/1‑174 (MQ=255)
                                                                                                                                                                                               ccaccGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCGGGCAATGa                           <  1:387288/236‑1 (MQ=255)
                                                                                                                                                                                                caccGCGCAGTTACCTTCCGGTGCCACCATGGTTAACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCGGGCAATGa                           <  1:248316/236‑1 (MQ=255)
                                                                                                                                                                                                  ccGCGCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCGGGCAATGAtt                         <  1:491891/236‑1 (MQ=255)
                                                                                                                                                                                                   cgcgCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGat                                                                                                                                                                                                    >  2:150388/1‑64 (MQ=255)
                                                                                                                                                                                                   cgcgCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGat                                                                                                                                                                                                    <  1:150388/64‑1 (MQ=255)
                                                                                                                                                                                                   cgcgCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCg                                                                                     <  2:714532/175‑1 (MQ=255)
                                                                                                                                                                                                   cgcgCAGTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCg                                                                                     >  1:714532/1‑175 (MQ=255)
                                                                                                                                                                                                         tttACGTGCCGGTACAACCATGGTTCACACAACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAAAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCGGGCAATGATTGCGTTaa                  <  1:358383/235‑1 (MQ=255)
                                                                                                                                                                                                         tttACCGTCCGGTGCCACCATGGTTCACACCAGGAAAGTGCTGCAACAGGTTACGGATTATTATCTGACTCAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCGGGCAATGATTGCGTTaa                  <  2:428193/235‑1 (MQ=255)
                                                                                                                                                                                                         gTTACCTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCACCAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCGGGCAATGATTGCGTTaa                  <  2:349761/236‑1 (MQ=255)
                                                                                                                                                                                                              ctgccGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAACGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCGGGCAATGATTGCGTTAAGc                <  2:434542/230‑1 (MQ=255)
                                                                                                                                                                                                              cTTCCGGTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCGGGCAATGATTGCGTTAAGCAGt             <  1:212106/236‑1 (MQ=255)
                                                                                                                                                                                                                    gTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGa                                                                             <  2:531917/166‑1 (MQ=255)
                                                                                                                                                                                                                    gTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGa                                                                             >  1:531917/1‑166 (MQ=255)
                                                                                                                                                                                                                    gTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGGGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCGGGCAATGATTGCGTTAAGCAGGATCaaa       >  2:324334/1‑235 (MQ=255)
                                                                                                                                                                                                                    gTGCCACCATGGTTCACACCACGAAAGTGCTGCAACAGGAGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCGGGCAATGATTGCGTTAAGCAGTATCAAt       >  2:376499/1‑236 (MQ=255)
                                                                                                                                                                                                                          ccgatgTTCACACCCCGAAAGTTCTGCAACAGGGTACGGATTAGTATCGAACTAAAGAGAAAGAGAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCGGGCAATGATGGCGTTAAGCAGTATCAATAAAGc  <  2:610260/230‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  
GCCCTGTTCGCACGCTTCAACACGCTGTTTGAAAAATCAACTCAACACTATACCGATAGCACCCGCTCGCTGTTGCGTTGTACCGGTCGCTACATGGTGGTCTACCTGCTGATTTGCGCCGGGATGGCGGTGCTGTTCCTGCGCACGCCGACCTCTTTCTTACCAGAAGAGGATCAGGGGGTATTTATGACCACCGCGCAGTTACCTTCCGGTGCCACCATGGTTAACACCACGAAAGTGCTGCAACAGGTGACGGATTATTATCTGACTAAAGAGAAAGATAATGTCCAGTCGGTGTTTACCGTTGGCGGCTTTGGCTTCAGCGGTCAGGGGCAAAACAACGGCCTGGCGTTTATCAGTCTCAAGCCGTGGTCTGAACGTGTCGGTGAGGAAAACTCGGTTACCGCGATCATTCAGCGGGCAATGATTGCGTTAAGCAGTATCAATAAAGC  >  NC_000913/3661938‑3662389

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