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.

GATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑C‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAA  >  NC_000913/2767403‑2767598
                                                                                                   |                                                                                                      
gATATTAACTTGTCCAATATGATCAAATAGCATTACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑c                                                                                                     <  1:2633826/100‑1 (MQ=255)
gATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑c                                                                                                     <  1:1927565/100‑1 (MQ=255)
 atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑c‑c                                                                                                   <  2:429770/100‑1 (MQ=255)
 atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑c‑c                                                                                                   <  2:1450165/100‑1 (MQ=255)
 atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑c‑c                                                                                                   <  1:1756397/100‑1 (MQ=255)
  tatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCTTCACAACGTCCTGCAGAGGCAACAAGTTTTCCCCACTGCACCCGCGACTCAAt‑ttatc                                                                                                      >  1:1165460/1‑96 (MQ=255)
   atTAACTTGTCCAATATGATCAAATAGCATTAAACACCCCTCACAACGTCCTGCATAGGGAAAAAGATATACTCTGTGAACCCACGACTAAAT‑AT‑c‑a‑c‑c‑c                                                                                               >  1:2561111/1‑100 (MQ=255)
    ttAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑c‑c‑c‑cc                                                                                              <  2:1612785/100‑1 (MQ=255)
      aaCTTGTCCAATATGATCAAATAGCATTACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑c‑c‑c‑cccc                                                                                            <  1:2384497/100‑2 (MQ=255)
         ttGTCCAATATGGTCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑C‑C‑C‑C‑CCCcttt                                                                                         <  2:1311008/100‑2 (MQ=255)
          tgcccAATATGATCAAATAGCATTCACCCCCCCTCACAACGTCCGGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑c‑c‑c‑cccc                                                                                            <  1:550059/93‑2 (MQ=255)
             ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAAc                                                                                      <  2:314881/100‑1 (MQ=255)
                atatGATCAAATAGCATTAACCCCCCCTCACAACATCCCGCATAGGGAACACGTCTTCCCCTGTGCACCCACGCCTAAAT‑TTAC‑C‑C‑C‑CCCTTTAAACt                                                                                     >  1:2373881/1‑98 (MQ=255)
                 tatGATCAAATAGCATTAACGCCCTCTCACAACGTAATGCATGAGGAACACGGTATCCCCAGTGCACGCACGACTAAAT‑TTCC‑C‑C‑C‑CCCGGTAAACTtct                                                                                  >  1:1440841/1‑100 (MQ=255)
                 tatGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTtct                                                                                  <  1:107018/100‑1 (MQ=255)
                  atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTtctc                                                                                 <  2:638849/100‑1 (MQ=255)
                   tGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCa                                                                                <  2:356936/100‑1 (MQ=255)
                   tGAGCAAATAGCAGTAGCCCCCCCTCACAACGTCCTGCATAGGGAACCCGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTtctc                                                                                 <  1:1409288/99‑1 (MQ=255)
                     aTCAAATAGCATTAACCCTCCTTCCCGGCTTTCTGCATCCGAAACACGATTACCCCGGAACCCCCAAGCCTAAAT‑TC‑C‑C‑CAC‑CCCTTCAATCTTCTTaat                                                                              >  1:1232461/1‑100 (MQ=255)
                     aTCAAATAGCATTAACCCACCCTCACAGCGTCCTGCATAGGGAACAAATTATCCCCTGTGAACCCACGACTAAAT‑TTAG‑C‑C‑T‑CCATTTAAACTgctg                                                                                 >  1:2741594/1‑96 (MQ=255)
                        aaaTAGCATTAACCCCCCATCACAACGTCCTGCGTAGGGAACACCTTTTCCCCAGTACAACCACGACTAAAT‑TT‑CAC‑C‑C‑CCCTTTAAACTCCTCaaaaat                                                                           >  1:1922879/1‑100 (MQ=255)
                         aaTAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCaataat                                                                           <  1:172858/99‑1 (MQ=255)
                         aaTAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATg                                                                          <  2:1798952/100‑1 (MQ=255)
                          aTAGCATTAACCCCCCCTCAAAACGCCCTGCATCAGGAACAAGTCCTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTTTCaataat                                                                           >  1:1986688/1‑98 (MQ=255)
                           tAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTc                                                                        <  2:747616/100‑1 (MQ=255)
                             gCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCCCCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATACTGTCAc                                                                      >  1:2843630/1‑100 (MQ=255)
                               aTTAACCCCCCCTCACAACGTACTGCATAGGGAAAACGTTTTCCCCTGTGCACCAACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATCATGTccc                                                                      >  1:1393637/1‑96 (MQ=255)
                                ttAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAC‑TT‑C‑CAC‑C‑CCCTTTACACTTCTAAAAAATCTCAccca                                                                   >  1:2152124/1‑98 (MQ=255)
                                   aCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAAttt                                                                <  2:2007628/100‑1 (MQ=255)
                                     ccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTc                                                              <  1:2207350/100‑1 (MQ=255)
                                      cctccTCACAACGTCCTGCATAGGGAACCCGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                             <  1:1385155/97‑1 (MQ=255)
                                      ccccccCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                             <  1:1642875/100‑1 (MQ=255)
                                      cccccTCCCACCGTCCTGCATACGTAACCCCTTTTCCCCTGTGCACCCACGACTAACT‑TTCC‑C‑C‑C‑CCCTTTAACCTTCTCAATAATGTCACCAATTTTCt                                                             <  1:1541809/100‑1 (MQ=255)
                                      cccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                             >  1:1052534/1‑100 (MQ=255)
                                        cccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg                                                           <  2:590703/100‑1 (MQ=255)
                                        cccTCACAACGCCCTGCATAGGCAACACGTCTTCCCCTGTGCACCCACGAATAACTCAC‑C‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTg                                                            <  1:2453800/99‑1 (MQ=255)
                                         ccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGa                                                          <  2:815033/100‑1 (MQ=255)
                                          cTCACGACGTCCTGCATAGGGAACACGTTTTCCCCGGTGCGCCCACGACTAAAT‑GTCC‑C‑C‑C‑CCCTTTAAATTTCTCAATAATGTCACCAATTTTCTGGaa                                                         >  1:1152165/1‑100 (MQ=255)
                                           tCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaaa                                                        >  2:1111822/1‑100 (MQ=255)
                                            cacaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTATAT‑TTCC‑C‑C‑C‑CCCCTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                       >  1:1290281/1‑100 (MQ=255)
                                            cacaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                       <  2:751225/100‑1 (MQ=255)
                                              caACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATcc                                                     >  2:1560352/1‑100 (MQ=255)
                                               atcGTCCTGCATAGGGAACACGTTTTCCCCCATGCACCCACGACTAAAT‑TT‑C‑C‑C‑CTCCCTTTATCCTTCTCCATAATGTTACCAATTTTCTGGaaagc                                                      <  1:2652907/96‑3 (MQ=255)
                                                aCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                                   <  2:2274057/100‑1 (MQ=255)
                                                 cGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTat                                                  <  2:843050/100‑1 (MQ=255)
                                                  gTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatg                                                 <  2:1503487/100‑1 (MQ=255)
                                                   tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga                                                <  1:1366931/100‑1 (MQ=255)
                                                   tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga                                                <  2:2806108/100‑1 (MQ=255)
                                                   tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga                                                >  1:2639653/1‑100 (MQ=255)
                                                    ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAAGTTCTCAATAATGCCACCAATTTTCTGGAAATCCTatgat                                               >  1:1326627/1‑100 (MQ=255)
                                                    ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgat                                               >  2:2666426/1‑100 (MQ=255)
                                                     cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                              >  1:2348852/1‑100 (MQ=255)
                                                      tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                             <  2:69592/100‑1 (MQ=255)
                                                      tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                             >  2:2292590/1‑100 (MQ=255)
                                                      tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCAGGAAATCCTatgatg                                              >  2:272252/1‑99 (MQ=255)
                                                      tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTGTAAACTTCTCAAGAATGTCACCAATGTTCTGGAAATCCTATGATGc                                             >  2:272260/1‑100 (MQ=255)
                                                      tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑CAC‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                             <  1:78977/100‑1 (MQ=255)
                                                       gCATGGGGAACACGTTTTCGCCTGTGCACTCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAGTTTTCTGGAAATCCTATGATGcc                                            >  1:2767262/1‑100 (MQ=255)
                                                       gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                            <  2:2736589/100‑1 (MQ=255)
                                                         caggggAAAACGTGTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCtt                                          >  1:25442/4‑100 (MQ=255)
                                                          tAGGGAACACGTTTT‑CCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                        >  1:2649641/1‑100 (MQ=255)
                                                             ggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                      >  2:377042/1‑100 (MQ=255)
                                                              gAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                      <  1:372406/99‑1 (MQ=255)
                                                              gAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                     >  2:1727356/1‑100 (MQ=255)
                                                              gAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                     <  1:2283435/100‑1 (MQ=255)
                                                                acacGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGt                                   <  2:1052536/100‑1 (MQ=255)
                                                                acacGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGt                                   <  2:1048829/100‑1 (MQ=255)
                                                                 cacGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                                  <  2:43487/100‑1 (MQ=255)
                                                                 cacGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                                  <  1:2862131/100‑1 (MQ=255)
                                                                 cacGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑C‑C‑CGT‑CGCTTTAGACTTCTCAAGAATGAAAACAAGATAAAGGAAATCCTATGATGCCTTTCCCGGTa                                  >  1:2043845/1‑100 (MQ=255)
                                                                   cGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa                                >  2:914573/1‑100 (MQ=255)
                                                                   cGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTCAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa                                >  2:1392039/1‑100 (MQ=255)
                                                                    gTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                               <  1:108617/100‑1 (MQ=255)
                                                                      tttCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgg                             <  1:272250/100‑1 (MQ=255)
                                                                      tttCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgg                             <  2:1032873/100‑1 (MQ=255)
                                                                      tttCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgg                             <  2:1632424/100‑1 (MQ=255)
                                                                        tCCCCTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTg                           >  2:119766/1‑100 (MQ=255)
                                                                         ccccTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATTATGCCTTTCCCGGTAAAt                               >  2:399154/1‑95 (MQ=255)
                                                                         ccccTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTg                           >  1:1733797/1‑99 (MQ=255)
                                                                         ccccTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                          >  2:35715/1‑100 (MQ=255)
                                                                          cccTGTGCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                          <  1:298162/99‑1 (MQ=255)
                                                                             tgtgCACCCACGACTAACT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                      <  1:2712480/100‑1 (MQ=255)
                                                                             tgtgCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                      <  2:2843640/100‑1 (MQ=255)
                                                                               tgCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                    <  2:1688988/100‑1 (MQ=255)
                                                                                gCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGt                   >  2:1923657/1‑100 (MQ=255)
                                                                                gCACCCACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGt                   >  1:234456/1‑100 (MQ=255)
                                                                                 cACCCACGACTAAAT‑TT‑C‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTc                 >  2:2320058/1‑100 (MQ=255)
                                                                                  aaccACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTc                 <  1:34497/98‑1 (MQ=255)
                                                                                   cccACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa                <  2:1244496/100‑1 (MQ=255)
                                                                                   cccACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa                >  2:1308249/1‑100 (MQ=255)
                                                                                   cccACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa                <  1:296597/100‑1 (MQ=255)
                                                                                   cccACGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa                <  2:1641025/100‑1 (MQ=255)
                                                                                       cGAGTAAAT‑TTCC‑C‑C‑C‑CCCTTTAATCTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa            <  1:949007/100‑1 (MQ=255)
                                                                                       cGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa            <  1:377040/100‑1 (MQ=255)
                                                                                       cGACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa            >  2:1496105/1‑100 (MQ=255)
                                                                                        gACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAg           <  2:1387155/100‑1 (MQ=255)
                                                                                        gACTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAg           >  2:687287/1‑100 (MQ=255)
                                                                                          cTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt         <  2:2117195/100‑1 (MQ=255)
                                                                                          cTAAAT‑TTCC‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt         >  2:1581892/1‑100 (MQ=255)
                                                                                                 tt‑C‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCaa  >  2:1624652/1‑100 (MQ=255)
                                                                                                   |                                                                                                      
GATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑C‑C‑C‑C‑CCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAA  >  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: ‑