Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,673,783 G→A 15.6% intergenic (+15/+33) yphA → / ← yphB DoxX family inner membrane protein/mutarotase superfamily protein, YphB family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,673,7830GA15.6% 448.8 / 54.8 177intergenic (+15/+33)yphA/yphBDoxX family inner membrane protein/mutarotase superfamily protein, YphB family
Reads supporting (aligned to +/- strand):  ref base G (78/70);  new base A (17/11);  total (96/81)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

TTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGA  >  NC_000913/2673684‑2673882
                                                                                                   |                                                                                                   
 ttCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgc                                                                                                    <  2:797255/100‑1 (MQ=255)
 ttCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgc                                                                                                    <  2:3377265/100‑1 (MQ=255)
  tCTGGAAGAATGTCAGTATCGCTGGCGCTTTCTTGCTATTGGTAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcg                                                                                                   <  2:1667456/100‑1 (MQ=255)
   cTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcgg                                                                                                  <  2:2602000/100‑1 (MQ=255)
   cTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcgg                                                                                                  >  1:829930/1‑100 (MQ=255)
   cTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCCATTGGCACTTCCCTGGCCGGGGGCAATTTCTCTCGATCCTCGTTAGGAAAGCTGCCGGATgcgg                                                                                                  >  2:2197125/1‑100 (MQ=255)
    tGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcggc                                                                                                 >  1:603840/1‑100 (MQ=255)
     ggAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcggcg                                                                                                <  2:512917/100‑1 (MQ=255)
     ggAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcggcg                                                                                                >  2:3035240/1‑100 (MQ=255)
      gaagaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGt                                                                                               <  1:2837084/100‑1 (MQ=255)
      gaagaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGt                                                                                               <  1:416677/100‑1 (MQ=255)
       aagaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTg                                                                                              <  1:2381098/100‑1 (MQ=255)
        agaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGa                                                                                             <  2:1424648/100‑1 (MQ=255)
         gaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGaa                                                                                            <  2:982108/100‑1 (MQ=255)
          aaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAAc                                                                                           >  1:528839/1‑100 (MQ=255)
          aaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAAc                                                                                           <  1:1071293/100‑1 (MQ=255)
             gTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                        <  2:1082470/100‑1 (MQ=255)
             gTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                        <  1:3249468/100‑1 (MQ=255)
             gTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                        <  1:2308115/100‑1 (MQ=255)
             cgatctATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                        <  1:1744739/95‑1 (MQ=255)
              tCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCt                                                                                       >  1:668111/1‑100 (MQ=255)
                        tGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCan                                                                             <  1:101675/100‑2 (MQ=25)
                         ggCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCAtt                                                                            >  1:99438/1‑100 (MQ=255)
                         ggCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAtc                                                                            <  1:3436056/100‑2 (MQ=25)
                         ggCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                             <  1:2026540/99‑1 (MQ=25)
                          gcgcGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCAttt                                                                           >  2:481668/1‑100 (MQ=255)
                           cgcgTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTta                                                                          >  2:3576045/1‑100 (MQ=255)
                             cgTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaat                                                                        >  2:3594971/1‑100 (MQ=255)
                               ttCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaataa                                                                      >  1:2195467/1‑100 (MQ=255)
                               ttCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaataa                                                                      <  2:1850135/100‑1 (MQ=255)
                                tCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATaaa                                                                     <  2:3596558/100‑1 (MQ=255)
                                 cTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAAt                                                                    >  2:1934281/1‑100 (MQ=255)
                                  ttGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAAtt                                                                   <  2:2888522/100‑1 (MQ=255)
                                   tGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAAtt                                                                   <  1:1486741/99‑1 (MQ=255)
                                    gCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTAc                                                                 >  1:1002195/1‑100 (MQ=255)
                                      tATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAg                                                               <  2:3308353/100‑1 (MQ=255)
                                      tATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAg                                                               <  1:2080165/100‑1 (MQ=255)
                                      tATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAg                                                               <  1:1786257/100‑1 (MQ=255)
                                       aTTGGCCATTACACGGACGGGGGCAATTTCTCTCTATCGGCGTTCCGAAAGATGCCGGATTCGGCGTGAACGCCTTAGCCGGCATTTAATAAATTACAg                                                               >  2:2879166/1‑99 (MQ=25)
                                       aTTGGCAATTACCGGGCCGGGGGCAATGTCTCTCGATTGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCGTTATCCGGCATTTAATAAATTACAg                                                               >  2:2319328/1‑99 (MQ=255)
                                        ttGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGcc                                                             >  2:2378393/1‑100 (MQ=255)
                                          ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                           >  2:1758007/1‑100 (MQ=255)
                                          ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                           >  2:1743055/1‑100 (MQ=255)
                                          ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                           >  1:2632090/1‑100 (MQ=255)
                                          ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                           >  1:205365/1‑100 (MQ=255)
                                           gCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAct                                                          <  1:3378928/100‑1 (MQ=255)
                                            cAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCActc                                                         >  1:1182657/1‑100 (MQ=255)
                                              aTTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTa                                                       <  1:3554789/100‑1 (MQ=255)
                                                tACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGCTGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                     >  2:3302396/1‑100 (MQ=255)
                                                tACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                     >  1:333429/1‑100 (MQ=255)
                                                tACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                     >  1:2673940/1‑100 (MQ=255)
                                                tACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                     >  2:3504511/1‑100 (MQ=255)
                                                tACCAGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATACGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                     <  2:227265/100‑1 (MQ=255)
                                                  ccgggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCg                                                   <  2:2103326/100‑1 (MQ=255)
                                                   cgggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGc                                                  <  1:91421/100‑1 (MQ=255)
                                                    gggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCa                                                 <  1:459222/100‑1 (MQ=255)
                                                       ccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACg                                              <  2:2540341/100‑1 (MQ=255)
                                                       ccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACg                                              <  1:113975/100‑1 (MQ=255)
                                                       ccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCAGGCATTTAATAAATTACAGCCACTCTACCCGCAACg                                              >  2:1392408/1‑100 (MQ=255)
                                                            ggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtc                                         >  2:1899567/1‑100 (MQ=255)
                                                                  ttCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGaa                                   <  2:3441282/100‑1 (MQ=255)
                                                                   tctctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAg                                  <  1:1981084/100‑1 (MQ=255)
                                                                   tcgctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAg                                  <  2:1251517/97‑1 (MQ=255)
                                                                    ctctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgt                                 >  1:2771813/1‑100 (MQ=255)
                                                                     tctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgtt                                >  2:408008/1‑100 (MQ=255)
                                                                      ctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgtt                                >  1:2284051/1‑99 (MQ=255)
                                                                       tcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttgt                              <  2:3271873/100‑1 (MQ=255)
                                                                        cGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttgtt                             >  1:3500099/1‑100 (MQ=255)
                                                                        cGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttgtt                             <  1:2968137/100‑1 (MQ=255)
                                                                         gATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGttgttg                            >  1:2731484/1‑100 (MQ=255)
                                                                         gATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGttgttg                            <  1:3035263/100‑1 (MQ=255)
                                                                         gATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGttgttg                            >  2:1461347/1‑100 (MQ=255)
                                                                         gATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGttgttg                            <  2:856408/100‑1 (MQ=255)
                                                                          aTCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGa                           <  1:498731/100‑1 (MQ=255)
                                                                           tCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGAt                          <  2:1137907/100‑1 (MQ=255)
                                                                                gTTAGGAAAGATGCCGCATGCGGCGTGAACGCCTTATCCGGCATTTAATATATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTcccc                     <  2:2130160/100‑1 (MQ=255)
                                                                                gTTAGGAAAGATGCCAGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTcccc                     <  1:3596094/100‑1 (MQ=255)
                                                                                  tttattaaaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTg                   >  1:1305916/9‑100 (MQ=21)
                                                                                  tAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTg                   >  2:2827830/1‑100 (MQ=255)
                                                                                  tAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTg                   <  2:714418/100‑1 (MQ=255)
                                                                                     attaaaTGCCGGCTAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTg                >  2:1530994/6‑100 (MQ=25)
                                                                                        aaaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTa             >  1:2032150/3‑100 (MQ=25)
                                                                                        aGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTa             >  1:738948/1‑100 (MQ=255)
                                                                                         aaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAg            >  2:3353089/2‑100 (MQ=25)
                                                                                          aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc           <  1:1654540/100‑1 (MQ=255)
                                                                                          aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc           >  1:2446112/1‑100 (MQ=255)
                                                                                          aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc           >  2:3171093/1‑100 (MQ=25)
                                                                                          aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTCAAGTTGTTGATTCCCCTGGTACTAgc           >  2:1098755/1‑100 (MQ=25)
                                                                                           tGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgcg          >  1:3193199/1‑100 (MQ=255)
                                                                                            gCCGGATGCGGCGTGAACGCCTTATCCGGCATTTACTAAATTACAGCAACTCTACCCGCACCGACATCTCTGAAGTTATTGATTCCCCTGGTGCTAgcgc         >  2:2525754/1‑100 (MQ=38)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCATTAAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGCTTCCCCTGGTGCTAgcgc         >  2:527507/1‑100 (MQ=25)
                                                                                             ccGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCa        >  1:3161011/1‑100 (MQ=255)
                                                                                             ccGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCa        >  2:2445834/1‑100 (MQ=255)
                                                                                             ccGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCATGGTGCTAGCGCa        >  2:1705191/1‑100 (MQ=25)
                                                                                             ccGGATAAGGCGTTAACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTGGTTGATTCCCCTGGTGCTAGCGCa        >  2:297578/1‑100 (MQ=37)
                                                                                              cGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCaa       >  1:1888998/1‑100 (MQ=255)
                                                                                              cGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCa        >  1:1364726/1‑99 (MQ=255)
                                                                                              cGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCaa       >  2:2299070/1‑100 (MQ=37)
                                                                                               ggATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAAt      <  2:133086/100‑1 (MQ=255)
                                                                                                   gcggcgTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATgaga  <  1:2078665/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGA  >  NC_000913/2673684‑2673882

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