Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,121,597 (T)5→4 100% coding (151/1332 nt) hslU ← molecular chaperone and ATPase component of HslUV protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,121,5930T.97.2% 758.8 / 8.2 173coding (155/1332 nt)hslUmolecular chaperone and ATPase component of HslUV protease
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base . (87/81);  total (90/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

TTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCC  >  NC_000913/4121494‑4121688
                                                                                                   |                                                                                               
ttGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCGGATCATCAGGatat                                                                                                 >  2:2604341/1‑100 (MQ=255)
  ggTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttt                                                                                               >  1:1334860/1‑100 (MQ=255)
   gTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAtttt                                                                                              >  1:2093508/1‑100 (MQ=255)
     cGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcg                                                                                            <  1:1406023/100‑3 (MQ=255)
      gCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGCTATTTTCCg                                                                                           >  2:3107518/1‑100 (MQ=255)
      gCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGGCGGGCCGATGATCAGGATAttttcgg                                                                                           >  1:1944624/1‑97 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          <  1:1176435/100‑5 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          >  2:463605/1‑96 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          >  1:3211156/1‑96 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          <  1:2015206/100‑5 (MQ=255)
        ttCAACTTTGATGCCCGGCGCCTTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCgggg                                                                                        <  1:1127337/100‑1 (MQ=255)
        ttCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCgggg                                                                                        <  1:2382041/100‑1 (MQ=255)
         tccaCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcgg                                                                                           >  2:3615082/4‑94 (MQ=255)
         tCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGt                                                                                       >  1:1400203/1‑100 (MQ=255)
            aCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCAc                                                                                    <  1:2665068/100‑1 (MQ=255)
            aCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCAc                                                                                    >  2:372808/1‑100 (MQ=255)
                tgatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCa                                                                                >  2:1785146/1‑100 (MQ=255)
                tgatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCa                                                                                >  2:3700557/1‑100 (MQ=255)
                 gatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCAt                                                                               <  2:2425823/100‑1 (MQ=255)
                  atgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCCCTTCAt                                                                               >  1:440589/1‑99 (MQ=255)
                   tgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATgg                                                                             >  2:2869158/1‑100 (MQ=255)
                    gaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGgc                                                                            >  2:575485/1‑100 (MQ=255)
                    gaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGgc                                                                            >  1:1333878/1‑100 (MQ=255)
                     aaCGGCGCATTCGCCAGCTTAGCCAGACGACGGTCGAATTAAGTTTTACCGACACCGGTCGGGTCTATCATCAGGATATTTTCGGGGTCACTTCATgg                                                                             <  2:2688109/98‑1 (MQ=255)
                      aCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCCGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGgcgc                                                                          <  1:3043680/100‑1 (MQ=255)
                         gcgcATTCGCCAGCTTAGCCAGACGACGGGCGCTTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCTGGGTCACTTCATGGcgccgc                                                                       >  2:1059571/1‑100 (MQ=255)
                         gcgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGc                                                                       <  2:1467508/100‑1 (MQ=255)
                         gcgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGc                                                                       >  2:1080067/1‑100 (MQ=255)
                          cgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGc                                                                       <  2:1224274/99‑1 (MQ=255)
                           gcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGctc                                                                     >  2:946744/1‑100 (MQ=255)
                            cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGctct                                                                    >  1:1083650/1‑100 (MQ=255)
                            cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGctct                                                                    >  1:1033251/1‑100 (MQ=255)
                            cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTGTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGctct                                                                    >  2:2033638/1‑100 (MQ=255)
                             aTTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCtt                                                                   >  2:886786/1‑100 (MQ=255)
                               tCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCg                                                                 <  1:1778630/100‑1 (MQ=255)
                                cGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGCTCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGt                                                                >  1:3234742/1‑100 (MQ=255)
                                  ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                              <  1:186508/100‑1 (MQ=255)
                                  ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                              <  2:152931/100‑1 (MQ=255)
                                   cAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGa                                                             <  1:1547850/100‑1 (MQ=255)
                                   cAGCTTAGCCAGACGACGCGCTATTTCATTTTTCCCGACACCGGTCGGGCCGAGCATCAGGATATTTACGGGGTCACTTCATGGCGCAGCTCTTCGt                                                                <  2:3225751/97‑1 (MQ=255)
                                    aGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAg                                                            >  1:3196032/1‑100 (MQ=255)
                                    aGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAg                                                            >  1:1292860/1‑100 (MQ=255)
                                         aGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCa                                                       <  1:3328076/100‑1 (MQ=255)
                                          gCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCAt                                                      >  1:754738/1‑100 (MQ=255)
                                          cccAGACGACGGGCGCTTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATCGTTTCGGGGTCACTTCATGCCGCAGCTCTTCGTTCAGCTGcct                                                      >  1:1761774/2‑98 (MQ=255)
                                           ccAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATg                                                     <  2:2788040/100‑1 (MQ=255)
                                            cAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGc                                                    >  1:3095584/1‑100 (MQ=255)
                                            cAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGc                                                    >  1:2296523/1‑100 (MQ=255)
                                             aGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCg                                                   >  2:1448417/1‑100 (MQ=255)
                                                cgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACg                                                >  1:2486032/1‑100 (MQ=255)
                                                 gacgGGCGATTTCAGTTTTACCGGCACCGGTCGGGCCGATCATCAGGATGTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGAc                                                 <  2:361903/98‑1 (MQ=255)
                                                  acgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGcc                                              >  2:2494640/1‑100 (MQ=255)
                                                  acgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGc                                               <  2:1869380/99‑1 (MQ=255)
                                                   cgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCa                                             >  1:3408909/1‑100 (MQ=255)
                                                   cgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCa                                             >  2:938921/1‑100 (MQ=255)
                                                    gggCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAg                                            <  1:2853873/100‑1 (MQ=255)
                                                    gggCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAg                                            >  2:521318/1‑100 (MQ=255)
                                                     ggCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGCCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGc                                           >  1:675166/1‑100 (MQ=255)
                                                        gATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGt                                        <  2:3403080/100‑1 (MQ=255)
                                                          tttCAGTTTTGCCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTGCATGGCGCGGCTCTTCGTTGAGCTTCATGGGACGCCAGGGGTTa                                      >  1:3627734/1‑100 (MQ=255)
                                                             cAGTTTTACCGACACCGGTCGGGTCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACg                                    <  1:2986609/99‑1 (MQ=255)
                                                             cAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGc                                   <  2:1373767/100‑1 (MQ=255)
                                                              aGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                  <  1:3675327/100‑1 (MQ=255)
                                                               gTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGc                                   <  2:3599710/98‑1 (MQ=255)
                                                               gTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 >  1:1427368/1‑100 (MQ=255)
                                                               gTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 >  1:3340762/1‑100 (MQ=255)
                                                               gTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 >  1:3136457/1‑100 (MQ=255)
                                                                 tttACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                               <  1:902540/100‑1 (MQ=255)
                                                                 tttACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                               >  2:2381821/1‑100 (MQ=255)
                                                                  ttACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGc                              <  1:3410567/100‑1 (MQ=255)
                                                                  ttACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGc                              <  1:968863/100‑1 (MQ=255)
                                                                   tACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCa                             <  1:2712171/100‑1 (MQ=255)
                                                                    aCCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCaa                            >  1:3312472/1‑100 (MQ=255)
                                                                    aCCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCaa                            >  2:2893496/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  1:414622/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  2:1988447/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  1:498698/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  2:1201352/1‑100 (MQ=255)
                                                                      cGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATc                          >  2:2715248/1‑100 (MQ=255)
                                                                       gACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCg                         >  2:2256666/1‑100 (MQ=255)
                                                                         cacCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGcc                       >  2:381787/1‑100 (MQ=255)
                                                                         cacCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGcc                       >  1:3175449/1‑100 (MQ=255)
                                                                          acCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCca                      >  2:343478/1‑100 (MQ=255)
                                                                            cGGTCGGGCCGATCATCAGGATATTTGCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcaca                    <  2:2995072/100‑1 (MQ=255)
                                                                             ggTCGGGCAGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCACATCGTTGAGATGAACGCTACGCCATCGGTTACGCAGAGCAATCGCCCCAg                   >  1:1616668/1‑100 (MQ=255)
                                                                              gTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGa                  <  1:1047704/100‑1 (MQ=255)
                                                                               tCGGGCCGATCATCAGGTTATTTTTCGGGGTCACTTCATGGCGCTGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcac                     <  2:3185825/97‑1 (MQ=255)
                                                                                  ggCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGc              <  2:1794753/100‑1 (MQ=255)
                                                                                     cGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           <  2:1334851/100‑1 (MQ=255)
                                                                                     cGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           <  1:3511953/100‑1 (MQ=255)
                                                                                     cGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           <  1:209191/100‑1 (MQ=255)
                                                                                      gATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGCCGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           >  2:2712001/1‑99 (MQ=255)
                                                                                         catcaGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGt       >  1:3053430/1‑100 (MQ=255)
                                                                                          atcaGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGtt      >  2:75171/1‑100 (MQ=255)
                                                                                           tcaGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGCGCAATCGCCACAGAACGCTTcgcg        >  2:1267291/1‑97 (MQ=255)
                                                                                           tcaGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTg     <  1:2031550/100‑1 (MQ=255)
                                                                                           tcaGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTg     <  1:1042778/100‑1 (MQ=255)
                                                                                             aGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTc   <  1:847069/100‑1 (MQ=255)
                                                                                              ggATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc  >  1:60689/1‑100 (MQ=255)
                                                                                              ggATATTTTCGGGGTCACGTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc  >  2:1389466/1‑100 (MQ=255)
                                                                                                   |                                                                                               
TTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCC  >  NC_000913/4121494‑4121688

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: