Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 870,233 G→T 7.6% R65S (CGC→AGC)  PP_5456 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947870,2330GT 7.6% 531.1 / 20.9 163R65S (CGC→AGC) PP_5456hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (83/67);  new base T (6/7);  total (89/74)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

ATCGGCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGC  >  NC_002947/870084‑870380
                                                                                                                                                     |                                                                                                                                                   
aTCGGCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGcgcg                                                                                                                                                     >  1:842354/1‑150 (MQ=255)
aTCGGCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGcgcg                                                                                                                                                     >  1:842266/1‑150 (MQ=255)
 tCGGCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGa                                                                                                                                                    <  1:3432054/150‑1 (MQ=255)
      gAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCa                                                                                                                                               >  1:2403300/1‑150 (MQ=255)
        gTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCTCGCGCaaa                                                                                                                                             <  1:1090150/150‑1 (MQ=255)
        gTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCTCGCGCaaa                                                                                                                                             <  2:4821036/150‑1 (MQ=255)
        gTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCaaa                                                                                                                                             <  2:3134433/150‑1 (MQ=255)
        gTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCaaa                                                                                                                                             <  1:4540112/150‑1 (MQ=255)
         tGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAAc                                                                                                                                            <  2:2452718/150‑1 (MQ=255)
          gCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCTCGCGCAAAcc                                                                                                                                           <  1:3199246/150‑1 (MQ=255)
            ggggCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCg                                                                                                                                         >  1:106276/1‑150 (MQ=255)
            ggggCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCg                                                                                                                                         <  2:299242/150‑1 (MQ=255)
                              gTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTc                                                                                                                       >  2:4945657/1‑150 (MQ=255)
                              gTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTc                                                                                                                       >  1:3601095/1‑150 (MQ=255)
                               tCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTcg                                                                                                                      >  2:5113440/1‑150 (MQ=255)
                                 cTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTcgcg                                                                                                                    >  1:2228285/1‑150 (MQ=255)
                                  ttATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCgcgc                                                                                                                   >  2:2848223/1‑150 (MQ=255)
                                  ttATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCgcgc                                                                                                                   <  1:1060297/150‑1 (MQ=255)
                                  ttATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCgcgc                                                                                                                   >  1:3851429/1‑150 (MQ=255)
                                  ttATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCgcgc                                                                                                                   <  2:4386775/150‑1 (MQ=255)
                                  ttATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCgcgc                                                                                                                   <  2:4386516/150‑1 (MQ=255)
                                  ttATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCgcgc                                                                                                                   >  1:3851408/1‑150 (MQ=255)
                                   tataATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGcc                                                                                                                  >  1:3012041/1‑150 (MQ=255)
                                        tCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCccagc                                                                                                             >  2:4599096/1‑150 (MQ=255)
                                        tCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCccagc                                                                                                             >  2:4095726/1‑150 (MQ=255)
                                        tCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCccagc                                                                                                             >  1:3783109/1‑150 (MQ=255)
                                        tCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGAGGTTTGCGCGAGCGCAAACCCGGGCATCGGCGGGCAGCTCGCGCCccagc                                                                                                             >  1:4253621/1‑150 (MQ=255)
                                             gggAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAgccaggc                                                                                                        >  1:4465611/1‑150 (MQ=255)
                                                 gaTTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCATCCAGGCCAGc                                                                                                    <  2:3645723/150‑1 (MQ=255)
                                                  aTTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCt                                                                                                   <  2:2228285/150‑1 (MQ=255)
                                                      taGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAAt                                                                                               >  1:1373843/1‑150 (MQ=255)
                                                        gCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAAttt                                                                                             >  1:171874/1‑150 (MQ=255)
                                                           aTGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTGGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGc                                                                                          <  2:3711642/150‑1 (MQ=255)
                                                            tGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCt                                                                                         <  2:2354103/150‑1 (MQ=255)
                                                              tCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGa                                                                                       >  2:2962737/1‑150 (MQ=255)
                                                               cgcgGGCAGAAATGTCGCCAGTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGAt                                                                                      <  2:4562635/150‑1 (MQ=255)
                                                                gcgGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATc                                                                                     >  2:4338443/1‑150 (MQ=255)
                                                                     cAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCGGCTGAATTTCGCTGATCTGCAc                                                                                >  1:4852654/1‑150 (MQ=255)
                                                                     cAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCAc                                                                                <  2:3163249/150‑1 (MQ=255)
                                                                     cAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCAc                                                                                <  2:3163307/150‑1 (MQ=255)
                                                                           tGTCGCCATTCGTGGATTATATAGCACCGTTTCTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTc                                                                          <  2:4467763/150‑1 (MQ=255)
                                                                             tCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGa                                                                        >  1:4309703/1‑150 (MQ=255)
                                                                              cGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGAt                                                                       <  2:3020554/150‑1 (MQ=255)
                                                                              cGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGAt                                                                       <  1:148873/150‑1 (MQ=255)
                                                                               gCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTTATCTGCACCCAGTCGATg                                                                      >  2:2123681/1‑150 (MQ=255)
                                                                               gCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATg                                                                      <  1:4252793/150‑1 (MQ=255)
                                                                                 cATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGcc                                                                    >  2:127869/1‑150 (MQ=255)
                                                                                  aTTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGccc                                                                   >  2:4083446/1‑150 (MQ=255)
                                                                                    tCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCg                                                                 <  1:1569528/150‑1 (MQ=255)
                                                                                      gTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtg                                                               >  2:3872962/1‑150 (MQ=255)
                                                                                      gTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtg                                                               <  1:265912/150‑1 (MQ=255)
                                                                                        ggTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtggt                                                             <  2:3647206/150‑1 (MQ=255)
                                                                                        ggTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtggt                                                             <  2:3976842/150‑1 (MQ=255)
                                                                                           ttttATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGcc                                                          >  1:2737519/1‑150 (MQ=255)
                                                                                            tttATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCt                                                         >  1:712702/1‑150 (MQ=255)
                                                                                             ttATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTg                                                        >  1:2429033/1‑150 (MQ=255)
                                                                                                    aGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCAt                                                 >  2:2014778/1‑150 (MQ=255)
                                                                                                       gTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGc                                              >  2:728526/1‑150 (MQ=255)
                                                                                                        tttGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCt                                             <  2:3386649/150‑1 (MQ=255)
                                                                                                         ttGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCtt                                            <  2:171874/150‑1 (MQ=255)
                                                                                                            tCAGGCAGCGAGGAGGCGCTGGTGCCACCGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGaa                                         <  1:1898071/150‑1 (MQ=255)
                                                                                                            tCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGaa                                         <  2:941411/150‑1 (MQ=255)
                                                                                                            tCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGaa                                         <  2:1043331/150‑1 (MQ=255)
                                                                                                             cAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAg                                        <  2:3491786/150‑1 (MQ=255)
                                                                                                             cAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAg                                        <  2:3819712/150‑1 (MQ=255)
                                                                                                             cAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAg                                        <  2:4767607/150‑1 (MQ=255)
                                                                                                              aGGCAGCGCGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCCGCCAGGCCAGCTGAATTTCGCTGAGCTGCACCCAGTCGATGCCCCGGGGTGCCTGGCTGCATTGCTTGAAgg                                       >  2:5036560/1‑150 (MQ=255)
                                                                                                              aGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAgg                                       >  1:2760555/1‑150 (MQ=255)
                                                                                                              aGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAgg                                       >  2:674965/1‑150 (MQ=255)
                                                                                                                gCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGcc                                     >  1:2571249/1‑150 (MQ=255)
                                                                                                                 cAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCa                                    >  1:2815325/1‑150 (MQ=255)
                                                                                                                  aGCGAGGAGGCGCTGGGGCCACCGGCGGGTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAg                                   <  1:4157495/150‑1 (MQ=255)
                                                                                                                      aggaggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGTGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCAt                               >  1:3069568/1‑150 (MQ=255)
                                                                                                                           gcgcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGcc                          >  1:500124/1‑150 (MQ=255)
                                                                                                                                  tGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCcgcg                   <  1:1655215/150‑1 (MQ=255)
                                                                                                                                   gCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGt                  <  2:3783109/150‑1 (MQ=255)
                                                                                                                                    ccAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTc                 >  2:2164604/1‑150 (MQ=255)
                                                                                                                                      aaCGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCaa               >  1:3698088/1‑150 (MQ=255)
                                                                                                                                         ggcggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTg                                           >  1:1966891/1‑119 (MQ=255)
                                                                                                                                         ggcggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTg                                           >  1:1966668/1‑119 (MQ=255)
                                                                                                                                         ggcggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTg                                           <  2:1966668/119‑1 (MQ=255)
                                                                                                                                         ggcggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGc            >  2:2184651/1‑150 (MQ=255)
                                                                                                                                         ccgggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGc            >  1:4590180/4‑150 (MQ=255)
                                                                                                                                         ccgggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGc            >  2:4741908/4‑150 (MQ=255)
                                                                                                                                          cgggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCg           >  2:1776074/3‑150 (MQ=255)
                                                                                                                                           gggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg          <  1:4741908/149‑1 (MQ=255)
                                                                                                                                           cggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg          <  1:162694/150‑1 (MQ=255)
                                                                                                                                           cggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg          <  1:4095726/150‑1 (MQ=255)
                                                                                                                                           cggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg          <  2:2737519/150‑1 (MQ=255)
                                                                                                                                           cggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg          <  1:4902712/150‑1 (MQ=255)
                                                                                                                                            ggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCGGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCggg         >  2:120981/1‑150 (MQ=255)
                                                                                                                                            ggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCggg         <  2:4409636/150‑1 (MQ=255)
                                                                                                                                            ggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCggg         <  2:2863088/150‑1 (MQ=255)
                                                                                                                                              tttGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCa       >  2:5022343/1‑150 (MQ=255)
                                                                                                                                              tttGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCa       >  1:5049995/1‑150 (MQ=255)
                                                                                                                                                 gcgcGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATaa    <  2:106276/150‑1 (MQ=255)
                                                                                                                                                 gcgcGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATaa    <  2:842354/150‑1 (MQ=255)
                                                                                                                                                   gcgAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGCATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATGGCGTGAAGGCCAGGCATTTGCCATCCGCGTCACGGCGGGCATAAGc  >  2:2813906/1‑150 (MQ=255)
                                                                                                                                                   gcgAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGc  >  2:4657422/1‑150 (MQ=255)
                                                                                                                                                   gcgAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGc  >  1:4404547/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATCGGCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGC  >  NC_002947/870084‑870380

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