Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,663,555 G→T 8.9% intergenic (+28/+335) mdtF → / ← gadW anaerobic multidrug efflux transporter, ArcA‑regulated/transcriptional activator of gadA and gadBC; repressor of gadX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,663,5550GT 8.9% 273.1 / 16.1 102intergenic (+28/+335)mdtF/gadWanaerobic multidrug efflux transporter, ArcA‑regulated/transcriptional activator of gadA and gadBC; repressor of gadX
Reads supporting (aligned to +/- strand):  ref base G (48/45);  new base T (3/6);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01

Only 100 of 104 total aligned reads displayed.

AAAACGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACAT  >  NC_000913/3663406‑3663703
                                                                                                                                                     |                                                                                                                                                    
aaaaCGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTGCTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATg                                                                                                                                                      >  1:1867303/1‑150 (MQ=255)
aaaaCGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATg                                                                                                                                                      >  2:952730/1‑150 (MQ=255)
    cgcgGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAgg                                                                                                                                                  <  1:472035/150‑1 (MQ=255)
         taggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAg                                                                                                                                             >  1:1332451/1‑150 (MQ=255)
           ggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGcc                                                                                                                                           >  1:1379315/1‑150 (MQ=255)
            gtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCt                                                                                                                                          <  2:1495250/150‑1 (MQ=255)
             taCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCtt                                                                                                                                         <  1:542189/150‑1 (MQ=255)
               ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTTAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAc                                                                                                                                       >  1:439107/1‑150 (MQ=255)
               ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAc                                                                                                                                       >  1:695261/1‑150 (MQ=255)
               ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAc                                                                                                                                       <  1:1802416/150‑1 (MQ=255)
                 ggCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTActc                                                                                                                                     <  1:619835/150‑1 (MQ=255)
                 ggCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTActc                                                                                                                                     <  2:125071/150‑1 (MQ=255)
                  gCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGGCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTActct                                                                                                                                    <  1:1015379/150‑1 (MQ=39)
                   ggttaTGGGCTGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTActctc                                                                                                                                   <  1:1688737/146‑2 (MQ=38)
                      gATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTcc                                                                                                                                <  2:851118/150‑1 (MQ=255)
                      gATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTcc                                                                                                                                <  1:925374/150‑1 (MQ=255)
                      gATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGGCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTcc                                                                                                                                <  2:923163/150‑1 (MQ=255)
                      gATGGGCGGGATGTGTGCCGCAACAGTGCTGGCAATTTAGTTCGTTGCGGTCTTTTTCGTTGTAGGGGAACATCTCTTTGCCCTCTTTAAAAAAGCGTAAGGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTcc                                                                                                                                <  2:267554/150‑1 (MQ=255)
                         gggCGTGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTGCGTTCTTTTTGGTTGTAGTGGAACATCTCTTTGCCCGCTTTTAAAAAGCGTATCGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTActctcatttac                                                                                                                             <  2:687789/150‑8 (MQ=25)
                              ggATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAAcc                                                                                                                        >  2:502199/1‑150 (MQ=255)
                               gATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCa                                                                                                                       <  2:596960/150‑1 (MQ=255)
                                  gTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTc                                                                                                                    >  2:1624084/1‑150 (MQ=255)
                                  gTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTc                                                                                                                    >  2:1558028/1‑150 (MQ=255)
                                      gCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTACCCATTcacca                                                                                                                >  2:704990/1‑150 (MQ=255)
                                        cGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCAGTcaccaca                                                                                                              >  1:1433662/1‑149 (MQ=255)
                                         gCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTc                                                                                                             >  1:731078/1‑150 (MQ=255)
                                                   tGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGc                                                                                                   >  1:1753105/1‑150 (MQ=255)
                                                     gCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTa                                                                                                 <  1:408667/150‑1 (MQ=255)
                                                     gCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTa                                                                                                 <  2:1744994/150‑1 (MQ=255)
                                                      cAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTaa                                                                                                >  1:1535724/1‑150 (MQ=255)
                                                       aaTTTGCTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAGCAACCAGCTAAc                                                                                               <  2:1867303/150‑1 (MQ=255)
                                                            aCTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCa                                                                                          <  1:427946/150‑1 (MQ=255)
                                                             cTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAg                                                                                         >  2:620539/1‑150 (MQ=255)
                                                                  ttCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATa                                                                                    <  2:1338217/150‑1 (MQ=255)
                                                                     cGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCa                                                                                 >  1:1021814/1‑150 (MQ=255)
                                                                         cTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTACCTGCTAACCATTCACAACACAACAACCAGCTAACCGTCAGTAATACCTGGCt                                                                             >  2:535134/1‑150 (MQ=255)
                                                                          tttttCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCtt                                                                            >  2:826267/1‑150 (MQ=255)
                                                                          tttttCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCtt                                                                            <  2:465299/150‑1 (MQ=255)
                                                                           ttttCGTTGTAGTGGAACATCTCTTTGCCCGCTTTATAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTa                                                                           <  1:252096/150‑1 (MQ=255)
                                                                           ttttCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAATGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCGTACTCTGTCCTGCTAACCCTTCACCACTCTCCACCCAGCTAACCGTCAGGAATTCCAGGCTTa                                                                           >  1:77318/1‑150 (MQ=255)
                                                                           ttttCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAAGATGATGGTTCAGCCTTACTCTTTCATGCTTACCATTGACCACTAATCAACCAGCTTAGCGTCAGTCATACCTGGCttt                                                                           >  1:956305/1‑149 (MQ=255)
                                                                           ttttCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTa                                                                           >  1:1290955/1‑150 (MQ=255)
                                                                            tttCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAg                                                                          >  1:847031/1‑150 (MQ=255)
                                                                            tttCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGAGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAg                                                                          <  2:848919/150‑1 (MQ=255)
                                                                             ttCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGt                                                                         <  2:1521409/150‑1 (MQ=255)
                                                                               cGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTac                                                                       <  1:1151880/150‑1 (MQ=255)
                                                                                gTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATGCCAGGCTTAGTaca                                                                      >  1:1757998/1‑150 (MQ=255)
                                                                                gTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTacc                                                                      >  2:753076/1‑149 (MQ=255)
                                                                                gTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTaca                                                                      >  1:152303/1‑150 (MQ=255)
                                                                                 ttGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCCCCACTCAACAACCAGCTAACCGTCAGGAAAACCAGGCTTAGTacac                                                                     >  2:729310/1‑150 (MQ=255)
                                                                                  tGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTAcaca                                                                    >  2:630391/1‑150 (MQ=255)
                                                                                   gTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAg                                                                   <  1:507536/150‑1 (MQ=255)
                                                                                   gTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAg                                                                   >  2:667656/1‑150 (MQ=255)
                                                                                    tAGTGGAACATCTCTTTGCCCGCTTTTAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGc                                                                  <  2:439107/150‑1 (MQ=255)
                                                                                    tAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGc                                                                  <  2:238914/150‑1 (MQ=255)
                                                                                    tAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGc                                                                  <  2:1194446/150‑1 (MQ=255)
                                                                                    tAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGACTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTCCACAGc                                                                  >  1:1062170/1‑150 (MQ=255)
                                                                                      gTGGAACATCTCTTTGCCCGCTTTATAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCAc                                                                >  1:1225568/1‑150 (MQ=255)
                                                                                        ggAACATCTCTTTTCCCGCATTAAAAAAGCGTAACGTTTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTg                                                              <  1:1600670/150‑1 (MQ=255)
                                                                                          aaCATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGaa                                                            >  2:337367/1‑150 (MQ=255)
                                                                                           aCATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGaaa                                                           >  2:1230869/1‑150 (MQ=255)
                                                                                             aTCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGt                                                         <  1:1077143/150‑1 (MQ=255)
                                                                                               ctctTTGCCCGCTTTAAAAAAGCCTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAAACATTCACAACTCAAAAACCAGATAACCATCAGGAATACAAGGCTTAGTACAAAGCACTGAAAGTAg                                                       >  1:1022326/1‑150 (MQ=255)
                                                                                                  tttGCCAGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGaaa                                                    <  1:132871/150‑1 (MQ=255)
                                                                                                   ttGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAAc                                                   <  2:1283863/150‑1 (MQ=255)
                                                                                                   ttGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAAc                                                   <  1:1097598/150‑1 (MQ=255)
                                                                                                    tGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCCACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAAcc                                                  >  1:1416427/1‑150 (MQ=255)
                                                                                                    tGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAAcc                                                  <  2:97538/150‑1 (MQ=255)
                                                                                                        cGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCCACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTc                                              >  1:947488/1‑150 (MQ=255)
                                                                                                        cGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTc                                              >  1:1472867/1‑150 (MQ=255)
                                                                                                         gCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTcc                                             >  1:1867920/1‑150 (MQ=255)
                                                                                                          cTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCa                                            <  1:809404/150‑1 (MQ=255)
                                                                                                            ttAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGc                                          <  1:502199/150‑1 (MQ=255)
                                                                                                            ttAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGc                                          <  2:1760432/150‑1 (MQ=255)
                                                                                                              aaaaaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTACCTGCTAACCATTAACCACTCAACAACCAGCTAACCGTCAGGAATACCACGCTTAGTACACAGCACTAATAGTAGAAACCAATCCAGCCa                                        >  2:1581770/1‑150 (MQ=255)
                                                                                                               aaaaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACCCAGCACTGAAAGTAGAAACCACTCCAGCCAt                                       >  2:1273805/1‑150 (MQ=255)
                                                                                                                   aGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGt                                   >  1:1750938/1‑150 (MQ=255)
                                                                                                                    gCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCTGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtc                                  <  2:271196/150‑1 (MQ=255)
                                                                                                                    gCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtc                                  <  1:668481/150‑1 (MQ=255)
                                                                                                                    ccGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtc                                  <  2:302808/149‑1 (MQ=255)
                                                                                                                      gTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCCCCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtcat                                >  1:1001012/1‑150 (MQ=255)
                                                                                                                      gTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAACCATTCGtcac                                <  2:1535724/150‑2 (MQ=255)
                                                                                                                       tAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCCCCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGCACCCACTCCAGCCATTCGtcatc                               >  1:1020960/1‑150 (MQ=255)
                                                                                                                       tAACGGGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtcatc                               <  1:945390/150‑1 (MQ=255)
                                                                                                                         aCGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCACCAACCAGCTACCCGTCAGGAATACCAGGCTTAGTACACAGCCCTGAAAGTAGAAACCACTCCAGCCTTTCGTCATCaa                             >  1:686220/1‑150 (MQ=255)
                                                                                                                          cGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAAt                            >  1:316567/1‑150 (MQ=255)
                                                                                                                             gtAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTg                         <  2:1332451/150‑1 (MQ=255)
                                                                                                                              caggaaagagTAAGGCTGAACCTTCATGTTCAACCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGa                        >  1:1503185/7‑150 (MQ=38)
                                                                                                                               aaaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaa                       <  1:707786/150‑1 (MQ=255)
                                                                                                                               aaaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaa                       <  1:661989/150‑1 (MQ=255)
                                                                                                                               aaaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaa                       <  1:261755/150‑1 (MQ=255)
                                                                                                                                ggaaagagTAAGGCTGAACCTTCATGTTCAACCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAAATTGaac                      >  2:336992/5‑150 (MQ=38)
                                                                                                                                ggaaagagTAAGGCTGAACCTTCATGATCAACCTTAATCTTTCCTGCTAACAATTCACCACTAAACAACCCGCTAACCGTCAGGAATACAAAGCTTAGTGCCGAGCACTGAAAGGAGAAACCACTCCAGCCATTCGTCATAAATTTGaaa                      >  1:1470959/5‑149 (MQ=255)
                                                                                                                                aaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaac                      <  1:1399418/150‑1 (MQ=255)
                                                                                                                                aaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAATCTTTACTGCTAACAATTCACCACTCAACAACCAGCTAACCGTGAGGAATACCAGGCTAAGTACACAGCACTGAAAGAAGAAAACCCTCCAGCCATTAGTCATCAATTTGata                      >  1:1776673/1‑148 (MQ=255)
                                                                                                                                 aTGAGAGTAAGGTTGAACATGAAGGTGCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaaca                     <  1:1305209/150‑1 (MQ=255)
                                                                                                                                   gagaGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTATAAACCACTCCAGCCATTCGTCAGCAATTTCaacacc                   >  1:1230844/1‑148 (MQ=255)
                                                                                                                                     gagTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAacac                 <  2:1702338/150‑1 (MQ=255)
                                                                                                                                                    tGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACGCACAt  >  2:810366/1‑150 (MQ=255)
                                                                                                                                                    tGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACAt  >  2:1524998/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAACGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACAT  >  NC_000913/3663406‑3663703

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