Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,673,784 C→A 15.1% intergenic (+16/+32) 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,7840CA15.1% 520.8 / 63.3 177intergenic (+16/+32)yphA/yphBDoxX family inner membrane protein/mutarotase superfamily protein, YphB family
Reads supporting (aligned to +/- strand):  ref base C (79/71);  new base A (17/10);  total (96/81)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 195 total aligned reads displayed.

TTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGAT  >  NC_000913/2673685‑2673883
                                                                                                   |                                                                                                   
ttCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgc                                                                                                     <  2:3377265/100‑1 (MQ=255)
 tCTGGAAGAATGTCAGTATCGCTGGCGCTTTCTTGCTATTGGTAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcg                                                                                                    <  2:1667456/100‑1 (MQ=255)
  cTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcgg                                                                                                   <  2:2602000/100‑1 (MQ=255)
  cTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCCATTGGCACTTCCCTGGCCGGGGGCAATTTCTCTCGATCCTCGTTAGGAAAGCTGCCGGATgcgg                                                                                                   >  2:2197125/1‑100 (MQ=255)
   tGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcggc                                                                                                  <  1:3294800/100‑1 (MQ=255)
   tGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcgg                                                                                                   >  2:2107995/1‑99 (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:1071293/100‑1 (MQ=255)
            cgatctATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                         <  1:1744739/95‑1 (MQ=255)
               aGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTa                                                                                      >  1:2303642/1‑100 (MQ=255)
                gTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTAt                                                                                     <  2:80874/100‑1 (MQ=255)
                     gCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                <  1:2272269/100‑1 (MQ=37)
                       tGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                                              <  1:981872/100‑1 (MQ=255)
                       tGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCan                                                                              <  1:101675/100‑2 (MQ=25)
                        ggCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAtc                                                                             <  1:2149358/100‑2 (MQ=25)
                        ggCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAtc                                                                             <  1:549596/100‑2 (MQ=25)
                        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)
                               tCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATaaa                                                                      <  2:3596558/100‑1 (MQ=255)
                               tCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATaaa                                                                      <  2:1909869/100‑1 (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:3298160/100‑1 (MQ=255)
                                     tATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAg                                                                <  2:3308353/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                                                            >  1:205365/1‑100 (MQ=255)
                                         ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                            >  2:1743055/1‑100 (MQ=255)
                                         ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                            >  1:2632090/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)
                                               tACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                      >  1:2673940/1‑100 (MQ=255)
                                               tACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                      >  2:3504511/1‑100 (MQ=255)
                                               tACCAGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATACGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                      <  2:227265/100‑1 (MQ=255)
                                                aCCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAccc                                                     >  1:714421/1‑100 (MQ=255)
                                                 ccgggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCg                                                    >  1:49970/1‑100 (MQ=255)
                                                   gggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCa                                                  <  1:459222/100‑1 (MQ=255)
                                                      ccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACg                                               <  1:113975/100‑1 (MQ=255)
                                                      ccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACg                                               <  2:2540341/100‑1 (MQ=255)
                                                      ccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCAGGCATTTAATAAATTACAGCCACTCTACCCGCAACg                                               >  2:1392408/1‑100 (MQ=255)
                                                        gggggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGAc                                             >  2:768264/1‑100 (MQ=255)
                                                           ggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtc                                          >  2:1097014/1‑100 (MQ=255)
                                                           ggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtc                                          <  2:1851811/100‑1 (MQ=255)
                                                                 ttCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGaa                                    <  2:3441282/100‑1 (MQ=255)
                                                                  tctctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAg                                   <  2:1053697/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)
                                                                      tcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttgt                               >  1:3193637/1‑100 (MQ=255)
                                                                       cGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttgtt                              <  1:2968137/100‑1 (MQ=255)
                                                                       cGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttgtt                              >  1:3500099/1‑100 (MQ=255)
                                                                        gATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGttgttg                             >  1:2731484/1‑100 (MQ=255)
                                                                        gATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGttgttg                             >  2:1461347/1‑100 (MQ=255)
                                                                        gATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGttgttg                             <  1:3035263/100‑1 (MQ=255)
                                                                        gATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGttgttg                             <  2:856408/100‑1 (MQ=255)
                                                                         aTCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGa                            <  1:498731/100‑1 (MQ=255)
                                                                          tCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGAt                           <  2:1137907/100‑1 (MQ=255)
                                                                             gCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTcc                        >  1:62921/1‑100 (MQ=255)
                                                                              cGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTccc                       <  1:859627/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)
                                                                                   ggAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGt                  >  1:2329396/1‑100 (MQ=255)
                                                                                    attaaaTGCCGGCTAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTg                 >  2:1530994/6‑100 (MQ=25)
                                                                                       aaaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTa              >  1:2032150/3‑100 (MQ=25)
                                                                                        aaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAg             >  2:1767858/2‑100 (MQ=25)
                                                                                        aaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAg             >  2:3353089/2‑100 (MQ=25)
                                                                                         aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc            <  1:1654540/100‑1 (MQ=255)
                                                                                         aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc            >  1:2446112/1‑100 (MQ=255)
                                                                                         aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTTATTCCCCTGGTGCTAgc            >  2:200873/1‑100 (MQ=25)
                                                                                         aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc            >  2:3171093/1‑100 (MQ=25)
                                                                                          tGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgcg           <  1:787/100‑1 (MQ=255)
                                                                                           gCCGGATGCGGCGTGAACGCCTTATCCGGCATTTACTAAATTACAGCAACTCTACCCGCACCGACATCTCTGAAGTTATTGATTCCCCTGGTGCTAgcgc          >  2:2525754/1‑100 (MQ=38)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCATTAAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGCTTCCCCTGGTGCTAgcgc          >  2:527507/1‑100 (MQ=25)
                                                                                            ccGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCa         >  2:3251664/1‑100 (MQ=255)
                                                                                            ccGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCa         >  2:2445834/1‑100 (MQ=255)
                                                                                            ccGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCa         >  1:3161011/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)
                                                                                                   cggcgTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGAt  >  2:60833/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGAT  >  NC_000913/2673685‑2673883

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: