Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 870,234 A→C 6.3% A64A (GCT→GCG PP_5456 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947870,2340AC 6.3% 651.4 / 23.0 208A64A (GCT→GCGPP_5456hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (107/87);  new base C (9/5);  total (116/92)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

TCGGCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGCGTT  >  NC_002947/870085‑870383
                                                                                                                                                     |                                                                                                                                                     
tCGGCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGa                                                                                                                                                       <  1:1019352/150‑1 (MQ=255)
  ggCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAgc                                                                                                                                                     <  2:4863358/150‑1 (MQ=255)
     gAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCa                                                                                                                                                  >  1:2859651/1‑150 (MQ=255)
     gAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCa                                                                                                                                                  >  2:3520999/1‑150 (MQ=255)
        tGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATTTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAAc                                                                                                                                               <  2:4231891/150‑1 (MQ=255)
         gCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAAcc                                                                                                                                              >  2:3431236/1‑150 (MQ=255)
           ggggCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCg                                                                                                                                            >  1:3182588/1‑150 (MQ=255)
            gggCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGACATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCgg                                                                                                                                           >  2:5170858/1‑150 (MQ=255)
            gggCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCTCGCGCAAAccgcc                                                                                                                                           <  2:4237126/150‑4 (MQ=255)
              gCCCGATCAGGAGACGTCCTTATAATGGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCTGGCAGCGAGGTGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGc                                                                                                                                         <  1:2502885/150‑1 (MQ=255)
              gCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGc                                                                                                                                         >  2:3483271/1‑150 (MQ=255)
                   aTCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTcgg                                                                                                                                    >  1:3004945/1‑150 (MQ=255)
                       ggAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCggg                                                                                                                                <  2:1121046/150‑1 (MQ=255)
                       ggAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCggg                                                                                                                                <  2:1435748/150‑1 (MQ=255)
                        gagaCGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGc                                                                                                                               >  2:5393188/1‑150 (MQ=255)
                         agaCGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCa                                                                                                                              >  2:1676322/1‑150 (MQ=255)
                             gTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTc                                                                                                                          >  1:4277264/1‑150 (MQ=255)
                                cTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTcgcg                                                                                                                       >  1:4461501/1‑150 (MQ=255)
                                 ttATAATCGGGGGGAGATTATAGCGATGTCGCTGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGATTTAGTCAGGCAGCGAGGAGGCGATGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGAGGGCAGCTCgcgc                                                                                                                      >  2:1918341/1‑150 (MQ=255)
                                 ttATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCgcgc                                                                                                                      >  2:102963/1‑150 (MQ=255)
                                  tataATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGcc                                                                                                                     <  2:2414660/150‑1 (MQ=255)
                                  tataATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGcc                                                                                                                     >  1:2097669/1‑150 (MQ=255)
                                      aTCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCccag                                                                                                                 >  1:1049177/1‑150 (MQ=255)
                                       tCGGGGGGAGATTATATCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTATTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCacagc                                                                                                                >  2:897924/1‑150 (MQ=255)
                                       tCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCccagc                                                                                                                <  1:2378543/150‑1 (MQ=255)
                                       tCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCccagc                                                                                                                >  2:3610927/1‑150 (MQ=255)
                                        cGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCCGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCccagcc                                                                                                               <  2:3994301/150‑1 (MQ=255)
                                          gggggAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAg                                                                                                                                                      <  2:912878/109‑1 (MQ=255)
                                              gagaTTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAgccaggcca                                                                                                         >  2:2138048/1‑150 (MQ=255)
                                              gagaTTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAgccaggcca                                                                                                         <  1:591365/150‑1 (MQ=255)
                                               agaTTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAgccaggccag                                                                                                        <  2:726541/150‑1 (MQ=255)
                                                gaTTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCACCGGCGGTGTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGc                                                                                                       <  1:522080/150‑1 (MQ=255)
                                                 aTTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCt                                                                                                      <  2:2097669/150‑1 (MQ=255)
                                                  ttATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGc                                                                                                                                         >  2:2821268/1‑114 (MQ=255)
                                                  ttATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTg                                                                                                     >  2:3247535/1‑150 (MQ=255)
                                                     taGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAAt                                                                                                  <  1:2090304/150‑1 (MQ=255)
                                                      aGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTcggcgg                                                                                                                                 >  1:1647889/1‑118 (MQ=255)
                                                      aGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTATGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCgg                                                                                                                                           >  1:673104/1‑108 (MQ=255)
                                                          aTGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGc                                                                                             >  2:778550/1‑150 (MQ=255)
                                                           tGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCt                                                                                            <  1:4951020/150‑1 (MQ=255)
                                                             tCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGa                                                                                          >  1:903350/1‑150 (MQ=255)
                                                                cgGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCt                                                                                       >  2:2433156/1‑150 (MQ=255)
                                                                  ggCAGAAATGTCGCCATTGGTGGTTTTTATTGCAGGGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCTACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGc                                                                                     <  1:3842526/150‑1 (MQ=255)
                                                                     aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCt                                                                                            >  1:937436/1‑140 (MQ=255)
                                                                       aaaTGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCa                                                                                >  1:5331313/1‑150 (MQ=255)
                                                                        aaTGTCGCCATTCGTGGTTTTTATTGCATCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGGGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAg                                                                               <  2:889161/150‑1 (MQ=255)
                                                                        aaTGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAg                                                                               >  1:2820555/1‑150 (MQ=255)
                                                                         aTGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGACAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCACGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGt                                                                              >  1:3594127/1‑150 (MQ=255)
                                                                          tGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTc                                                                             <  2:4856841/150‑1 (MQ=255)
                                                                             cGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGAt                                                                          <  2:3050340/150‑1 (MQ=255)
                                                                              gCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATg                                                                         >  1:4729609/1‑150 (MQ=255)
                                                                              gCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATg                                                                         <  1:3149410/150‑1 (MQ=255)
                                                                              gCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATg                                                                         >  2:5206597/1‑150 (MQ=255)
                                                                               ccATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGATGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGc                                                                        >  1:4861209/1‑150 (MQ=255)
                                                                                  ttCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGcccc                                                                     >  1:5246125/1‑150 (MQ=255)
                                                                                     gTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtg                                                                  >  1:2311837/1‑150 (MQ=255)
                                                                                     gTGGTTTTTAGTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtg                                                                  <  2:4649370/150‑1 (MQ=255)
                                                                                        gTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtggtg                                                               >  2:1738084/1‑150 (MQ=255)
                                                                                         tttttATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGc                                                              >  2:3970164/1‑150 (MQ=255)
                                                                                          ttttATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGcc                                                             <  1:4087794/150‑1 (MQ=255)
                                                                                           tttATTGCAGCTTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGTTGCCt                                                            <  2:2131295/150‑1 (MQ=255)
                                                                                           tttATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCt                                                            >  2:3311920/1‑150 (MQ=255)
                                                                                               atGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTTGCt                                                        <  2:604392/149‑1 (MQ=255)
                                                                                                  ccgCGTTGGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGGTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCa                                                     <  2:4951990/148‑1 (MQ=255)
                                                                                                  cAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCa                                                     >  2:4699097/1‑150 (MQ=255)
                                                                                                  cAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCa                                                     >  1:181181/1‑150 (MQ=255)
                                                                                                   aGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCAt                                                    >  1:3938239/1‑150 (MQ=255)
                                                                                                   aGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCAt                                                    >  2:2042719/1‑150 (MQ=255)
                                                                                                   aGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCAt                                                    >  1:1431739/1‑150 (MQ=255)
                                                                                                      gTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGc                                                 <  2:5273892/150‑1 (MQ=255)
                                                                                                             aGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTcggcgg                                                                                                                                 <  2:4812931/63‑1 (MQ=255)
                                                                                                             aGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAgg                                          >  2:1084687/1‑150 (MQ=255)
                                                                                                             aGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAgg                                          >  2:3825312/1‑150 (MQ=255)
                                                                                                              ggCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGc                                         <  2:3002406/150‑1 (MQ=255)
                                                                                                              ggCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGc                                         <  2:2084504/150‑1 (MQ=255)
                                                                                                              ggCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGc                                         <  1:2138048/150‑1 (MQ=255)
                                                                                                               gCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGTGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGcc                                        >  2:3349961/1‑150 (MQ=255)
                                                                                                               gCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGcc                                        >  2:3710880/1‑150 (MQ=255)
                                                                                                                 aGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAg                                      <  2:2673647/150‑1 (MQ=255)
                                                                                                                 aGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAg                                      <  1:1952581/150‑1 (MQ=255)
                                                                                                                   cGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCgg                                                                                                                                           <  1:5235551/47‑1 (MQ=255)
                                                                                                                    gaggagGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCa                                   >  1:822626/1‑150 (MQ=255)
                                                                                                                     aggaggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGc                                    >  1:1794926/1‑148 (MQ=255)
                                                                                                                      ggaggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCAtt                                 >  2:2160889/1‑150 (MQ=255)
                                                                                                                          gcgcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGcc                             >  1:2064674/1‑150 (MQ=255)
                                                                                                                              tGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATcc                         <  1:2269680/150‑1 (MQ=255)
                                                                                                                              tGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATcc                         <  1:4657636/150‑1 (MQ=255)
                                                                                                                                   ccAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCTCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTc                    >  2:438798/1‑150 (MQ=255)
                                                                                                                                    cAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCa                   >  1:2374924/1‑150 (MQ=255)
                                                                                                                                          cggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg             <  1:102963/150‑1 (MQ=255)
                                                                                                                                          cggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg             <  2:1092813/150‑1 (MQ=255)
                                                                                                                                           ggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCggg            >  1:3135235/1‑150 (MQ=255)
                                                                                                                                           ggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCggg            <  1:2199015/150‑1 (MQ=255)
                                                                                                                                            gTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGc           >  1:1238223/1‑150 (MQ=255)
                                                                                                                                                gcgcGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATaa       >  2:1723050/1‑150 (MQ=255)
                                                                                                                                                gcgcGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTCGCCATCCGCGTCAAGGCGGGCATaa       <  1:3398644/150‑1 (MQ=255)
                                                                                                                                                  gcgAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCAGTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGc     >  2:4945760/1‑150 (MQ=255)
                                                                                                                                                  gcgAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTACATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGc     >  2:2719327/1‑150 (MQ=255)
                                                                                                                                                     aGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGCGtt  <  1:5388007/150‑1 (MQ=255)
                                                                                                                                                     aGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGCGtt  <  2:1865317/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGGCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGCGTT  >  NC_002947/870085‑870383

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