Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,472,857 A→G intergenic (+117/‑28) bdcR → / → tabA transcriptional repressor for divergent bdcA/biofilm modulator regulated by toxins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,8570AG99.7% 1356.7 / NA 315intergenic (+117/‑28)bdcR/tabAtranscriptional repressor for divergent bdcA/biofilm modulator regulated by toxins
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (165/149);  minor base C (1/0);  total (166/149)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 320 total aligned reads displayed.

AGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATA  >  NC_000913/4472759‑4472954
                                                                                                  |                                                                                                 
gggCGGAAAAGCCCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCg                                                                                                 <  2:4756786/100‑1 (MQ=255)
aGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCg                                                                                                 >  1:1808999/1‑101 (MQ=255)
  gCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgtt                                                                                               >  1:6766687/1‑101 (MQ=255)
   cGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttg                                                                                              >  1:1213342/1‑101 (MQ=255)
   cGGAAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttg                                                                                              <  2:5515214/101‑1 (MQ=255)
    gggAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttgt                                                                                             <  2:1270901/101‑1 (MQ=255)
    gggAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttgt                                                                                             >  1:5031613/1‑101 (MQ=255)
     ggAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttgtt                                                                                            >  2:3115687/1‑101 (MQ=255)
       aaaGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGtttt                                                                                          >  1:6236517/1‑101 (MQ=255)
       aaaGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGtttt                                                                                          >  1:6893961/1‑101 (MQ=255)
         aGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaa                                                                                        <  2:6737266/101‑1 (MQ=255)
         aGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaa                                                                                        <  1:221610/101‑1 (MQ=255)
          gcgcCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaat                                                                                       <  1:1093655/101‑1 (MQ=255)
           cgcCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaatt                                                                                      >  2:1772930/1‑101 (MQ=255)
           cgcCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaatt                                                                                      >  2:1219161/1‑101 (MQ=255)
            gcCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGAATAACGATAGAAATGAAGGGTTGTTttaatta                                                                                     >  1:6412294/1‑101 (MQ=255)
            gcCTGCCGATAAGTTGTGATAAGACAATAACTCACGCGTTAAGGCTACCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTctaatta                                                                                     <  2:1035659/101‑1 (MQ=255)
             ccTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaattaa                                                                                    <  1:244825/101‑1 (MQ=255)
             ccTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaattaa                                                                                    >  1:2885110/1‑101 (MQ=255)
              cTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAg                                                                                   <  1:6522906/101‑1 (MQ=255)
               tGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAgg                                                                                  >  1:3295121/1‑101 (MQ=255)
                gCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGa                                                                                 <  1:2571464/101‑1 (MQ=255)
                gCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGa                                                                                 <  1:2715382/101‑1 (MQ=255)
                 ccGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAg                                                                                >  2:2309371/1‑101 (MQ=255)
                   gATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTg                                                                              >  1:6582557/1‑101 (MQ=255)
                   gATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTg                                                                              <  2:5634477/101‑1 (MQ=255)
                     tAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGa                                                                            <  1:4298073/101‑1 (MQ=255)
                        gTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaaa                                                                         <  2:1155842/101‑1 (MQ=255)
                        gTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaaa                                                                         >  1:6655708/1‑101 (MQ=255)
                        gTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaaa                                                                         >  1:5288841/1‑101 (MQ=255)
                         ttGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAg                                                                        >  1:3229911/1‑101 (MQ=255)
                           gtgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTgg                                                                             >  2:1985659/1‑94 (MQ=255)
                           gtgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgat                                                                      >  2:20899/1‑101 (MQ=255)
                            tgATAAGACATTAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatg                                                                     <  1:216165/101‑1 (MQ=255)
                                 aGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATACCGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcag                                                                <  1:4853212/101‑2 (MQ=255)
                                  gACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcatc                                                               <  1:6982648/101‑1 (MQ=255)
                                  gACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcatc                                                               <  1:6300946/101‑1 (MQ=255)
                                   aCAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCg                                                              <  1:2309371/101‑1 (MQ=255)
                                      ataatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaa                                                           <  1:5902708/101‑1 (MQ=255)
                                      ataatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaa                                                           <  1:5941029/101‑1 (MQ=255)
                                       taatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaaa                                                          >  2:5973921/1‑101 (MQ=255)
                                        aatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAat                                                         >  1:536020/1‑101 (MQ=255)
                                              cgcgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGCTAGAAATGAAGCGTTGTTTTCATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCAt                                                   <  2:3257963/101‑1 (MQ=255)
                                               gcgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATa                                                  >  2:1868949/1‑101 (MQ=255)
                                               gcgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATa                                                  >  1:3808120/1‑101 (MQ=255)
                                               gcgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATa                                                  >  1:3729058/1‑101 (MQ=255)
                                                cgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATaa                                                 <  1:176526/101‑1 (MQ=255)
                                                cgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATaa                                                 >  1:5177614/1‑101 (MQ=255)
                                                 gTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAAt                                                >  1:1992163/1‑101 (MQ=255)
                                                 gTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAAt                                                >  1:1952987/1‑101 (MQ=255)
                                                 gTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAAt                                                <  1:4239597/101‑1 (MQ=255)
                                                  ttAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATc                                               <  1:3451734/101‑1 (MQ=255)
                                                    aaGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCtt                                             >  1:1626260/1‑101 (MQ=255)
                                                     aGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTc                                            >  2:534089/1‑101 (MQ=255)
                                                       gCTAGCGGAATTGATTATCTTTTCGTATCACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAg                                          <  2:2192694/101‑1 (MQ=255)
                                                       gCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAg                                          >  1:5437218/1‑101 (MQ=255)
                                                        cTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGc                                         <  1:5674582/101‑1 (MQ=255)
                                                         tAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGcc                                        >  2:6446270/1‑101 (MQ=255)
                                                            cGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTg                                     >  1:817093/1‑101 (MQ=255)
                                                             ggAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTgg                                    >  2:503019/1‑101 (MQ=255)
                                                              gAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGc                                   <  2:2967639/101‑1 (MQ=255)
                                                              gAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGc                                   >  1:1022653/1‑101 (MQ=255)
                                                               aaTTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCt                                  <  1:5765099/101‑1 (MQ=255)
                                                                aTTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTa                                 <  1:3024496/101‑1 (MQ=255)
                                                                  tGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAcc                               >  2:5387505/1‑101 (MQ=255)
                                                                  tGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAcc                               >  1:5108695/1‑101 (MQ=255)
                                                                     ttATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCa                            >  1:3929141/1‑101 (MQ=255)
                                                                     ttATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCa                            <  2:6181034/101‑1 (MQ=255)
                                                                      tATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAg                           >  1:3186381/1‑101 (MQ=255)
                                                                       aTCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAgg                          <  2:2486830/101‑1 (MQ=255)
                                                                       aTCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAgg                          >  1:1418763/1‑101 (MQ=255)
                                                                            ttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTa                     >  1:475563/1‑101 (MQ=255)
                                                                             tCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAATGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTAc                    <  1:4977660/101‑1 (MQ=255)
                                                                             tCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTAc                    <  2:4659666/101‑1 (MQ=255)
                                                                              cGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACg                   >  2:2397564/1‑101 (MQ=255)
                                                                                  taACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAg               <  2:307892/101‑1 (MQ=255)
                                                                                   aaCGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAgg              <  1:4094016/101‑1 (MQ=255)
                                                                                   aaCGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAgg              <  1:6740001/101‑1 (MQ=255)
                                                                                    aCGATAGAAATGAAGCGTTGTTTTAATTAATTAGTTGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGc             >  2:6926327/1‑101 (MQ=255)
                                                                                     cGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCg            <  2:4932544/101‑1 (MQ=255)
                                                                                     cGATAGAAATGAAGCGTTGTTTTAATTAAGGAGGGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCg            >  2:4220994/1‑101 (MQ=255)
                                                                                      gATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGa           >  1:3135850/1‑101 (MQ=255)
                                                                                      gATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGa           >  1:6153844/1‑101 (MQ=255)
                                                                                       aTAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGAt          >  1:2116857/1‑101 (MQ=255)
                                                                                        tAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGAtt         >  2:4085537/1‑101 (MQ=255)
                                                                                        tAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGAtt         <  2:4619663/101‑1 (MQ=255)
                                                                                        tAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGAtt         >  2:3754969/1‑101 (MQ=255)
                                                                                        tAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGAtt         <  2:1626260/101‑1 (MQ=255)
                                                                                        tAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGAtt         >  2:3535776/1‑101 (MQ=255)
                                                                                        tAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATATTAAGGCGTGGTTACCCCAGGAGTTACGCCAGGCGAtt         <  2:6893961/101‑1 (MQ=255)
                                                                                          gAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGa       >  1:4530568/1‑101 (MQ=255)
                                                                                          gAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGa       >  2:2732758/1‑101 (MQ=255)
                                                                                           aaaTGAAGCGTTGTTTTAAGTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAg      <  1:2036525/101‑1 (MQ=255)
                                                                                             aTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCa    >  2:6350813/1‑101 (MQ=255)
                                                                                             aTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCa    >  2:3879855/1‑101 (MQ=255)
                                                                                              tGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCat   <  2:5216630/101‑1 (MQ=255)
                                                                                              tGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCat   >  2:1881526/1‑101 (MQ=255)
                                                                                              tGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCat   >  2:761951/1‑101 (MQ=255)
                                                                                              tGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCat   <  2:851225/101‑1 (MQ=255)
                                                                                               gAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCata  <  2:6306890/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
AGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATA  >  NC_000913/4472759‑4472954

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: