Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 654,566 G→T 11.9% intergenic (+24/+17) dpiA → / ← dcuC DNA‑binding transcriptional dual regulator DpiA/anaerobic C4‑dicarboxylate transporter DcuC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913654,5660GT11.9% 283.1 / 21.5 92intergenic (+24/+17)dpiA/dcuCDNA‑binding transcriptional dual regulator DpiA/anaerobic C4‑dicarboxylate transporter DcuC
Reads supporting (aligned to +/- strand):  ref base G (42/39);  new base T (7/4);  total (49/43)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

GCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCAT  >  NC_000913/654417‑654711
                                                                                                                                                     |                                                                                                                                                  
gCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTg                                                                                                                                                    >  2:513665/1‑150 (MQ=255)
     gACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTcc                                                                                                                                               >  1:854403/1‑150 (MQ=255)
       ccATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCgg                                                                                                                                             <  2:714246/150‑1 (MQ=255)
       ccATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCgg                                                                                                                                             <  2:664548/150‑1 (MQ=255)
        cATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                                                                            <  2:1343475/150‑1 (MQ=255)
         aTCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCa                                                                                                                                           <  2:1109975/150‑1 (MQ=255)
            aGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTa                                                                                                                                        >  1:1467476/1‑150 (MQ=255)
            aGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGATGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTa                                                                                                                                        >  2:1561752/1‑150 (MQ=255)
             gCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATGTACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTAc                                                                                                                                       <  1:1072800/150‑1 (MQ=255)
             gCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTAc                                                                                                                                       >  1:1275638/1‑150 (MQ=255)
             gCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTAc                                                                                                                                       >  2:1392005/1‑150 (MQ=255)
             gCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTAc                                                                                                                                       >  1:987335/1‑150 (MQ=255)
               cGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGCTTTatca                                                                                                                                     <  2:461524/150‑4 (MQ=255)
                gCACCACTGCCAGGCGTTATCTTCAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTCGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAt                                                                                                                                    >  2:277054/1‑150 (MQ=255)
                 caccacTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGCGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAtt                                                                                                                                   >  2:245964/1‑150 (MQ=255)
                 caccacTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGNNNNCAGACCACAACGCATATACCACAGTGGGCGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAtt                                                                                                                                   >  2:251976/1‑150 (MQ=255)
                  accacTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGCTTTatcatgt                                                                                                                                  <  1:1473475/150‑7 (MQ=255)
                  accacTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                                                                                  <  1:1402427/150‑1 (MQ=255)
                       tGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTg                                                                                                                             >  2:815642/1‑150 (MQ=255)
                        gCCAGGCGTTATCTTGAATAATGCGCCAGCCGCCATCTGATTATTGCTGAAATAGTTCACTGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTAGTATTTACTTGc                                                                                                                            >  1:1583227/1‑150 (MQ=255)
                         ccAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGcc                                                                                                                           <  1:527291/150‑1 (MQ=255)
                            ggCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGATCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgt                                                                                                                        <  2:553491/150‑1 (MQ=255)
                                     ttGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgc                                                                                                               <  2:368048/150‑1 (MQ=255)
                                          tatTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTa                                                                                                                                     >  1:307288/1‑123 (MQ=255)
                                          tatTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTa                                                                                                                                     <  2:307288/123‑1 (MQ=255)
                                              gcgcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTcc                                                                                                      <  1:223443/150‑1 (MQ=255)
                                               cgcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCt                                                                                                     >  2:420874/1‑150 (MQ=255)
                                               cgcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCt                                                                                                     >  2:414925/1‑150 (MQ=255)
                                                gcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTg                                                                                                    <  1:134905/150‑1 (MQ=255)
                                                   aGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCa                                                                                                 <  1:462924/150‑1 (MQ=255)
                                                   aGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCa                                                                                                 <  2:524833/150‑1 (MQ=255)
                                                     ccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAcc                                                                                               >  2:503692/1‑150 (MQ=255)
                                                     ccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAcc                                                                                               <  2:322448/150‑1 (MQ=255)
                                                        ccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATc                                                                                            <  2:894965/150‑1 (MQ=255)
                                                        ccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATc                                                                                            <  2:1046246/150‑1 (MQ=255)
                                                         cATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCa                                                                                           >  2:267507/1‑150 (MQ=255)
                                                         cATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCa                                                                                           >  2:1002022/1‑150 (MQ=255)
                                                         cATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCa                                                                                           <  2:139145/150‑1 (MQ=255)
                                                          aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAg                                                                                          >  1:1327814/1‑150 (MQ=255)
                                                           tCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGAACAATCAGc                                                                                         >  1:391103/1‑150 (MQ=255)
                                                             tGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctc                                                                                       >  1:1371805/1‑150 (MQ=255)
                                                                  attGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAg                                                                                  <  1:554956/150‑1 (MQ=255)
                                                                   ttGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGGATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGc                                                                                 >  1:567299/1‑150 (MQ=255)
                                                                       tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACg                                                                             <  2:579941/150‑1 (MQ=255)
                                                                         aaaTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGAt                                                                           <  2:1467476/150‑1 (MQ=255)
                                                                              gTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACaa                                                                      >  1:620303/1‑150 (MQ=255)
                                                                               ttCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAAt                                                                     >  1:1348132/1‑150 (MQ=255)
                                                                                tCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc                                                                    >  1:1043404/1‑150 (MQ=255)
                                                                                   cGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacc                                                                 <  2:492370/150‑1 (MQ=255)
                                                                                   cGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacc                                                                 >  1:1468822/1‑150 (MQ=255)
                                                                                   cGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacc                                                                 >  1:530300/1‑150 (MQ=255)
                                                                                         aGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaacca                                                           >  1:591198/1‑150 (MQ=255)
                                                                                           ttGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaac                                                         <  2:1058721/150‑1 (MQ=255)
                                                                                            tGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaaca                                                        <  1:768545/150‑1 (MQ=255)
                                                                                               cAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGc                                                     <  2:1275638/150‑1 (MQ=255)
                                                                                                   ccACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCg                                                 >  1:1019130/1‑150 (MQ=255)
                                                                                                    cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGGATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCgg                                                <  2:567299/150‑1 (MQ=255)
                                                                                                    cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAAt                                                                     >  2:113576/1‑129 (MQ=255)
                                                                                                    cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAAt                                                                     <  1:113576/129‑1 (MQ=255)
                                                                                                    cacaACGCATATACCACAGTGGGAGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCgg                                                <  2:744914/150‑1 (MQ=255)
                                                                                                         cGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACAATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCg                                           >  1:85437/1‑150 (MQ=255)
                                                                                                            atatACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAgg                                        <  2:1234818/150‑1 (MQ=255)
                                                                                                             tataCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAAc                                                                              >  1:203056/1‑111 (MQ=255)
                                                                                                             tataCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAAc                                                                              <  2:203056/111‑1 (MQ=255)
                                                                                                              ataCCAAACTGGTTGACATGATAAAGCCGGAGGAAACTTGCCTAAGGCATTACTATTTACTTGACTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTATAGCAACGATAACAATCACCAAACCAACAAGAACCGGTACCGAGGTg                                      >  2:1216961/1‑150 (MQ=255)
                                                                                                               taCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGc                                     >  2:848192/1‑150 (MQ=255)
                                                                                                                aCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCg                                    <  2:854403/150‑1 (MQ=255)
                                                                                                                      gggggTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTAc                              <  2:836780/148‑1 (MQ=255)
                                                                                                                      gTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTAc                              >  2:998067/1‑150 (MQ=255)
                                                                                                                      gTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTAc                              >  1:350583/1‑150 (MQ=255)
                                                                                                                        gggTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGa                            <  1:503692/150‑1 (MQ=255)
                                                                                                                         ggTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGAc                           <  1:267507/150‑1 (MQ=255)
                                                                                                                         ggTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGAc                           <  2:1043404/150‑1 (MQ=255)
                                                                                                                         ggTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGAc                           >  1:754624/1‑150 (MQ=255)
                                                                                                                          gTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACt                          <  2:85912/150‑1 (MQ=255)
                                                                                                                             aaatagtAATGCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCa                       >  2:1473136/7‑150 (MQ=255)
                                                                                                                             aCATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCa                       <  1:1561752/150‑1 (MQ=255)
                                                                                                                              cATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCaa                      >  1:600352/1‑150 (MQ=255)
                                                                                                                                tagtAATGCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAAc                    >  2:386395/4‑150 (MQ=255)
                                                                                                                                 gATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACg                   >  2:848866/1‑150 (MQ=255)
                                                                                                                                  gtAATGCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACgg                  >  1:1614313/2‑150 (MQ=255)
                                                                                                                                  aTAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAACCgg                  >  2:56862/1‑150 (MQ=255)
                                                                                                                                    aatgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCg                >  2:912751/4‑150 (MQ=255)
                                                                                                                                     atgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCga              >  1:1169734/3‑150 (MQ=255)
                                                                                                                                      tgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCCCCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCgag             >  1:1573634/2‑150 (MQ=255)
                                                                                                                                       gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCgaga            >  2:1450481/1‑150 (MQ=255)
                                                                                                                                         cGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATc          >  1:1601304/1‑150 (MQ=255)
                                                                                                                                         cGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATc          <  1:719204/150‑1 (MQ=255)
                                                                                                                                           gaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCtt        >  2:1121681/1‑150 (MQ=255)
                                                                                                                                             ggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCg      <  1:822736/150‑1 (MQ=255)
                                                                                                                                               aaaCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGcc    <  2:530300/150‑1 (MQ=255)
                                                                                                                                                 aCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCAt  >  1:819371/1‑150 (MQ=255)
                                                                                                                                                  cTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAGACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAAGATCtcgccgtt  <  2:1377692/150‑8 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCAT  >  NC_000913/654417‑654711

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