Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,161,702 C→G A162G (GCA→GGA)  btuB → BtuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,7020CG99.6% 826.1 / NA 234A162G (GCA→GGA) btuBBtuB
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (113/120);  total (113/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

GGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATT  >  NC_000913/4161612‑4161791
                                                                                          |                                                                                         
ggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTcaggag                                                                                         >  2:3505068/1‑93 (MQ=255)
    gcgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCCGGAggggg                                                                                     >  1:2281347/1‑91 (MQ=255)
    gcgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTg                                                                                     <  1:4876656/93‑1 (MQ=255)
    gcgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGAGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTg                                                                                     >  2:1776584/1‑93 (MQ=255)
      gcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTggg                                                                                   >  2:1363824/1‑93 (MQ=255)
      gcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTggg                                                                                   <  1:4496598/93‑1 (MQ=255)
         ccGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGATGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                                                >  2:946469/1‑93 (MQ=255)
         ccGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                                                >  2:3349176/1‑93 (MQ=255)
          cGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAg                                                                               >  1:1361029/1‑93 (MQ=255)
           gCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGc                                                                              <  1:2131891/93‑1 (MQ=255)
            cTGTTTATGGTTCCGATGCAGTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCa                                                                             <  2:1475389/93‑1 (MQ=255)
             tGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCaa                                                                            <  2:2827209/93‑1 (MQ=255)
             tGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCaa                                                                            >  1:5580689/1‑93 (MQ=255)
              gTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCaaa                                                                           >  2:4195391/1‑92 (MQ=255)
              gTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAAt                                                                           >  1:461457/1‑93 (MQ=255)
              gTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAAt                                                                           <  1:6072545/93‑1 (MQ=255)
              gTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAAt                                                                           >  1:2856977/1‑93 (MQ=255)
              gTTTATGGTTCCGATGCAATAGGCGGGGTGGAGAATATCATCAAGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGCGGGGGAGCAAt                                                                           >  1:5091778/1‑93 (MQ=255)
               tttATGGTTCCGATGCAATAGGCGGGGTTGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATa                                                                          <  2:6215434/93‑1 (MQ=255)
               tttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATa                                                                          >  2:5198005/1‑93 (MQ=255)
               tttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATa                                                                          <  2:1109928/93‑1 (MQ=255)
                ttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                                         <  1:1483173/93‑1 (MQ=255)
                ttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                                         <  1:4980140/93‑1 (MQ=255)
                 tATGGTTCCGATGCAATAGGCGTGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCATGAGGGTGGGGAAGCAATAGt                                                                        <  1:2985082/93‑1 (MQ=255)
                   tGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTa                                                                      <  1:2486675/93‑1 (MQ=255)
                   tGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTa                                                                      >  1:531998/1‑93 (MQ=255)
                   tGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTa                                                                      >  1:2433488/1‑93 (MQ=255)
                    ggTTCCGATGCAATAGGCGGGTTGGTGAATATCATCACGACGCGCGATGAACCCGGACCGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTAt                                                                     <  1:56196/93‑1 (MQ=255)
                     gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATc                                                                    <  2:3115929/93‑1 (MQ=255)
                     gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATc                                                                    >  1:3424613/1‑93 (MQ=255)
                         cGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGaa                                                                <  2:1462727/93‑1 (MQ=255)
                            tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTa                                                             >  2:5208915/1‑93 (MQ=255)
                            tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTa                                                             >  2:2985684/1‑93 (MQ=255)
                             gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTat                                                            <  2:1170239/93‑1 (MQ=255)
                             gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTat                                                            <  1:1532479/93‑1 (MQ=255)
                              cAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatg                                                           >  1:1697879/1‑93 (MQ=255)
                               aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatga                                                          <  2:2100591/93‑1 (MQ=255)
                                aTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgat                                                         >  2:786224/1‑93 (MQ=255)
                                aTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgat                                                         >  2:4497848/1‑93 (MQ=255)
                                   ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                                      >  1:1237996/1‑93 (MQ=255)
                                   ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                                      <  2:3680202/93‑1 (MQ=255)
                                     cGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtctc                                                    <  1:5963264/93‑1 (MQ=255)
                                     cGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtctc                                                    <  2:2237121/93‑1 (MQ=255)
                                      ggggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTctct                                                   <  2:4090514/93‑1 (MQ=255)
                                      ggggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTctct                                                   <  1:3890360/93‑1 (MQ=255)
                                      ggggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTctct                                                   <  1:4119432/93‑1 (MQ=255)
                                       gggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                                  >  2:4610393/1‑93 (MQ=255)
                                       gggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                                  >  2:2541939/1‑93 (MQ=255)
                                        ggtggtGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAc                                                 >  2:5318920/1‑93 (MQ=255)
                                        ggtggtGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAc                                                 <  2:3011253/93‑1 (MQ=255)
                                         gtggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACg                                                >  1:5658892/1‑93 (MQ=255)
                                           ggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgca                                              >  1:5879541/1‑93 (MQ=255)
                                           ggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgca                                              >  2:3803887/1‑93 (MQ=255)
                                           ggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgca                                              <  2:10000/93‑1 (MQ=255)
                                            gtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGCGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcag                                             >  1:1967284/1‑93 (MQ=255)
                                                atatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaac                                         >  2:5713149/1‑93 (MQ=255)
                                                    catcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACt                                     >  1:3050994/1‑93 (MQ=255)
                                                     atcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                    >  1:812754/1‑93 (MQ=255)
                                                     atcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                    >  2:1623961/1‑93 (MQ=255)
                                                     atcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                    <  1:4464333/93‑1 (MQ=255)
                                                     atcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                    <  1:5554143/93‑1 (MQ=255)
                                                      tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgg                                   >  2:2762929/1‑93 (MQ=255)
                                                          gacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGa                               <  2:4633280/93‑1 (MQ=255)
                                                          gacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGa                               >  2:2400620/1‑93 (MQ=255)
                                                           acgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAt                              <  1:6330689/93‑1 (MQ=255)
                                                             gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa                            <  2:4413444/93‑1 (MQ=255)
                                                               gcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGa                          >  2:3775619/1‑93 (MQ=255)
                                                                cgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac                         >  1:2953554/1‑93 (MQ=255)
                                                                cgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac                         >  2:5464998/1‑93 (MQ=255)
                                                                cgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac                         <  2:5834694/93‑1 (MQ=255)
                                                                 gATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGaca                        >  2:6038926/1‑93 (MQ=255)
                                                                 gATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGaca                        <  1:5176086/93‑1 (MQ=255)
                                                                  aTGAACCCGGAGCGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGacac                       <  1:6212672/93‑1 (MQ=255)
                                                                  aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGacac                       >  2:6448415/1‑93 (MQ=255)
                                                                   tGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACg                      >  2:4891622/1‑93 (MQ=255)
                                                                   tGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACg                      >  2:3857882/1‑93 (MQ=255)
                                                                    gaccccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACgg                     <  2:1846548/90‑1 (MQ=255)
                                                                    gAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACgg                     >  2:3913428/1‑93 (MQ=255)
                                                                    gAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACgg                     >  1:1030409/1‑93 (MQ=255)
                                                                    gAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCGACGCAGCAACAACTGGGGGATAAGACACgg                     <  1:3173465/93‑1 (MQ=255)
                                                                     aaCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACggg                    >  1:2171073/1‑93 (MQ=255)
                                                                      aCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGt                   >  1:6281882/1‑93 (MQ=255)
                                                                       cccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTa                  <  1:854829/93‑1 (MQ=255)
                                                                       cccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTa                  <  1:5473864/93‑1 (MQ=255)
                                                                         cggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAAc                >  1:490952/1‑93 (MQ=255)
                                                                         cggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAAc                <  2:222858/93‑1 (MQ=255)
                                                                           gaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGc              >  2:125751/1‑93 (MQ=255)
                                                                           gaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGc              >  1:431092/1‑93 (MQ=255)
                                                                           gaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGc              <  1:3603912/93‑1 (MQ=255)
                                                                            aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCt             <  1:5803034/93‑1 (MQ=255)
                                                                             acggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTg            >  1:4886172/1‑93 (MQ=255)
                                                                             acggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTg            <  2:4161139/93‑1 (MQ=255)
                                                                                gaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTg         <  2:4828217/93‑1 (MQ=255)
                                                                                gaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTg         <  1:4836881/93‑1 (MQ=255)
                                                                                   aTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGc      >  2:590890/1‑93 (MQ=255)
                                                                                   aTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGc      >  1:5050463/1‑93 (MQ=255)
                                                                                    tttCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCg     <  2:2632105/93‑1 (MQ=255)
                                                                                      tCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAt   <  1:1598424/93‑1 (MQ=255)
                                                                                       caggagGGTGGGGGAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAtt  >  1:3858430/1‑93 (MQ=255)
                                                                                       caggagGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAtt  <  2:5537003/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATT  >  NC_000913/4161612‑4161791

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