Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,122,367 (G)7→6 100% intergenic (‑80/+13) hslV ← / ← ftsN peptidase component of the HslUV protease/essential cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,122,3610G.78.2% 377.4 / 145.5 147intergenic (‑74/+19)hslV/ftsNpeptidase component of the HslUV protease/essential cell division protein
Reads supporting (aligned to +/- strand):  ref base G (18/14);  new base . (60/55);  total (78/69)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 158 total aligned reads displayed.

CGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTT  >  NC_000913/4122262‑4122459
                                                                                                   |                                                                                                  
cGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATg                                                                                                    >  1:2436926/1‑100 (MQ=255)
 ggCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATgg                                                                                                   >  2:3528298/1‑100 (MQ=255)
  gCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATggg                                                                                                  <  1:3373403/100‑1 (MQ=255)
  gCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATggg                                                                                                  >  1:3424628/1‑100 (MQ=255)
   cGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATgggg                                                                                                 >  2:2260676/1‑100 (MQ=255)
    ttccgatctcTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATATATggggg                                                                                                <  1:1910964/91‑1 (MQ=255)
    gTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATggggg                                                                                                <  1:2032389/100‑1 (MQ=255)
    gTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATggggg                                                                                                <  2:2786897/100‑1 (MQ=255)
         cTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATggggggattt                                                                                           >  2:2025353/1‑96 (MQ=255)
             cTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAg                                                                                      >  1:598115/1‑100 (MQ=255)
              tataGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAgg                                                                                     <  1:3100558/100‑1 (MQ=255)
              tataGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAgg                                                                                     >  1:3401650/1‑100 (MQ=255)
               ataGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGTTTTTGAggg                                                                                    >  2:1295181/1‑100 (MQ=255)
               ataGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAggg                                                                                    <  1:2162421/100‑1 (MQ=255)
                taGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAggg                                                                                    >  2:243074/1‑100 (MQ=255)
                 aGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGtt                                                                                  >  1:1127367/1‑100 (MQ=255)
                  tctGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGttt                                                                                 <  2:2093771/98‑1 (MQ=255)
                  gTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGttt                                                                                 <  2:3601865/100‑1 (MQ=255)
                  gTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGttt                                                                                 <  2:1758390/100‑1 (MQ=255)
                  gTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGttt                                                                                 >  2:1882725/1‑100 (MQ=255)
                     gTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATATATGGGGCGATTTTGATGGTTTCa                                                                               >  2:386874/1‑99 (MQ=255)
                      tCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAAc                                                                             <  1:2021969/100‑1 (MQ=255)
                       cACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAAcc                                                                            >  1:2313360/1‑100 (MQ=255)
                        aCGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAAcc                                                                            <  2:2445667/100‑1 (MQ=255)
                        aCGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAAcc                                                                            <  1:1343028/100‑1 (MQ=255)
                         cGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAAccc                                                                           >  2:1414225/1‑99 (MQ=255)
                          gAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAAccccc                                                                         <  1:640355/100‑1 (MQ=255)
                           aGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCcg                                                                        >  1:3382915/1‑100 (MQ=255)
                           aGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCcg                                                                        >  1:2859214/1‑100 (MQ=255)
                             cTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCcggc                                                                      >  1:1088234/1‑100 (MQ=255)
                             cTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCcggc                                                                      >  2:1102990/1‑100 (MQ=255)
                                aCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCggcggc                                                                   <  1:1087762/100‑1 (MQ=255)
                                 ccccTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGCTGCGGAGGATTTTGAGGGTTTCAACCCCcgccgg                                                                    >  2:1961784/1‑99 (MQ=255)
                                 ccccTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGgcggcg                                                                  >  1:1391843/1‑100 (MQ=255)
                                  ctcTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGAGGGGGGGGATTTTGAGGGTTTCAACCCCcggcgg                                                                    <  2:3208140/96‑1 (MQ=255)
                                  cccTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGa                                                                 <  1:1450707/100‑1 (MQ=255)
                                   ccTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAg                                                                <  2:1081577/100‑1 (MQ=255)
                                   ccTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAg                                                                <  1:314732/100‑1 (MQ=255)
                                   ccTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAg                                                                <  1:1149275/100‑1 (MQ=255)
                                     ttGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGcc                                                              <  1:3048708/100‑1 (MQ=255)
                                     ttGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCACGGCGGCGAGcc                                                              <  1:2398799/100‑1 (MQ=255)
                                      tGGTTACGAATACAAAGTACAAACCCCGTACAAAAGTACGGGGCTTAATGCAGTTATAGATGGGGGGATTTTGAGTGTTTCAACCACCGGCGGCGAGca                                                              <  2:2575725/99‑2 (MQ=255)
                                       ggTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCg                                                             <  1:674884/100‑1 (MQ=255)
                                       ggTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGa                                                            >  1:404726/1‑100 (MQ=255)
                                       ggTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGa                                                            >  1:70782/1‑100 (MQ=255)
                                          tACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAAt                                                          >  2:3031990/1‑100 (MQ=255)
                                           aCGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGc                                                        <  1:2136926/100‑1 (MQ=255)
                                            cGAATACAGAGTACAAACCACGTACAAAAGAACAGGTCATAATGCAATTATAGATGGTGGTGAGTTTGAGGGTTGCAACCCCCTGCGGCGatc                                                               >  2:3038500/1‑91 (MQ=255)
                                              aaTACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAg                                                      >  2:3229442/1‑99 (MQ=255)
                                              aaTACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGt                                                     >  1:3069122/1‑100 (MQ=255)
                                               aTACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGtt                                                    <  1:2937810/100‑1 (MQ=255)
                                               aTACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGtt                                                    <  2:901995/100‑1 (MQ=255)
                                                 aCAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTg                                                  <  1:221916/100‑1 (MQ=255)
                                                 aCAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTg                                                  <  2:3134681/100‑1 (MQ=255)
                                                  cAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTg                                                  <  1:2573616/100‑1 (MQ=255)
                                                  cAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGt                                                 <  1:3312090/100‑1 (MQ=255)
                                                  cAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGt                                                 >  2:674154/1‑100 (MQ=255)
                                                  cAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGt                                                 >  2:3124486/1‑100 (MQ=255)
                                                   agagTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCGCGGCGGCGAGCCGAATGCAGttt                                                   >  2:286236/1‑97 (MQ=255)
                                                   agagTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTa                                                <  1:2283679/100‑1 (MQ=255)
                                                    gagTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTa                                                <  1:538811/100‑1 (MQ=255)
                                                    gagTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTAt                                               <  2:2730792/100‑1 (MQ=255)
                                                    gagTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTAt                                               <  1:2300418/100‑1 (MQ=255)
                                                    gagTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTAt                                               >  2:2300402/1‑100 (MQ=255)
                                                     cgTACACCCCACGTACAATAGTACGGACCACGATACAACTACAGATCGGCGTATTCTGATGATTTCCCCCCCCGCCGGCGAGCCGAATGCAGCTCGTAtc                                              <  2:2458827/99‑2 (MQ=255)
                                                       tACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGa                                             >  2:861623/1‑99 (MQ=255)
                                                         cAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGAcc                                           >  1:1556237/1‑100 (MQ=255)
                                                           aaCCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGc                                        <  1:3293570/100‑1 (MQ=255)
                                                            aCCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGcc                                       >  2:3586532/1‑100 (MQ=255)
                                                              cccGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCa                                      >  1:456543/1‑100 (MQ=255)
                                                              cccGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCAt                                     >  2:1529882/1‑100 (MQ=255)
                                                                cGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCt                                   <  1:2503831/100‑1 (MQ=255)
                                                                 gTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCtt                                  <  1:1258581/100‑1 (MQ=255)
                                                                   aCAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCa                                <  2:205476/100‑1 (MQ=255)
                                                                    cAAAGGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCaa                               >  1:1128640/1‑100 (MQ=255)
                                                                    cAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGGGAGCCGAATGCAGTTTGTATGACCCGCCCTCTTCaa                               >  2:3207743/1‑100 (MQ=255)
                                                                    cAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCaa                               >  2:2500870/1‑100 (MQ=255)
                                                                    cAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCa                                >  2:1845800/1‑99 (MQ=255)
                                                                      aaaGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCg                            >  1:2031880/1‑100 (MQ=255)
                                                                       aaGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCg                            <  1:2382571/100‑1 (MQ=255)
                                                                       aaGGACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGACATCTTCAACCg                            <  2:2314349/100‑1 (MQ=255)
                                                                             ggggCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAg                      >  2:44974/1‑100 (MQ=255)
                                                                             ggggCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAg                      >  1:2336868/1‑100 (MQ=255)
                                                                              gggCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAgg                     >  1:941921/1‑100 (MQ=255)
                                                                              gggCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAgg                     >  1:3003980/1‑100 (MQ=255)
                                                                               ggCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAggg                    <  2:328358/100‑1 (MQ=255)
                                                                                gCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAggg                    >  1:969160/1‑100 (MQ=255)
                                                                                     aTGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTg               >  2:1832233/1‑100 (MQ=255)
                                                                                     aTGCAATTATAGATGGGGGGATTTTGAGGGTTTACACCCCCGGCGGCGCGCCGAATGCAGTTTGTATGCCCCGCCATCTTCAACCGATTGAGGGTGCTGt              <  1:2500886/100‑1 (MQ=255)
                                                                                      tGCCATTATAGATGGGGGGATTTTGAGGGCTTCAACCCCCGCCGGCGAGCCGAATGCAGTTCGTATGCCCCGCCCTCTTCCACCGCTCGAGGGCGCTg               >  2:2896784/1‑98 (MQ=255)
                                                                                      tGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTc             >  1:3577774/1‑100 (MQ=255)
                                                                                      tGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTc             >  1:2645999/1‑100 (MQ=255)
                                                                                        cAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCt            >  2:3645664/1‑100 (MQ=255)
                                                                                        cAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCCTCTTCAACCGATTGAGGGTGCTGTCTg           >  2:3068863/1‑100 (MQ=255)
                                                                                           ttAGAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTg           <  2:1195199/97‑1 (MQ=255)
                                                                                            tataGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGt        <  1:2873055/100‑1 (MQ=255)
                                                                                             ataGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTtc      <  2:885241/100‑1 (MQ=255)
                                                                                             ataGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTtc      >  2:966981/1‑100 (MQ=255)
                                                                                              taGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTtct     >  1:1846436/1‑100 (MQ=255)
                                                                                                  tGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCtt  >  1:1414958/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTT  >  NC_000913/4122262‑4122459

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