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.

AGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑TTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACA  >  NC_000913/2333745‑2333938
                                                                                                     |                                                                                               
agaTGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGat‑a                                                                                                  >  1:972613/1‑100 (MQ=255)
 gaTGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGat‑a‑t                                                                                                >  2:419768/1‑100 (MQ=255)
   tGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt                                                                                              >  1:1004343/1‑100 (MQ=255)
   tGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt                                                                                              <  1:2753351/100‑1 (MQ=255)
   tGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGACTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt                                                                                              <  2:2395886/100‑1 (MQ=255)
    gAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt‑t                                                                                            <  2:2615259/100‑1 (MQ=255)
      aaTTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt‑ttt                                                                                          <  1:2793573/100‑1 (MQ=255)
      aaTTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt‑ttt                                                                                          >  2:2412172/1‑100 (MQ=255)
      aaTTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt‑ttt                                                                                          <  1:3049601/100‑1 (MQ=255)
      aaTTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt‑ttt                                                                                          >  2:2863708/1‑100 (MQ=255)
       aTTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt‑tttt                                                                                         <  2:365156/100‑2 (MQ=255)
        ttGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt‑ttttc                                                                                        <  1:869790/100‑3 (MQ=255)
         tGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt‑ttttca                                                                                       >  1:268518/1‑97 (MQ=255)
          gCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt‑ttttcag                                                                                      <  2:912474/100‑5 (MQ=255)
          gCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt‑ttttcag                                                                                      >  2:2635873/1‑96 (MQ=255)
           cTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑ttt‑ttttcagc                                                                                     <  1:419769/100‑6 (MQ=255)
              gATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCggg                                                                                   <  2:1813502/100‑1 (MQ=255)
               aTGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTTAGCgggg                                                                                  >  1:1635048/1‑100 (MQ=255)
               aTGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCgggg                                                                                  <  2:534540/100‑1 (MQ=255)
               aTGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCgggg                                                                                  <  2:1004342/100‑1 (MQ=255)
               aTGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCgggg                                                                                  <  1:2851905/100‑1 (MQ=255)
                 gcgcCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGaa                                                                                >  2:1913664/1‑100 (MQ=255)
                  cgcCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCgggg                                                                                  >  2:3275095/1‑97 (MQ=255)
                  cgcCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGaaa                                                                               >  2:2918625/1‑100 (MQ=255)
                    ccTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAAc                                                                             >  2:195737/1‑100 (MQ=255)
                    ccTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAAc                                                                             >  1:1721029/1‑100 (MQ=255)
                     cTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACt                                                                            <  1:3321329/100‑1 (MQ=255)
                     cTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACt                                                                            <  2:2541366/100‑1 (MQ=255)
                        ccAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGgt                                                                         <  1:2182531/100‑1 (MQ=255)
                        ccAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGgt                                                                         <  1:195738/100‑1 (MQ=255)
                         cAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCTTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGgtt                                                                        <  1:3234893/100‑1 (MQ=255)
                         cAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGgtt                                                                        >  2:1179598/1‑100 (MQ=255)
                         cAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGgtt                                                                        >  2:1424554/1‑100 (MQ=255)
                         cAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGgtt                                                                        >  1:3095486/1‑100 (MQ=255)
                           gTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGgttgt                                                                      >  1:911197/1‑100 (MQ=255)
                             tttACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTc                                                                    <  1:1274626/100‑1 (MQ=255)
                             tttACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTc                                                                    <  2:2562320/100‑1 (MQ=255)
                              ttACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCt                                                                   >  1:2104325/1‑100 (MQ=255)
                              ttACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCt                                                                   >  1:590084/1‑100 (MQ=255)
                               tACCCCTTGCTGCGGGTCGAACCGTTTAGCCCGCCACAGTGGATAAAGGGACTTTTTCACCGAGGGATAA‑TTT‑TTTTAAGGGGGAAAAATGggtg                                                                       >  2:1582446/1‑95 (MQ=255)
                               tACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTg                                                                  >  1:2217873/1‑100 (MQ=255)
                                nccccTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGt                                                                 >  1:2594276/2‑100 (MQ=255)
                                aCCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGt                                                                 <  1:2678125/100‑1 (MQ=255)
                                 ccccTTGCTGCGGCTCGAACGGTTTAGCCTGCCATATTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCGGTa                                                                >  2:280701/1‑100 (MQ=255)
                                 ccccTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTa                                                                >  1:1481610/1‑100 (MQ=255)
                                  cccTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTAc                                                               <  1:68873/100‑1 (MQ=255)
                                   ccTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACa                                                              >  2:1162245/1‑100 (MQ=255)
                                    cTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTAc                                                               <  2:2839151/98‑1 (MQ=255)
                                    cTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACAt                                                             <  1:2705087/100‑1 (MQ=255)
                                    cTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACAt                                                             >  1:2512181/1‑100 (MQ=255)
                                       gctgcGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAAc                                                          <  1:2863383/100‑1 (MQ=255)
                                         tgcGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACgg                                                        >  2:488048/1‑100 (MQ=255)
                                         tgcGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACgg                                                        >  2:2633426/1‑100 (MQ=255)
                                         tgcGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACgg                                                        <  2:411000/100‑1 (MQ=255)
                                          gcGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACggg                                                       <  1:1051931/100‑1 (MQ=255)
                                              cTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAgg                                                   >  1:1059291/1‑100 (MQ=255)
                                              cTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAgg                                                   <  2:1671264/100‑1 (MQ=255)
                                                cGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGtt                                                 <  1:3180066/100‑1 (MQ=255)
                                                 ttccgatctAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGttt                                                <  1:58854/92‑1 (MQ=255)
                                                 gAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGttt                                                <  1:729150/100‑1 (MQ=255)
                                                  aaCGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTta                                               <  2:590083/100‑1 (MQ=255)
                                                  aaCGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTta                                               >  2:1936463/1‑100 (MQ=255)
                                                    cGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTtaat                                             <  2:1764833/100‑1 (MQ=255)
                                                      gTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTtaataa                                           >  1:148581/1‑100 (MQ=255)
                                                        ttAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTaataat                                          <  2:584639/99‑1 (MQ=255)
                                                         tAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATgg                                        <  1:1127192/100‑1 (MQ=255)
                                                           gCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCa                                      >  2:1432044/1‑100 (MQ=255)
                                                           gCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCa                                      >  2:2267910/1‑100 (MQ=255)
                                                             cTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTTTTTTTTTTTAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCAtt                                    >  2:200143/1‑100 (MQ=255)
                                                              tGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTg                                   >  2:1733366/1‑100 (MQ=255)
                                                               gcgaTAGTGGAGAACGGGACTGTTGCACCAGGGGAG‑A‑TTTGTTGCAGCGGGGAAAAGTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATg                                         <  2:3235642/90‑1 (MQ=255)
                                                               gCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGc                                  >  1:493110/1‑100 (MQ=255)
                                                                ccATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGcc                                 >  1:1744507/1‑100 (MQ=255)
                                                                ccATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGcc                                 <  2:512683/100‑1 (MQ=255)
                                                                 cATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑TTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGa                               <  1:226132/100‑1 (MQ=255)
                                                                  aTAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCg                                >  1:1617532/1‑99 (MQ=255)
                                                                  aTAGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGa                               <  2:3090690/100‑1 (MQ=255)
                                                                    aGTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGc                             <  2:138812/100‑1 (MQ=255)
                                                                     gTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCt                            <  1:2846550/100‑1 (MQ=255)
                                                                     gTGGATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCt                            <  1:980505/100‑1 (MQ=255)
                                                                        gATAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTa                         <  1:1700739/100‑1 (MQ=255)
                                                                         aTAACGGAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGtn                         <  2:546889/99‑2 (MQ=255)
                                                                              ggAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCgg                   <  1:3277129/100‑1 (MQ=255)
                                                                               gAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCggg                  <  2:1947895/100‑1 (MQ=255)
                                                                               gAACTGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCggg                  <  2:1057983/100‑1 (MQ=255)
                                                                                   tGTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAAc              <  2:2512173/100‑1 (MQ=255)
                                                                                    gTTTCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAAc              >  2:1255008/1‑99 (MQ=255)
                                                                                      ttCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGcc           <  1:2507445/100‑1 (MQ=255)
                                                                                      ttCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGcc           <  1:2633434/100‑1 (MQ=255)
                                                                                       tCACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCt          >  1:1510160/1‑100 (MQ=255)
                                                                                        cACCAGGGGAT‑ATTTT‑TTTTAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTg         >  2:161567/1‑100 (MQ=255)
                                                                                        cACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTg         <  2:2217867/100‑1 (MQ=255)
                                                                                        cACCAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTg         >  2:70556/1‑100 (MQ=255)
                                                                                          ccAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAg       <  1:2918636/100‑1 (MQ=255)
                                                                                          ccAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAg       <  1:568200/100‑1 (MQ=255)
                                                                                           cAGGGGAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGt      <  1:2805156/100‑1 (MQ=255)
                                                                                             ggggAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTca    <  1:280702/100‑1 (MQ=255)
                                                                                              gggAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTcac   <  2:268517/100‑1 (MQ=255)
                                                                                              gggAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTcac   <  1:161566/100‑1 (MQ=255)
                                                                                               ggAT‑ATTTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTcaca  <  2:1059290/100‑1 (MQ=255)
                                                                                                     |                                                                                               
AGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGAT‑A‑TTT‑TTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACA  >  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: ‑