Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,445,494 T→G 19.9% F763C (TTT→TGT)  PP_3062 → tail length determination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,445,4940TG19.9% 147.8 / 4.7 93F763C (TTT→TGT) PP_3062tail length determination protein
Reads supporting (aligned to +/- strand):  ref base T (55/16);  new base G (16/5);  total (71/22)
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

AGGGCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCGGTC  >  NC_002947/3445347‑3445642
                                                                                                                                                   |                                                                                                                                                    
aGGGCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTg                                                                                                                                                    <  2:4303899/150‑1 (MQ=255)
aGGGCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTg                                                                                                                                                    <  1:2736360/150‑1 (MQ=255)
 gggCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTgg                                                                                                                                                   <  1:4026712/150‑1 (MQ=255)
  ggCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTggg                                                                                                                                                  >  1:3480577/1‑150 (MQ=255)
   gCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTgggg                                                                                                                                                 <  1:4629522/150‑1 (MQ=255)
   gCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTgggg                                                                                                                                                 <  1:2735472/150‑1 (MQ=255)
    cTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTggggg                                                                                                                                                <  1:3983675/150‑1 (MQ=255)
       cggcggCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTg                                                                                                                                             >  2:3903878/1‑150 (MQ=255)
        ggcggcGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGGTTGGGGGGTgg                                                                                                                                            >  2:4989701/1‑150 (MQ=255)
         gcggcgCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGc                                                                                                                                           >  1:5037623/1‑150 (MQ=255)
          cggcgCTGCCGGCAGACTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGGGTGGGGGGTGGcg                                                                                                                                          >  2:4218822/1‑149 (MQ=255)
              gcTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCGGGGGGGGCGAGTCGGGTGGGGGGTGGCTGGGc                                                                                                                                      >  2:3908271/1‑150 (MQ=255)
              gcTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGc                                                                                                                                      >  1:2936251/1‑150 (MQ=255)
                tGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCaa                                                                                                                                    >  1:2218388/1‑150 (MQ=255)
                         ccTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGtt                                                                                                                           >  1:3483115/1‑150 (MQ=255)
                          cTGGCGGGCCCCCCCGCCGGGGCTGCCGCGGGGGCGGCTATGGGCTCCGTGGTGCCGGTGAGCGGACCCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGCGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGttt                                                                                                                          <  2:2196158/150‑1 (MQ=255)
                          cTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCACGCGATGGttt                                                                                                                          >  1:947158/1‑150 (MQ=255)
                           gggCGGGCCCCCGGGCCGGGGCGGCAGCGGGGGCGGCCAGGGGCCCCGGGGGGCCGGGGCCCGGGCCCGCGGGGGGGGGCGCGGGCGGCGCGGCCCCCGGTGGCAGGGGGGGCGAGCCGTGGGGGGGGGGGCTGGGCAAGCGATGGTTTg                                                                                                                         <  2:2570762/149‑1 (MQ=255)
                           cggCGGGCCCCCTCGCCGGTGCGGCAGCGGGGGCGGCTCCTGGCCCCGTGGAGCCGGTGACCGCAACCGCCGGGGGGGGCGCGGCCGGCGCGGCCCCGGGCGCCAGGGGGGGCGCGCCGTTTGGGGGGGGGCTGGGCAAGCGATGGTTTg                                                                                                                         <  2:2218388/149‑1 (MQ=255)
                            cgCGGGCACCCGGGCCGGGGCTGCAGCGGGGGCGGCTATTGGCTCCGCGGTGCCGGTGATCGGTACCGCTGTGGGCGGCGCGGCCGGCGCTGACCTTGGCGGCACGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGTGATGGTTTgg                                                                                                                        <  2:1331770/149‑1 (MQ=255)
                               gggCCCCCCGGCCGGGGCGGCCGCGGGGGCGGCGAGGGGCTCCGTGGGGCCGGGGCGCGGGCCCGCGGGGGGGGGCGCGGGCGGCGCGGCCCGGGGGGGCCGGGGGGGCGCGGCGCGTGGGGGGGGGCGGGGCAAGCGATGGTTTGGTGa                                                                                                                     <  2:2902273/150‑1 (MQ=255)
                               gggCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCGTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGa                                                                                                                     >  2:4526931/1‑150 (MQ=255)
                                  cACCCGTGCCGGCGCCGCAGCGCGGGCGCCTCGGGGCTCCGTGGCGCCGCTGACCGCCCCCGCCGCGGGCGGCGCGGGCGCCGCCGCCCCCGGCGCCAAGGGCGGCGAGCCGGGGGGGGGGGGGCTGGGCAAGCGATGGTTTGGTGACGa                                                                                                                  <  2:4120751/150‑1 (MQ=255)
                                  cACCCCCGCCGGCGCCGCCGCGGGGGCGGCGCCGGGCCCCCTGGCGCCGGCGACCGGGACCGCCGCGCGCGGCGCGCCCGCCGCGGCCCCCGCCCCCACGGGGGGCGAGCCGCTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGa                                                                                                                  <  2:2932839/150‑1 (MQ=255)
                                    ccccggCCGGGGCGGCCGCGGGGGCGGCGAGGGGCGCCGGGCGGCCGGTGCGCGGCACCGCGGGGGGCGGCGCGCGCGCCGCCGTCCTCGGGGGCCAGGGGGGCGCGCCGCCTGGGGGGGGGCTGGGCAAGCGATGGTTTGGTGACGAgc                                                                                                                <  2:3862340/145‑1 (MQ=255)
                                     ccTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTggg                                                                                                                                       >  1:888802/1‑126 (MQ=255)
                                     ccTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTggg                                                                                                                                       <  2:888802/126‑1 (MQ=255)
                                     ccTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGTCGAgcc                                                                                                               >  1:667579/1‑150 (MQ=255)
                                     ccTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAgcc                                                                                                               >  2:2076546/1‑150 (MQ=255)
                                     ccTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAgcc                                                                                                               >  1:4177361/1‑150 (MQ=255)
                                     ccTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAgcc                                                                                                               >  2:2917518/1‑150 (MQ=255)
                                     ccTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAgcc                                                                                                               >  1:1541553/1‑150 (MQ=255)
                                       ttGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCAGGGGGGGCGAGTCGGGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAgccgg                                                                                                             >  2:4270995/1‑150 (MQ=255)
                                        cgCCGGGGCGGCAGCGGGGGCGGCCACGGGCGCCGTGGGGCCGGGGCGCGGGACCGCGGTGGGTGGCGCGGTCGGCGCGGCCCCTGGGGGCCTGGGGGGCGCGCCGTGTGGGGGGCGGCTGGGCAAGCGATGGTTTGGTGACGAgccggc                                                                                                            <  2:4409296/149‑1 (MQ=255)
                                                 ggccgCGGGGGCGGCCAGGGGCGCCGGGGTGCCGGTGAGCGGGCCCGCGGGGGGGGGCGCGGTCGCCGCGGTCCGTGGTGGCCGGGGGGCCGCGCCGTGTGGGGGGGGGCTGGGCAAGCAATGGTTTGGTGACGAGCCGGCCGACGATCa                                                                                                   <  2:4923373/146‑1 (MQ=255)
                                                   cAGCGGGGGCGGCGATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATAGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAag                                                                                                 <  2:4542579/150‑1 (MQ=255)
                                                             ggcTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGTTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGa                                                                                       >  2:2995017/1‑150 (MQ=255)
                                                             ggcTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGGGGCTGGGCCAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGa                                                                                       >  1:256532/1‑150 (MQ=255)
                                                               cTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAAt                                                                                     >  1:3924331/1‑150 (MQ=255)
                                                                tATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGCATc                                                                                    >  1:4298294/1‑150 (MQ=255)
                                                                tATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATc                                                                                    >  2:5302181/1‑150 (MQ=255)
                                                                tATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATc                                                                                    >  2:5302077/1‑150 (MQ=255)
                                                                tATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATc                                                                                    >  2:2697057/1‑150 (MQ=255)
                                                                tATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATAAAGCCGAGCCGGaaac                                                                                    >  1:199146/1‑148 (MQ=255)
                                                                tATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGACGAGCCGGAATc                                                                                    >  2:5302032/1‑150 (MQ=255)
                                                                     gCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCCAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATcgccgc                                                                               >  1:1652832/1‑150 (MQ=255)
                                                                      ttgaccggtgCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCACGCGATGGGTTGGTGGCGAGCCGGGCGACGATCAAGCCGAGCCGGACTCgccgcc                                                                              >  2:2848210/7‑150 (MQ=255)
                                                                       tCCGTGGTGCCGGTGATCGGTACCGCTGTGGTTGGCTCGGTCGGCGCTGTCCTTGGTGGCATGGTGGGCTAGTCGTTTGGGGGGTGGATGGGCAAGCGATTGTTTGGTGACGAGCCGGCCGACGATAAAGCCGAGCCGGAATCGccgccg                                                                             >  2:2551950/1‑150 (MQ=255)
                                                                       tCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGccgccg                                                                             >  1:2885879/1‑150 (MQ=255)
                                                                       tCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCCAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGccgccg                                                                             >  2:2790532/1‑150 (MQ=255)
                                                                       tCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCAGGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCCAGCGATGGGTTGGTGACGCGCCGGCCGACGATCAAGCCGAGCCGGAATCGccgccg                                                                             >  2:2823853/1‑150 (MQ=255)
                                                                              tgCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCCAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTgg                                                                      >  2:533351/1‑150 (MQ=255)
                                                                                ccGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGt                                                                    >  1:5373549/1‑150 (MQ=255)
                                                                                   gTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTaaa                                                                 >  1:5078974/1‑150 (MQ=255)
                                                                                   gTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCCGAATCGCCCCCGGCACTGGGTaaa                                                                 >  1:224211/1‑150 (MQ=255)
                                                                                          gTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCCAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCCCGGGGTAAAGCGTTAc                                                          >  1:2313878/1‑150 (MQ=255)
                                                                                          gTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGCTGGGTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACAGGGTAAAGCGGTAc                                                          >  1:557443/1‑150 (MQ=255)
                                                                                           tACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTTGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGTTCAAGCCGTGCCGGAATCGCCGCCGGCAGTGGGTAATGCGTTagg                                                         >  2:2193354/1‑148 (MQ=255)
                                                                                           tACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCCAGCGATGGTTTGGTGACGAGCCGGCCGACGGTCAAGCCGAGCCGGAATCGCCGCCGGCACGGGGTAAAGCGTTAcg                                                         >  2:83616/1‑150 (MQ=255)
                                                                                           tACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTAcg                                                         >  2:3564875/1‑150 (MQ=255)
                                                                                           tACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTAcg                                                         >  1:2895651/1‑150 (MQ=255)
                                                                                             ccGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTAcgcg                                                       >  2:5338985/1‑150 (MQ=255)
                                                                                              cGCTGGGGGCGGCGCGGTCGGCGCTGCCCCTGGCGGCCCGGGGGGCGCGCCGCTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGt                                                      <  2:4991268/150‑1 (MQ=255)
                                                                                               gCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCCAGCGATGGTTTGGTGACGAGCCGGCCGACCATCAAGCCGAGCCGGAATCCCCGCCGGCACAGGGTAAAGCGTTACGCGTa                                                     >  1:491083/1‑150 (MQ=255)
                                                                                               gCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGCTGGTTTGGTGACGAGCCGGCCGACGAACAAGCCCAGCCGGAATCGCCCCCGGCACAGGGTAAAGCGTTACGCGTa                                                     >  1:490851/1‑150 (MQ=255)
                                                                                                cTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTGGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTaa                                                    >  2:934227/1‑150 (MQ=255)
                                                                                                        ggCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATg                                            >  2:2216462/1‑150 (MQ=255)
                                                                                                         gcgcGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGGAACGGTGATgc                                           >  2:3872832/1‑150 (MQ=255)
                                                                                                         gcgcGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATgc                                           >  2:3872739/1‑150 (MQ=255)
                                                                                                              gTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATgccagcc                                      >  2:776347/1‑150 (MQ=255)
                                                                                                               acGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATgccagcca                                     >  2:35549/2‑150 (MQ=255)
                                                                                                                      tGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATTGTTTGGTGACGAGCCCGCCGACGATAAAGCCGAGCCGGAATCGCCCCCCGCACAGGGTAAAGCGTTAAGCGTAACCGTGATGCCAGCCAAAGAAGt                              >  1:5056096/1‑150 (MQ=255)
                                                                                                                       gTCCTTGGTGGCATGGGGGGCGAGTCGTGTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAATCCGAGCCGGAATCGCCGCCCGCACAGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTc                             >  2:1341181/1‑150 (MQ=255)
                                                                                                                          gttGGTGGCATGGGGGGCGTGTCGGTTGGGGGGTGGCTGGGCACGCGCTGGTTTGGTGTCGAGCCCGCCGGCGTTCAAGCCGGGCCGGAATCGCCGGCGGCCCTGGGTAATGCGTTACGGGTAACGGTGATGCCAGCCAAGGAAGTCgtg                          >  2:1546713/2‑148 (MQ=255)
                                                                                                                          cTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAg                          >  1:4508126/1‑150 (MQ=255)
                                                                                                                          cTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAAACGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAg                          >  1:2337822/1‑150 (MQ=255)
                                                                                                                          cTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCCAGCCGGAATCGCCGCCGGCACGGGGTAAAGCGATACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAg                          >  2:1546560/1‑150 (MQ=255)
                                                                                                                           atggtggCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGc                         >  2:768952/2‑150 (MQ=255)
                                                                                                                             ggtggCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGAt                                             >  1:3981260/1‑128 (MQ=255)
                                                                                                                             ggtggCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGAt                                             <  2:3981260/128‑1 (MQ=255)
                                                                                                                             ggtggCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCt                       >  2:4014692/1‑150 (MQ=255)
                                                                                                                             ggtggCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCt                       >  2:4014340/1‑150 (MQ=255)
                                                                                                                                ggCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGATCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAg                    >  1:2668288/1‑150 (MQ=255)
                                                                                                                                ggCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAg                    >  1:2668270/1‑150 (MQ=255)
                                                                                                                                   aTGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCCGCCGACGATCAAGCCGAGCCGGAATCGCCCCCGGCACAGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGccc                 >  1:1717773/1‑150 (MQ=255)
                                                                                                                                          gCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACgg          <  2:4508126/150‑1 (MQ=255)
                                                                                                                                          gCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACgg          >  2:2046927/1‑150 (MQ=255)
                                                                                                                                           cGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAATCCCAGCCGGAAACGCCGCCCGCACAGGGTAAATCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCAGAGCCCGCCACGGc         >  1:3516688/1‑150 (MQ=255)
                                                                                                                                            gAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCa                                                                                                   <  1:2602867/59‑1 (MQ=255)
                                                                                                                                            gAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCa                                                                                                   >  2:2602867/1‑59 (MQ=255)
                                                                                                                                            gAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGAACAAGCCGAGCCCGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGca        >  2:2506805/1‑149 (MQ=255)
                                                                                                                                               tCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCg     >  1:791673/1‑150 (MQ=255)
                                                                                                                                               tCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCg     >  2:4002537/1‑150 (MQ=255)
                                                                                                                                                  tttGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCGGTc  >  2:5186970/1‑150 (MQ=255)
                                                                                                                                                  tttGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCGGTc  >  2:539579/1‑150 (MQ=255)
                                                                                                                                                  tttGGGGGGTGGCTGGGCAAGCGAGGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCGGTc  >  2:3613517/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGGGCTACGGCGGCGCTGCCGGCAGCCTGGCGGGCACCCTTGCCGGTGCTGCAGCGGGGGCGGCTATTGGCTCCGTGGTGCCGGTGATCGGTACCGCTGTGGGTGGCGCGGTCGGCGCTGTCCTTGGTGGCATGGGGGGCGAGTCGTTTGGGGGGTGGCTGGGCAAGCGATGGTTTGGTGACGAGCCGGCCGACGATCAAGCCGAGCCGGAATCGCCGCCGGCACTGGGTAAAGCGTTACGCGTAACGGTGATGCCAGCCAAGGAAGTCGAGCCTGAGCCCGCCACGGCTCCGGTC  >  NC_002947/3445347‑3445642

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