Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,333,850 (T)6→7 pseudogene (553/4605 nt) yfaS ← pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,333,8441.T95.4% 379.1 / 9.7 108pseudogene (559/4605 nt)yfaSpseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base T (45/58);  total (47/61)
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 127 total aligned reads displayed.

AGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACA  >  NC_000913/2333745‑2333938
                                                                                                     |                                                                                               
agaTGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGata                                                                                                  >  1:972613/1‑100 (MQ=255)
 gaTGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGatat                                                                                                >  2:419768/1‑100 (MQ=255)
   tGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAttt                                                                                              >  1:1004343/1‑100 (MQ=255)
   tGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAttt                                                                                              <  1:2753351/100‑1 (MQ=255)
   tGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGACTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAttt                                                                                              <  2:2395886/100‑1 (MQ=255)
    gAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAtttt                                                                                            <  2:2615259/100‑1 (MQ=255)
      aaTTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAtttttt                                                                                          <  1:2793573/100‑1 (MQ=255)
      aaTTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAtttttt                                                                                          >  2:2412172/1‑100 (MQ=255)
      aaTTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAtttttt                                                                                          <  1:3049601/100‑1 (MQ=255)
      aaTTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAtttttt                                                                                          >  2:2863708/1‑100 (MQ=255)
       aTTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAttttttt                                                                                         <  2:365156/100‑2 (MQ=255)
        ttGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAtttttttc                                                                                        <  1:869790/100‑3 (MQ=255)
         tGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAtttttttca                                                                                       >  1:268518/1‑97 (MQ=255)
          gCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAtttttttcag                                                                                      <  2:912474/100‑5 (MQ=255)
          gCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAtttttttcag                                                                                      >  2:2635873/1‑96 (MQ=255)
           cTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATAtttttttcagc                                                                                     <  1:419769/100‑6 (MQ=255)
              gATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCggg                                                                                   <  2:1813502/100‑1 (MQ=255)
               aTGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTTAGCgggg                                                                                  >  1:1635048/1‑100 (MQ=255)
               aTGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCgggg                                                                                  <  2:534540/100‑1 (MQ=255)
               aTGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCgggg                                                                                  <  2:1004342/100‑1 (MQ=255)
               aTGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCgggg                                                                                  <  1:2851905/100‑1 (MQ=255)
                 gcgcCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGaa                                                                                >  2:1913664/1‑100 (MQ=255)
                  cgcCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCgggg                                                                                  >  2:3275095/1‑97 (MQ=255)
                  cgcCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGaaa                                                                               >  2:2918625/1‑100 (MQ=255)
                    ccTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAAc                                                                             >  2:195737/1‑100 (MQ=255)
                    ccTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAAc                                                                             >  1:1721029/1‑100 (MQ=255)
                     cTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACt                                                                            <  1:3321329/100‑1 (MQ=255)
                     cTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACt                                                                            <  2:2541366/100‑1 (MQ=255)
                        ccAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGgt                                                                         <  1:2182531/100‑1 (MQ=255)
                        ccAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGgt                                                                         <  1:195738/100‑1 (MQ=255)
                         cAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCTTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGgtt                                                                        <  1:3234893/100‑1 (MQ=255)
                         cAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGgtt                                                                        >  2:1179598/1‑100 (MQ=255)
                         cAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGgtt                                                                        >  2:1424554/1‑100 (MQ=255)
                         cAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGgtt                                                                        >  1:3095486/1‑100 (MQ=255)
                           gTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGgttgt                                                                      >  1:911197/1‑100 (MQ=255)
                             tttACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTc                                                                    <  1:1274626/100‑1 (MQ=255)
                             tttACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTc                                                                    <  2:2562320/100‑1 (MQ=255)
                              ttACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCt                                                                   >  1:2104325/1‑100 (MQ=255)
                              ttACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCt                                                                   >  1:590084/1‑100 (MQ=255)
                               tACCCCTTGCTGCGGGTCGAACCGTTTAGCCCGCCACAGTGGATAAAGGGACTTTTTCACCGAGGGATAATTTTTTTAAGGGGGAAAAATGggtg                                                                       >  2:1582446/1‑95 (MQ=255)
                               tACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTg                                                                  >  1:2217873/1‑100 (MQ=255)
                                nccccTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGt                                                                 >  1:2594276/2‑100 (MQ=255)
                                aCCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGt                                                                 <  1:2678125/100‑1 (MQ=255)
                                 ccccTTGCTGCGGCTCGAACGGTTTAGCCTGCCATATTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCGGTa                                                                >  2:280701/1‑100 (MQ=255)
                                 ccccTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTa                                                                >  1:1481610/1‑100 (MQ=255)
                                  cccTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTAc                                                               <  1:68873/100‑1 (MQ=255)
                                   ccTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACa                                                              >  2:1162245/1‑100 (MQ=255)
                                    cTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTAc                                                               <  2:2839151/98‑1 (MQ=255)
                                    cTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACAt                                                             <  1:2705087/100‑1 (MQ=255)
                                    cTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACAt                                                             >  1:2512181/1‑100 (MQ=255)
                                       gctgcGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAAc                                                          <  1:2863383/100‑1 (MQ=255)
                                         tgcGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACgg                                                        >  2:488048/1‑100 (MQ=255)
                                         tgcGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACgg                                                        >  2:2633426/1‑100 (MQ=255)
                                         tgcGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACgg                                                        <  2:411000/100‑1 (MQ=255)
                                          gcGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACggg                                                       <  1:1051931/100‑1 (MQ=255)
                                              cTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAgg                                                   >  1:1059291/1‑100 (MQ=255)
                                              cTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAgg                                                   <  2:1671264/100‑1 (MQ=255)
                                                cGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGtt                                                 <  1:3180066/100‑1 (MQ=255)
                                                 ttccgatctAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGttt                                                <  1:58854/92‑1 (MQ=255)
                                                 gAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGttt                                                <  1:729150/100‑1 (MQ=255)
                                                  aaCGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTta                                               <  2:590083/100‑1 (MQ=255)
                                                  aaCGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTta                                               >  2:1936463/1‑100 (MQ=255)
                                                    cGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTtaat                                             <  2:1764833/100‑1 (MQ=255)
                                                      gTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTtaataa                                           >  1:148581/1‑100 (MQ=255)
                                                        ttAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTaataat                                          <  2:584639/99‑1 (MQ=255)
                                                         tAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATgg                                        <  1:1127192/100‑1 (MQ=255)
                                                           gCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCa                                      >  2:1432044/1‑100 (MQ=255)
                                                           gCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCa                                      >  2:2267910/1‑100 (MQ=255)
                                                             cTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTTTTTTTTTTTAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCAtt                                    >  2:200143/1‑100 (MQ=255)
                                                              tGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTg                                   >  2:1733366/1‑100 (MQ=255)
                                                               gcgaTAGTGGAGAACGGGACTGTTGCACCAGGGGAGATTTGTTGCAGCGGGGAAAAGTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATg                                         <  2:3235642/90‑1 (MQ=255)
                                                               gCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGc                                  >  1:493110/1‑100 (MQ=255)
                                                                ccATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGcc                                 >  1:1744507/1‑100 (MQ=255)
                                                                ccATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGcc                                 <  2:512683/100‑1 (MQ=255)
                                                                 cATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGa                               <  1:226132/100‑1 (MQ=255)
                                                                  aTAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCg                                >  1:1617532/1‑99 (MQ=255)
                                                                  aTAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGa                               <  2:3090690/100‑1 (MQ=255)
                                                                    aGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGc                             <  2:138812/100‑1 (MQ=255)
                                                                     gTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCt                            <  1:2846550/100‑1 (MQ=255)
                                                                     gTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCt                            <  1:980505/100‑1 (MQ=255)
                                                                        gATAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTa                         <  1:1700739/100‑1 (MQ=255)
                                                                         aTAACGGAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGtn                         <  2:546889/99‑2 (MQ=255)
                                                                              ggAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCgg                   <  1:3277129/100‑1 (MQ=255)
                                                                               gAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCggg                  <  2:1947895/100‑1 (MQ=255)
                                                                               gAACTGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCggg                  <  2:1057983/100‑1 (MQ=255)
                                                                                   tGTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAAc              <  2:2512173/100‑1 (MQ=255)
                                                                                    gTTTCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAAc              >  2:1255008/1‑99 (MQ=255)
                                                                                      ttCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGcc           <  1:2507445/100‑1 (MQ=255)
                                                                                      ttCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGcc           <  1:2633434/100‑1 (MQ=255)
                                                                                       tCACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCt          >  1:1510160/1‑100 (MQ=255)
                                                                                        cACCAGGGGATATTTTTTTTAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTg         >  2:161567/1‑100 (MQ=255)
                                                                                        cACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTg         <  2:2217867/100‑1 (MQ=255)
                                                                                        cACCAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTg         >  2:70556/1‑100 (MQ=255)
                                                                                          ccAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAg       <  1:2918636/100‑1 (MQ=255)
                                                                                          ccAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAg       <  1:568200/100‑1 (MQ=255)
                                                                                           cAGGGGATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGt      <  1:2805156/100‑1 (MQ=255)
                                                                                             ggggATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTca    <  1:280702/100‑1 (MQ=255)
                                                                                              gggATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTcac   <  2:268517/100‑1 (MQ=255)
                                                                                              gggATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTcac   <  1:161566/100‑1 (MQ=255)
                                                                                               ggATATTTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTcaca  <  2:1059290/100‑1 (MQ=255)
                                                                                                     |                                                                                               
AGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACA  >  NC_000913/2333745‑2333938

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