Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,662,481 C→A 5.1% H690N (CAC→AAC)  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,4810CA 5.1% 270.7 / 2.8 78H690N (CAC→AAC) mdtFanaerobic multidrug efflux transporter, ArcA‑regulated
Reads supporting (aligned to +/- strand):  ref base C (35/39);  new base A (2/2);  total (37/41)
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.94e-01

TCGGTTACCGCGATCATTCAGCGGGCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCAAAGTGAACGTCAACGCTGCGAAAGCTGAAGCGATG  >  NC_000913/3662334‑3662624
                                                                                                                                                   |                                                                                                                                               
tCGGTTACCGCGATCATTCAGCGGGCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCa                                                                                                                                                >  3:160865/1‑149 (MQ=255)
      aCCGCGATCATTCAGCGGGCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGaaa                                                                                                                                           >  1:416792/1‑148 (MQ=255)
          cgATCATTCAGCGGGCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGTTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTa                                                                                                                                      <  6:212169/149‑1 (MQ=255)
           gATCATTCAGCGGGCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTaa                                                                                                                                     >  6:181745/1‑149 (MQ=255)
            aTCATTCAGCGGGCAATGATTGCGTTAAGCAGTATCAATAAAGCTGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAAc                                                                                                                                    >  2:252093/1‑149 (MQ=255)
            aTCATTCAGCGGGCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAAc                                                                                                                                    <  7:151841/149‑1 (MQ=255)
              cATTCAGCGGGCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAAcc                                                                                                                                   >  6:36821/1‑148 (MQ=255)
               aTTCAGCGGGCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTTGACAACGGTAACCTGGGGCACGAAAAACTAAccc                                                                                                                                  <  1:210732/148‑1 (MQ=255)
               aTTCAGCGGGCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAAccc                                                                                                                                  <  7:113381/148‑1 (MQ=255)
                     cGGGCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGgcgc                                                                                                                            >  8:230767/1‑148 (MQ=255)
                     cGGGCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGcgcg                                                                                                                           <  4:129153/149‑1 (MQ=255)
                      gggCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGgc                                                                                                                              <  7:237789/145‑1 (MQ=255)
                      gggCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGgc                                                                                                                              >  8:237789/1‑145 (MQ=255)
                           aTGATTGCGTTAAGCAGTATCAATAAAGCAGTCGTCTTACCGTTAAACTTACCCGCGGTGGCTGAACTGGTTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGa                                                                                                                     >  3:249981/1‑149 (MQ=255)
                            tGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAg                                                                                                                    >  3:30330/1‑149 (MQ=255)
                            tGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAg                                                                                                                    <  4:240466/149‑1 (MQ=255)
                                         cAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTgg                                                                                                       <  5:181745/149‑1 (MQ=255)
                                           gTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTgg                                                                                                       >  3:332747/1‑147 (MQ=255)
                                            tATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAg                                                                                                    >  3:370121/1‑149 (MQ=255)
                                                aaTAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAgcgc                                                                                                 >  2:76851/1‑148 (MQ=255)
                                                aaTAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAgcgc                                                                                                 <  2:395324/148‑1 (MQ=255)
                                                 aTAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTaa                                                                                                                                     <  7:111290/111‑1 (MQ=255)
                                                 aTAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTaa                                                                                                                                     >  8:111290/1‑111 (MQ=255)
                                                  tAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGTTTTTGATATGGAACTGCGGGTCAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCCGCGCAAt                                                                                              >  3:245965/1‑149 (MQ=255)
                                                    aaGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAgcgc                                                                                                 >  5:156788/1‑144 (MQ=255)
                                                    aaGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAgcgc                                                                                                 <  6:156788/144‑1 (MQ=255)
                                                      gCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCAcc                                                                                          >  2:66993/1‑149 (MQ=255)
                                                      gCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCAcc                                                                                          >  4:23060/1‑149 (MQ=255)
                                                            gtcTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAAt                                                                                      >  5:186079/1‑147 (MQ=255)
                                                            gtcTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATc                                                                                     <  7:230767/148‑1 (MQ=255)
                                                             tcTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAg                                                                                   <  1:252093/149‑1 (MQ=255)
                                                             tcTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAg                                                                                   <  4:336520/149‑1 (MQ=255)
                                                               ttCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCa                                                                                    <  1:185575/146‑1 (MQ=255)
                                                               ttCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCa                                                                                    >  2:185575/1‑146 (MQ=255)
                                                                 cccGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTc                                                                                <  8:74763/148‑1 (MQ=255)
                                                                 cccGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCa                                                                               >  1:401321/1‑149 (MQ=255)
                                                                  ccGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCAc                                                                              >  3:191883/1‑149 (MQ=255)
                                                                    gTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGgc                                                                                                                              <  8:274432/99‑1 (MQ=255)
                                                                    gTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGgc                                                                                                                              >  7:274432/1‑99 (MQ=255)
                                                                        aaCTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCgggg                                                                         <  1:369440/148‑1 (MQ=255)
                                                                        aaCTTACCCGCGGTGGCTGAACTGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGAACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTa                                                                       >  2:275247/1‑149 (MQ=255)
                                                                        aaCTTACCCGCGGTGGCTGAACTGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGAACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTa                                                                       >  4:48368/1‑149 (MQ=255)
                                                                             aCCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGcc                                                                   >  8:269659/1‑149 (MQ=255)
                                                                              cccGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGcc                                                                   <  4:327314/148‑1 (MQ=255)
                                                                               ccGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGccc                                                                  >  2:408833/1‑148 (MQ=255)
                                                                                 gcgGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGaa                                                               <  4:332747/149‑1 (MQ=255)
                                                                                  cgGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAAc                                                              >  8:71291/1‑149 (MQ=255)
                                                                                     tGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGc                                                           <  7:269659/149‑1 (MQ=255)
                                                                                           aaCTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTgg                                                       <  1:76851/147‑1 (MQ=255)
                                                                                            aCTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAg                                                    <  6:186079/149‑1 (MQ=255)
                                                                                             cTGGGTACCGCGTCATGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGa                                                   <  4:370121/149‑1 (MQ=255)
                                                                                                ggTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATa                                                 <  1:408833/148‑1 (MQ=255)
                                                                                                ggTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATa                                                 <  7:71291/148‑1 (MQ=255)
                                                                                                ggTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATa                                                 <  2:416792/148‑1 (MQ=255)
                                                                                                  tACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGc                                              <  3:258908/149‑1 (MQ=255)
                                                                                                  tACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGc                                              <  4:245965/149‑1 (MQ=255)
                                                                                                    ccGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGcc                                             >  7:38932/1‑148 (MQ=255)
                                                                                                      gcgTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGAt                                          >  5:54198/1‑149 (MQ=255)
                                                                                                      gcgTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGAt                                          >  5:302274/1‑149 (MQ=255)
                                                                                                           aGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTc                                      <  2:401321/148‑1 (MQ=255)
                                                                                                             gTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCaaa                                   <  5:36821/149‑1 (MQ=255)
                                                                                                                  gATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCAAAGTGaa                              >  4:361649/1‑149 (MQ=255)
                                                                                                                   atatGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCTCTGGCAGCGCATTCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCAAAGTGaa                              <  4:30330/148‑1 (MQ=255)
                                                                                                                   atatGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCAAAGTGAAc                             >  4:189943/1‑149 (MQ=255)
                                                                                                                       ggAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGa                                           >  3:143428/1‑131 (MQ=255)
                                                                                                                       ggAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGa                                           <  4:143428/131‑1 (MQ=255)
                                                                                                                        gAACTGCTGGACAACGGTAACCTGGGGAACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCAAAGTGAACGTCaa                        <  1:275247/149‑1 (MQ=255)
                                                                                                                        gAACTGCTGGACAACGGTAACCTGGGGAACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCAAAGTGAACGTCaa                        <  3:48368/149‑1 (MQ=255)
                                                                                                                         aaCTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGTCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATCTTCAAAGTGAACGTCAAc                       >  4:320575/1‑149 (MQ=255)
                                                                                                                              ctgGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCAAAGTGAACGTCAACGCTGc                  <  3:361649/149‑1 (MQ=255)
                                                                                                                                 cccAACGGTATCCTGCGGCACGAAAAACTAACCCAGGCCCGAAACGAGCTGTCATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCAAAGTGAACGTCAACTCTGCGa                <  4:249981/146‑1 (MQ=255)
                                                                                                                                    aaCGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCaaa                                   >  6:176515/1‑126 (MQ=255)
                                                                                                                                    aaCGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCaaa                                   <  5:176515/126‑1 (MQ=255)
                                                                                                                                       ggTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCAAAGTGAACGTCAACGCTGCGAAAGCTg          >  1:368329/1‑148 (MQ=255)
                                                                                                                                        gTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGGAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGATGATACGCCGATGTTCAAAGTGAACGTCAACGCTGCGAAAGCTGaa        <  4:160865/149‑1 (MQ=255)
                                                                                                                                             cTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCAAAGTGAACGTCAACGCTGCGAAAGCTGAAGCGAt   <  1:66993/149‑1 (MQ=255)
                                                                                                                                             cTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCAAAGTGAACGTCAACGCTGCGAAAGCTGAAGCGAt   >  7:408230/1‑149 (MQ=255)
                                                                                                                                                gggCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCAAAGTGAACGTCAACGCTGCGAAAGCTGAAGCGATg  <  4:191883/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TCGGTTACCGCGATCATTCAGCGGGCAATGATTGCGTTAAGCAGTATCAATAAAGCCGTCGTCTTCCCGTTCAACTTACCCGCGGTGGCTGAACTGGGTACCGCGTCAGGTTTTGATATGGAACTGCTGGACAACGGTAACCTGGGGCACGAAAAACTAACCCAGGCGCGAAACGAGCTGTTATCACTGGCAGCGCAATCACCGAATCAGGTCACCGGGGTACGCCCGAACGGCCTGGAAGATACGCCGATGTTCAAAGTGAACGTCAACGCTGCGAAAGCTGAAGCGATG  >  NC_000913/3662334‑3662624

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