Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,767,507 (C)7→8 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.C88.1% 342.9 / 41.6 110intergenic (+145/‑210)rnlB/yfjPCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (7/6);  new base C (37/59);  total (45/65)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

TGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTC  >  NC_000913/2767402‑2767593
                                                                                                    |                                                                                                 
tGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑c                                                                                                 <  1:377368/100‑1 (MQ=255)
tGATATTAACTTGTCCAATATGATCAAATAGCATTAACCACCCCTCACAACGTCCTGCATAGTGAACTCGATTTCCCCTGTGCACCAACGACTAAAatt                                                                                                    >  1:558785/1‑99 (MQ=255)
 gATATTAACTTGTCCAATATGATCAAATAGCATTANCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑cc                                                                                               <  1:1932155/100‑1 (MQ=255)
 gATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑cc                                                                                               <  2:2314965/100‑1 (MQ=255)
   tatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑cccc                                                                                            <  2:2961140/100‑1 (MQ=255)
     ttccgatctCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑cccccc                                                                                          <  2:1376051/92‑1 (MQ=255)
     ttAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTCCAAAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑cccccc                                                                                          >  2:2517324/1‑100 (MQ=255)
       aaCTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCACTGTACGGAACGCGGACACACCTGTTCACCCACGACTAACTTT‑‑c                                                                                                 >  1:1724056/1‑93 (MQ=255)
           tGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCcttta                                                                                   <  1:1006365/100‑3 (MQ=255)
           tGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCcttta                                                                                   <  1:1352224/100‑3 (MQ=255)
           tGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCcttta                                                                                   <  2:654050/100‑3 (MQ=255)
           tGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCAAAGGGCACACGTGTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCcttta                                                                                   >  2:79130/1‑98 (MQ=255)
            gTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTaa                                                                                   <  1:798970/100‑1 (MQ=255)
              ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAAc                                                                                 <  2:2196273/100‑1 (MQ=255)
              ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAAc                                                                                 <  1:2177144/100‑1 (MQ=255)
              ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAAc                                                                                 <  2:967116/100‑1 (MQ=255)
              ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTCTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAAc                                                                                 <  2:695843/100‑1 (MQ=255)
              ccAATATGATCAAATAGCATTAACACCCCCTCACAACCTCCCGCATAGAGAACACGTTTTCCCCTGTGCACCCACGCCTAAATCT‑‑caccccagatcggaa                                                                                 >  1:1376059/1‑91 (MQ=255)
               cAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACt                                                                                <  2:271548/100‑1 (MQ=255)
               cAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACt                                                                                <  1:662508/100‑1 (MQ=255)
               cAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑ccccccccctttaaac                                                                               <  2:2195752/100‑10 (MQ=255)
                 tcctGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCCTTTAAACtt                                                                               <  1:2604244/97‑1 (MQ=255)
                  tatGATCAAATAGCATTAACCACCCCTCCCAACGTCCTGCATATTGAACTCGTTTTCCCCTGTGCACCCACTTCTAAATTT‑‑CCGTCCCCCTTaaa                                                                                  >  1:105698/1‑95 (MQ=255)
                   atGATCAAATAGCATTAACCCCCCCTCACAGCGTCCTGCATAGTGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTtctn                                                                            <  2:297163/100‑2 (MQ=255)
                      aTCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCaat                                                                         <  2:254537/100‑1 (MQ=255)
                       tCAAATAGCATTAACCCGCCCTCACAACGTCCTGCAAAGGGAAGACGTTTTCCCCTGTTCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCacta                                                                        >  2:2891506/1‑100 (MQ=255)
                        cAAATAGCATTAACCCCCCCTCACGACGTCCTGCATAGCCAACACGCCTTCCTCTGTGCCCCCACGACGAAATCAC‑‑CCCCCCCTTTAAACTTCTCaata                                                                        <  2:1142523/99‑1 (MQ=255)
                        cAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCACCCCCTTTAAACTTCTAaataa                                                                       >  2:2198005/1‑100 (MQ=255)
                        cAAATAGCATTAACCCCCACTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCaataa                                                                       >  2:1747439/1‑100 (MQ=255)
                         aaaTAGCATTAACCTCACATCACAACATACTGCATAAGTAATTCTATTTATCATTTGCTCCCTCGACTAAATTT‑‑CCTCCCCCTTTCATCtttt                                                                             >  2:2698630/1‑91 (MQ=255)
                             aGCATTAACCCCCCCGCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCa                                                                  <  2:1622159/100‑1 (MQ=255)
                               cATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATg                                                                     >  2:1756263/1‑95 (MQ=255)
                                     cccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAAtttt                                                          <  1:1985745/100‑1 (MQ=255)
                                      ccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCCCCAATTTTc                                                         >  2:1427387/1‑100 (MQ=255)
                                      ccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTc                                                         <  1:2330164/100‑1 (MQ=255)
                                      ccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTc                                                         <  1:2436635/100‑1 (MQ=255)
                                       cccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCt                                                        <  2:2623367/100‑1 (MQ=255)
                                        ttccgatctaCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTg                                                       <  1:2356965/91‑1 (MQ=255)
                                          ccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGa                                                     >  2:239597/1‑100 (MQ=255)
                                           cTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaa                                                    <  2:127698/100‑1 (MQ=255)
                                             cacaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                  >  2:2835925/1‑100 (MQ=255)
                                              acaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCCTTAAACTTCCTAACAATGTCACCAATTTTCTGGAAATc                                                 >  1:924591/1‑100 (MQ=255)
                                                aaCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCt                                               <  2:1669507/100‑1 (MQ=255)
                                                 aCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgagatcggaa                                              >  2:2356951/1‑91 (MQ=255)
                                                 aCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                              <  1:1917979/100‑1 (MQ=255)
                                                 aCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                              <  1:2333658/100‑1 (MQ=255)
                                                 aCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                              <  1:98803/100‑1 (MQ=255)
                                                 aCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                              >  1:2335246/1‑100 (MQ=255)
                                                 aCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                              <  1:1395619/100‑1 (MQ=255)
                                                  cGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                              <  1:1567956/99‑1 (MQ=255)
                                                     ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgat                                          <  1:2145921/100‑1 (MQ=255)
                                                       tGCATAGGGAACACGTTTTCCCCTTTGCACCCACGTCTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                         <  2:137251/99‑1 (MQ=255)
                                                       tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                        >  2:793293/1‑100 (MQ=255)
                                                       tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                        >  1:2962149/1‑100 (MQ=255)
                                                       tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                        <  1:1747450/100‑1 (MQ=255)
                                                        gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                       >  1:2000004/1‑100 (MQ=255)
                                                        gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                       >  1:2332358/1‑100 (MQ=255)
                                                        gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                       >  1:757806/1‑100 (MQ=255)
                                                        gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                       >  2:261281/1‑100 (MQ=255)
                                                        gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACAAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                       >  1:2110147/1‑100 (MQ=255)
                                                           tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                    <  2:2297082/100‑1 (MQ=255)
                                                              ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CACGCCACTTTAAACCTCTCAAAAATCCTACCAATTTACTCGAAATCCTATGATGCCTTccc                                  >  2:281436/1‑99 (MQ=255)
                                                               gAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTATCAATAATGTCGCCAATTTTTTGGAAATCCTATGATGCCTTTCCCg                                >  1:1752298/1‑100 (MQ=255)
                                                                aaCACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCgg                               >  2:2006245/1‑100 (MQ=255)
                                                                aaCACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCgg                               <  1:2198017/100‑1 (MQ=255)
                                                                aaCACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                >  2:944864/1‑99 (MQ=255)
                                                                 acacGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGt                              <  1:2619647/100‑1 (MQ=255)
                                                                  cacGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                             >  2:96833/1‑100 (MQ=255)
                                                                  cacGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                             <  1:1410248/100‑1 (MQ=255)
                                                                  cacGTTTTCCCCTGTGCACACACGACTAAATTTCCCCCCCCTTTAAACTTCTCACTAATGTCACCAATTTCCTGAAAATCCTATGATCCCTTTCACGGtc                             >  2:2664517/1‑99 (MQ=255)
                                                                   acGTTTTCCCCTGAGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaa                            <  1:630088/100‑1 (MQ=255)
                                                                     gTTTTCCCCTGTGCACCCACGACTAACTTT‑‑CCCCCCTCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTACAt                          >  2:1405685/1‑100 (MQ=255)
                                                                     gTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CTCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                          <  2:2828778/100‑1 (MQ=255)
                                                                      ttttCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAATCTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                         <  1:2698646/100‑1 (MQ=255)
                                                                      ttttCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                         <  2:94742/100‑1 (MQ=255)
                                                                         tCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATggcg                      <  1:2026549/100‑3 (MQ=255)
                                                                          ccccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                     <  1:1036948/100‑1 (MQ=255)
                                                                           cccTGTGCACCCACGACTAAATTT‑‑CCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                  >  1:1175036/1‑100 (MQ=255)
                                                                            ccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtg                   <  2:1674370/100‑1 (MQ=255)
                                                                            ccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtg                   >  2:2975533/1‑100 (MQ=255)
                                                                             cTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                  <  1:2695385/100‑1 (MQ=255)
                                                                              tgtgCACCCACGACTAAATTTCCCCCCCCTTTATACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                 <  1:319578/100‑1 (MQ=255)
                                                                              tgtgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                 <  2:2547092/100‑1 (MQ=255)
                                                                              tgtgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                 <  2:17828/100‑1 (MQ=255)
                                                                              tgtgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                  >  2:1706107/1‑99 (MQ=255)
                                                                                tgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg               <  2:641314/100‑1 (MQ=255)
                                                                                 gCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGt              <  1:1273929/100‑1 (MQ=255)
                                                                                 gCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGt              <  1:2517339/100‑1 (MQ=255)
                                                                                    cccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa           <  2:2650752/100‑1 (MQ=255)
                                                                                    cccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa           >  1:2670073/1‑100 (MQ=255)
                                                                                    cccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAAGGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa           >  1:2964641/1‑100 (MQ=255)
                                                                                      cACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGa         >  2:195503/1‑100 (MQ=255)
                                                                                      cACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGa         >  1:2529553/1‑100 (MQ=255)
                                                                                        cGACTAAATTTCCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaa        >  1:2660137/1‑100 (MQ=255)
                                                                                        cGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa       <  1:1092843/100‑1 (MQ=255)
                                                                                         gACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAg      <  2:2725643/100‑1 (MQ=255)
                                                                                          aCTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa     >  2:1698196/1‑100 (MQ=255)
                                                                                          aCTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa     >  2:184371/1‑100 (MQ=255)
                                                                                          aCTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa     >  1:77312/1‑100 (MQ=255)
                                                                                           cTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt    >  2:533362/1‑100 (MQ=255)
                                                                                                    |                                                                                                 
TGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTC  >  NC_000913/2767402‑2767593

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: