Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,518,687 A→C E74A (GAA→GCA)  yfeC → DUF1323 family putative DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,518,6870AC97.2% 354.4 / 0.6 109E74A (GAA→GCA) yfeCDUF1323 family putative DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (54/52);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 109 total aligned reads displayed.

GTTCACGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGAAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTGTTAC  >  NC_000913/2518590‑2518782
                                                                                                 |                                                                                               
gTTCTCGTCAATGAACAGGTTCGTGAATATAGTCGCAATGCTGATCGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAg                                                                                               <  2:3201747/100‑1 (MQ=255)
gTTCACGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAg                                                                                               >  2:2460368/1‑100 (MQ=255)
  tCACGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGGGCCGCTTGCAGTa                                                                                             >  1:2441213/1‑100 (MQ=255)
     ngTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTg                                                                                          >  1:3249633/2‑100 (MQ=255)
     cGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTg                                                                                          >  2:905703/1‑100 (MQ=255)
     cGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTg                                                                                          >  2:868962/1‑100 (MQ=255)
     cGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTg                                                                                          >  1:871372/1‑100 (MQ=255)
     cGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTg                                                                                          >  1:938923/1‑100 (MQ=255)
        cAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTg                                                                                       <  1:3037754/100‑1 (MQ=255)
         aaTGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTgg                                                                                      <  2:1240434/100‑1 (MQ=255)
           tGAACAGCTTCGTGAATTTATTCGCAATGCTGTACGTCCAGAAGCTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTg                                                                                    <  1:3746166/100‑1 (MQ=255)
            gAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGa                                                                                   <  2:1052503/100‑1 (MQ=255)
            gAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGa                                                                                   <  2:2595900/100‑1 (MQ=255)
            gAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGa                                                                                   <  1:3783481/100‑1 (MQ=255)
            gAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGGCAGGGAGAAGCGACTGCCCTTTCTGGTGCTGCGCCGCTTGCAGTATTGCTGGTGa                                                                                   >  2:2624616/1‑100 (MQ=255)
              aCAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGaca                                                                                 >  2:287178/1‑100 (MQ=255)
               cAGGTCCGCGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGCAGAAGCACCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGacac                                                                                <  1:1018025/100‑1 (MQ=255)
                 ngTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTg                                                                              >  1:3234537/2‑100 (MQ=255)
                 ggTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTg                                                                              >  2:379975/1‑100 (MQ=255)
                 ggTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTg                                                                              <  2:1420037/100‑1 (MQ=255)
                  gTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTgg                                                                             >  2:3095873/1‑100 (MQ=255)
                   ttCGGGGATATATTGGCAATGCTGAACGTCCAGAGGGTGAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGc                                                                            <  1:3351936/100‑1 (MQ=255)
                    tCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCg                                                                           <  1:3870093/100‑1 (MQ=255)
                      gTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGGGaa                                                                         <  1:3143236/100‑1 (MQ=255)
                       tGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGaca                                                                        >  2:7916/1‑98 (MQ=255)
                        gAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAag                                                                       <  1:2960062/100‑1 (MQ=255)
                          atatatTCGCAATGCTGTACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAagag                                                                     <  1:3749367/100‑1 (MQ=255)
                          atatatTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAagag                                                                     <  1:1853535/100‑1 (MQ=255)
                             tatTCGCCATGCAGAACGTCCAGAATGTCAGGGAGAAGCGCCTGCCCTTACTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAACGAGGTg                                                                  >  2:3848119/1‑100 (MQ=255)
                             tatTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATg                                                                  <  1:291949/100‑1 (MQ=255)
                             tatTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATg                                                                  <  2:3985138/100‑1 (MQ=255)
                             tatTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATg                                                                  <  2:1507020/100‑1 (MQ=255)
                              atTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGa                                                                 >  1:2487513/1‑100 (MQ=255)
                              atTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGAAGTATTGCTGGTGACACTGGCGAAAGAGATGa                                                                 >  2:749355/1‑100 (MQ=255)
                               ttCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGAc                                                                <  1:72218/100‑1 (MQ=255)
                               ttCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGAc                                                                <  2:852912/100‑1 (MQ=255)
                               ttCGCAATGCTGAACGTCCAGAAGGTCAGGGAAAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGAc                                                                <  2:970649/100‑1 (MQ=255)
                                  gCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGcc                                                             <  1:797306/100‑1 (MQ=255)
                                   cAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCa                                                            >  2:1537753/1‑100 (MQ=255)
                                     aTGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGt                                                          >  1:931271/1‑100 (MQ=255)
                                      tGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGtt                                                         <  1:2353460/100‑1 (MQ=255)
                                       gCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGGTGACGCCAGTTg                                                        >  1:1288235/1‑100 (MQ=255)
                                         tGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGaa                                                      >  2:2679601/1‑100 (MQ=255)
                                          gAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAGAGAGATGACGCCAGTTGAAc                                                     >  2:1338487/1‑100 (MQ=255)
                                          gAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAAc                                                     <  1:1754851/100‑1 (MQ=255)
                                          gAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAAc                                                     >  2:1576552/1‑100 (MQ=255)
                                          gAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCACCGCTTGCAGTATGGCAAGTGACACTGGCGAAAGAGATGGCGCCAGTTGAAc                                                     >  2:3315904/1‑100 (MQ=255)
                                               tCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACaaaaa                                                <  2:514491/100‑1 (MQ=255)
                                               tCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACaaaaa                                                >  1:1964713/1‑100 (MQ=255)
                                                ccAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAAc                                               <  1:98158/100‑1 (MQ=255)
                                                ccAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCGGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAAc                                               >  2:617743/1‑100 (MQ=255)
                                                 cAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACa                                              >  1:2276079/1‑100 (MQ=255)
                                                 cAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACa                                              <  1:2414109/100‑1 (MQ=255)
                                                      ggTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTa                                         <  2:3630011/100‑1 (MQ=255)
                                                      ggTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCGGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTa                                         >  2:2772690/1‑100 (MQ=255)
                                                       gTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTAc                                        >  1:1972705/1‑100 (MQ=255)
                                                       gTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTAc                                        >  2:3033848/1‑100 (MQ=255)
                                                         cAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACAt                                      >  1:3856677/1‑100 (MQ=255)
                                                         cAGGGAGAAGCGCCTGCCCTGTCGGCTGATGCACCGCTTGCAGTATTGCTGTTGACACTGGCGAAAGAGATGACGCAAAAAGAACAAAAACAGTTTACAt                                      <  2:527469/100‑1 (MQ=255)
                                                           gggAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATcc                                    >  2:771259/1‑100 (MQ=255)
                                                             gagaAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCtt                                  >  1:998959/1‑100 (MQ=255)
                                                              agaAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTg                                 <  2:2094620/100‑1 (MQ=255)
                                                              agaAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTg                                 <  1:887689/100‑1 (MQ=255)
                                                                aaGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCt                               <  1:2530091/100‑1 (MQ=255)
                                                                 aGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCtt                              >  1:1107398/1‑100 (MQ=255)
                                                                  gcgcCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTc                             >  1:3192780/1‑100 (MQ=255)
                                                                  gcgcCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTc                             <  1:1420846/100‑1 (MQ=255)
                                                                  gcgcCTGCCCTTTCTGGTGATGCGCCGCTTGAAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTc                             <  1:765360/100‑1 (MQ=255)
                                                                    gcCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTg                           >  1:164540/1‑100 (MQ=255)
                                                                    gcCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTg                           >  1:2397739/1‑100 (MQ=255)
                                                                      cTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCg                         >  2:1677122/1‑100 (MQ=255)
                                                                      cTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCg                         >  1:43079/1‑100 (MQ=255)
                                                                       tGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCgg                        <  2:1822789/100‑1 (MQ=255)
                                                                       tGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCgg                        <  1:52784/100‑1 (MQ=255)
                                                                       tGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCgg                        <  1:16111/100‑1 (MQ=255)
                                                                        gCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCggg                       >  1:665829/1‑100 (MQ=255)
                                                                         cccTTTCTGGTGATGCGCCGCTTGAAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGa                      >  1:2497947/1‑100 (MQ=255)
                                                                          ccTTTCTGGGGATGCGCCGCTTGCGGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGCTTACATCCTTGCTTCTGCGGGaa                     >  2:2459964/1‑100 (MQ=255)
                                                                               cTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGAt                <  1:1088650/100‑1 (MQ=255)
                                                                                tGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGAtt               >  1:3178769/1‑100 (MQ=255)
                                                                                 ggTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTGCATCCTTGCTTCTGCGGGAAGGGATTa              >  2:1739951/1‑100 (MQ=255)
                                                                                 ggTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTa              >  2:1556280/1‑100 (MQ=255)
                                                                                   tGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATc            <  1:3089186/100‑1 (MQ=255)
                                                                                    gATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCg           <  2:3781227/100‑1 (MQ=255)
                                                                                     aTGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCgg          >  1:3845244/1‑100 (MQ=255)
                                                                                      tGCGCCGCTTGCAGTATTGCTGGTGACACTTGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGa         >  2:1465057/1‑100 (MQ=255)
                                                                                      tGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGa         >  2:2977145/1‑100 (MQ=255)
                                                                                      tGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGa         >  2:1306192/1‑100 (MQ=255)
                                                                                       gcgcCGCTTGCAGTAATGATGGTGACGCTGGCGAAAGAGATGATGAAAGATAAAAAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGan        <  2:75203/100‑2 (MQ=255)
                                                                                       gccccgcTTGCAGTATTGCTGGTGACACCGGAGAAAGACATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGAt        <  1:1206951/97‑1 (MQ=255)
                                                                                        ggccgcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGAtt       <  1:4084202/99‑1 (MQ=255)
                                                                                         gccgcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTg      >  1:2832795/1‑100 (MQ=255)
                                                                                         gccgcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTg      <  2:278931/100‑1 (MQ=255)
                                                                                         gccgcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTg      <  1:1254728/100‑1 (MQ=255)
                                                                                         aaagcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAATTAAACAAAAACAGTTTACATTCTTGCGTCTGCGGGAAGGGCTTCTCGGAttc      >  2:196065/4‑99 (MQ=255)
                                                                                          ccgcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTGt     <  1:3469390/100‑1 (MQ=255)
                                                                                          ccgcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTGt     >  1:1304622/1‑100 (MQ=255)
                                                                                          ccgcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTGt     <  2:1654546/100‑1 (MQ=255)
                                                                                            gcTTGCAGTATTGCGGGTGACACTGGCGAAAGAGGTGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTGTTa   <  1:441708/100‑1 (MQ=255)
                                                                                             cTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTGTTAc  <  1:1319183/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GTTCACGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGAAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTGTTAC  >  NC_000913/2518590‑2518782

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: