Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,507:1 (C)7→8 100% intergenic (+152/‑203) rnlB → / → yfjP CP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,5001.C89.6% 295.0 / 28.0 96intergenic (+145/‑210)rnlB/yfjPCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (6/4);  new base C (35/49);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

GATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAA  >  NC_000913/2767403‑2767598
                                                                                                   |                                                                                                      
gATATTAACTTGTCCAATATGATCAAATAGCATTACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTcc                                                                                                     <  1:2633826/100‑1 (MQ=255)
gATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTcc                                                                                                     <  1:1927565/100‑1 (MQ=255)
 atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                                                                   <  2:429770/100‑1 (MQ=255)
 atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                                                                   <  2:1450165/100‑1 (MQ=255)
 atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                                                                   <  1:1756397/100‑1 (MQ=255)
  tatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCTTCACAACGTCCTGCAGAGGCAACAAGTTTTCCCCACTGCACCCGCGACTCAAtttatc                                                                                                      >  1:1165460/1‑96 (MQ=255)
   atTAACTTGTCCAATATGATCAAATAGCATTAAACACCCCTCACAACGTCCTGCATAGGGAAAAAGATATACTCTGTGAACCCACGACTAAATATcaccc                                                                                               >  1:2561111/1‑100 (MQ=255)
    ttAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTcccccc                                                                                              <  2:1612785/100‑1 (MQ=255)
      aaCTTGTCCAATATGATCAAATAGCATTACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTcccccccc                                                                                            <  1:2384497/100‑2 (MQ=255)
         ttGTCCAATATGGTCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCcttt                                                                                         <  2:1311008/100‑2 (MQ=255)
          tgcccAATATGATCAAATAGCATTCACCCCCCCTCACAACGTCCGGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTcccccccc                                                                                            <  1:550059/93‑2 (MQ=255)
             ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAAc                                                                                      <  2:314881/100‑1 (MQ=255)
                atatGATCAAATAGCATTAACCCCCCCTCACAACATCCCGCATAGGGAACACGTCTTCCCCTGTGCACCCACGCCTAAATTTACCCCCCCTTTAAACt                                                                                     >  1:2373881/1‑98 (MQ=255)
                 tatGATCAAATAGCATTAACGCCCTCTCACAACGTAATGCATGAGGAACACGGTATCCCCAGTGCACGCACGACTAAATTTCCCCCCCCGGTAAACTtct                                                                                  >  1:1440841/1‑100 (MQ=255)
                 tatGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTtct                                                                                  <  1:107018/100‑1 (MQ=255)
                  atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTtctc                                                                                 <  2:638849/100‑1 (MQ=255)
                   tGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCa                                                                                <  2:356936/100‑1 (MQ=255)
                   tGAGCAAATAGCAGTAGCCCCCCCTCACAACGTCCTGCATAGGGAACCCGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTtctc                                                                                 <  1:1409288/99‑1 (MQ=255)
                     aTCAAATAGCATTAACCCTCCTTCCCGGCTTTCTGCATCCGAAACACGATTACCCCGGAACCCCCAAGCCTAAATTCCCCACCCCTTCAATCTTCTTaat                                                                              >  1:1232461/1‑100 (MQ=255)
                     aTCAAATAGCATTAACCCACCCTCACAGCGTCCTGCATAGGGAACAAATTATCCCCTGTGAACCCACGACTAAATTTAGCCTCCATTTAAACTgctg                                                                                 >  1:2741594/1‑96 (MQ=255)
                        aaaTAGCATTAACCCCCCATCACAACGTCCTGCGTAGGGAACACCTTTTCCCCAGTACAACCACGACTAAATTTCACCCCCCTTTAAACTCCTCaaaaat                                                                           >  1:1922879/1‑100 (MQ=255)
                         aaTAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCaataat                                                                           <  1:172858/99‑1 (MQ=255)
                         aaTAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATg                                                                          <  2:1798952/100‑1 (MQ=255)
                          aTAGCATTAACCCCCCCTCAAAACGCCCTGCATCAGGAACAAGTCCTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTTTCaataat                                                                           >  1:1986688/1‑98 (MQ=255)
                           tAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTc                                                                        <  2:747616/100‑1 (MQ=255)
                             gCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCCCCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATACTGTCAc                                                                      >  1:2843630/1‑100 (MQ=255)
                               aTTAACCCCCCCTCACAACGTACTGCATAGGGAAAACGTTTTCCCCTGTGCACCAACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATCATGTccc                                                                      >  1:1393637/1‑96 (MQ=255)
                                ttAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAACTTCCACCCCCTTTACACTTCTAAAAAATCTCAccca                                                                   >  1:2152124/1‑98 (MQ=255)
                                   aCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAAttt                                                                <  2:2007628/100‑1 (MQ=255)
                                     ccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTc                                                              <  1:2207350/100‑1 (MQ=255)
                                      cctccTCACAACGTCCTGCATAGGGAACCCGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                             <  1:1385155/97‑1 (MQ=255)
                                      ccccccCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                             <  1:1642875/100‑1 (MQ=255)
                                      cccccTCCCACCGTCCTGCATACGTAACCCCTTTTCCCCTGTGCACCCACGACTAACTTTCCCCCCCCTTTAACCTTCTCAATAATGTCACCAATTTTCt                                                             <  1:1541809/100‑1 (MQ=255)
                                      cccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                             >  1:1052534/1‑100 (MQ=255)
                                        cccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg                                                           <  2:590703/100‑1 (MQ=255)
                                        cccTCACAACGCCCTGCATAGGCAACACGTCTTCCCCTGTGCACCCACGAATAACTCACCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTg                                                            <  1:2453800/99‑1 (MQ=255)
                                         ccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGa                                                          <  2:815033/100‑1 (MQ=255)
                                          cTCACGACGTCCTGCATAGGGAACACGTTTTCCCCGGTGCGCCCACGACTAAATGTCCCCCCCCTTTAAATTTCTCAATAATGTCACCAATTTTCTGGaa                                                         >  1:1152165/1‑100 (MQ=255)
                                           tCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaaa                                                        >  2:1111822/1‑100 (MQ=255)
                                            cacaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTATATTTCCCCCCCCCTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                       >  1:1290281/1‑100 (MQ=255)
                                            cacaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                       <  2:751225/100‑1 (MQ=255)
                                              caACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATcc                                                     >  2:1560352/1‑100 (MQ=255)
                                               atcGTCCTGCATAGGGAACACGTTTTCCCCCATGCACCCACGACTAAATTTCCCCTCCCTTTATCCTTCTCCATAATGTTACCAATTTTCTGGaaagc                                                      <  1:2652907/96‑3 (MQ=255)
                                                aCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                                   <  2:2274057/100‑1 (MQ=255)
                                                 cGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTat                                                  <  2:843050/100‑1 (MQ=255)
                                                  gTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatg                                                 <  2:1503487/100‑1 (MQ=255)
                                                   tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga                                                <  1:1366931/100‑1 (MQ=255)
                                                   tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga                                                <  2:2806108/100‑1 (MQ=255)
                                                   tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga                                                >  1:2639653/1‑100 (MQ=255)
                                                    ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAAGTTCTCAATAATGCCACCAATTTTCTGGAAATCCTatgat                                               >  1:1326627/1‑100 (MQ=255)
                                                    ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgat                                               >  2:2666426/1‑100 (MQ=255)
                                                     cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                              >  1:2348852/1‑100 (MQ=255)
                                                      tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                             <  2:69592/100‑1 (MQ=255)
                                                      tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                             >  2:2292590/1‑100 (MQ=255)
                                                      tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCAGGAAATCCTatgatg                                              >  2:272252/1‑99 (MQ=255)
                                                      tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTGTAAACTTCTCAAGAATGTCACCAATGTTCTGGAAATCCTATGATGc                                             >  2:272260/1‑100 (MQ=255)
                                                      tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCACCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                             <  1:78977/100‑1 (MQ=255)
                                                       gCATGGGGAACACGTTTTCGCCTGTGCACTCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAGTTTTCTGGAAATCCTATGATGcc                                            >  1:2767262/1‑100 (MQ=255)
                                                       gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                            <  2:2736589/100‑1 (MQ=255)
                                                         caggggAAAACGTGTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCtt                                          >  1:25442/4‑100 (MQ=255)
                                                          tAGGGAACACGTTTTCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                        >  1:2649641/1‑100 (MQ=255)
                                                             ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                      >  2:377042/1‑100 (MQ=255)
                                                              gAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                      <  1:372406/99‑1 (MQ=255)
                                                              gAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                     >  2:1727356/1‑100 (MQ=255)
                                                              gAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                     <  1:2283435/100‑1 (MQ=255)
                                                                acacGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGt                                   <  2:1052536/100‑1 (MQ=255)
                                                                acacGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGt                                   <  2:1048829/100‑1 (MQ=255)
                                                                 cacGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                                  <  2:43487/100‑1 (MQ=255)
                                                                 cacGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                                  <  1:2862131/100‑1 (MQ=255)
                                                                 cacGTTTTCCCCTGTGCACCCACGACTAAATTTCCCGTCGCTTTAGACTTCTCAAGAATGAAAACAAGATAAAGGAAATCCTATGATGCCTTTCCCGGTa                                  >  1:2043845/1‑100 (MQ=255)
                                                                   cGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa                                >  2:914573/1‑100 (MQ=255)
                                                                   cGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTCAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa                                >  2:1392039/1‑100 (MQ=255)
                                                                    gTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                               <  1:108617/100‑1 (MQ=255)
                                                                      tttCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgg                             <  1:272250/100‑1 (MQ=255)
                                                                      tttCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgg                             <  2:1032873/100‑1 (MQ=255)
                                                                      tttCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgg                             <  2:1632424/100‑1 (MQ=255)
                                                                        tCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTg                           >  2:119766/1‑100 (MQ=255)
                                                                         ccccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATTATGCCTTTCCCGGTAAAt                               >  2:399154/1‑95 (MQ=255)
                                                                         ccccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTg                           >  1:1733797/1‑99 (MQ=255)
                                                                         ccccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                          >  2:35715/1‑100 (MQ=255)
                                                                          cccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                          <  1:298162/99‑1 (MQ=255)
                                                                             tgtgCACCCACGACTAACTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                      <  1:2712480/100‑1 (MQ=255)
                                                                             tgtgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                      <  2:2843640/100‑1 (MQ=255)
                                                                               tgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                    <  2:1688988/100‑1 (MQ=255)
                                                                                gCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGt                   >  2:1923657/1‑100 (MQ=255)
                                                                                gCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGt                   >  1:234456/1‑100 (MQ=255)
                                                                                 cACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTc                 >  2:2320058/1‑100 (MQ=255)
                                                                                  aaccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTc                 <  1:34497/98‑1 (MQ=255)
                                                                                   cccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa                <  2:1244496/100‑1 (MQ=255)
                                                                                   cccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa                >  2:1308249/1‑100 (MQ=255)
                                                                                   cccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa                <  1:296597/100‑1 (MQ=255)
                                                                                   cccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa                <  2:1641025/100‑1 (MQ=255)
                                                                                       cGAGTAAATTTCCCCCCCCTTTAATCTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa            <  1:949007/100‑1 (MQ=255)
                                                                                       cGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa            <  1:377040/100‑1 (MQ=255)
                                                                                       cGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa            >  2:1496105/1‑100 (MQ=255)
                                                                                        gACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAg           <  2:1387155/100‑1 (MQ=255)
                                                                                        gACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAg           >  2:687287/1‑100 (MQ=255)
                                                                                          cTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt         <  2:2117195/100‑1 (MQ=255)
                                                                                          cTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt         >  2:1581892/1‑100 (MQ=255)
                                                                                                 ttCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCaa  >  2:1624652/1‑100 (MQ=255)
                                                                                                   |                                                                                                      
GATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAA  >  NC_000913/2767403‑2767598

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