Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,003,786 G→A 8.0% R182Q (CGG→CAG)  PP_2625 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,003,7860GA 8.0% 573.3 / 27.4 187R182Q (CGG→CAG) PP_2625hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (93/79);  new base A (5/10);  total (98/89)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 213 total aligned reads displayed.

TTGCATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTG  >  NC_002947/3003637‑3003933
                                                                                                                                                     |                                                                                                                                                   
ttGCATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCcgggcg                                                                                                                                                     >  2:2600119/1‑150 (MQ=255)
  gCATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCcgggcggg                                                                                                                                                   <  2:1493827/150‑1 (MQ=255)
   cATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCgggg                                                                                                                                                  <  1:1871081/150‑1 (MQ=255)
   cATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCgggg                                                                                                                                                  >  2:2717166/1‑150 (MQ=255)
    aTACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGc                                                                                                                                                 >  2:1732810/1‑150 (MQ=255)
            cATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGc                                                                                                                                         >  1:4221934/1‑150 (MQ=255)
             aTCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGcc                                                                                                                                        <  2:2625995/150‑1 (MQ=255)
             aTCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCAGGCcccgcccggc                                                                                                                                        <  1:2634745/150‑3 (MQ=255)
              tctcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCa                                                                                                                                       <  2:3154313/150‑1 (MQ=255)
                tcCCGGCTACAGCTTATTTCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTCCGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCCAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATc                                                                                                                                     <  2:828095/150‑1 (MQ=255)
                tcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                                     >  2:671300/1‑150 (MQ=255)
                tcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                                     <  1:2252940/150‑1 (MQ=255)
                tcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                                     <  1:1994469/150‑1 (MQ=255)
                tcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                                     <  2:700138/150‑1 (MQ=255)
                tcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                                     >  2:151969/1‑150 (MQ=255)
                  ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATCGa                                                                                                                                   <  1:1541819/150‑1 (MQ=255)
                  ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGa                                                                                                                                   <  1:393652/150‑1 (MQ=255)
                   cGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGAGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAg                                                                                                                                  >  2:3505470/1‑150 (MQ=255)
                    ggCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATCGAGc                                                                                                                                 <  2:2367883/150‑1 (MQ=39)
                       tACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGa                                                                                                                              >  1:2862352/1‑150 (MQ=255)
                       tACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGa                                                                                                                              >  2:1723830/1‑150 (MQ=255)
                         cAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATc                                                                                                                            >  1:3560544/1‑150 (MQ=255)
                          aGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCa                                                                                                                           <  1:3755060/150‑1 (MQ=255)
                            cTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATCGAGCacttggcc                                                                                                                         <  2:556865/150‑8 (MQ=38)
                            cTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAgg                                                                                                                         <  1:834810/150‑1 (MQ=255)
                              gATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGaa                                                                                                                       <  2:3536175/150‑1 (MQ=255)
                                  gCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTg                                                                                                                   <  2:4245228/150‑1 (MQ=255)
                                   ccAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGa                                                                                                                  >  1:1525205/1‑150 (MQ=255)
                                    cAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGAt                                                                                                                 <  1:2680723/150‑1 (MQ=255)
                                    cAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGAt                                                                                                                 >  2:2650755/1‑150 (MQ=255)
                                      ggCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCg                                                                                                               <  2:3718218/150‑1 (MQ=255)
                                           ccGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGcc                                                                                                          <  2:3878792/150‑1 (MQ=255)
                                            cGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGccc                                                                                                         <  2:3191732/150‑1 (MQ=255)
                                            cGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGccc                                                                                                         >  1:1444293/1‑150 (MQ=255)
                                                  aGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTa                                                                                                   >  1:2814398/1‑150 (MQ=255)
                                                  aGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTa                                                                                                   >  2:4259220/1‑150 (MQ=255)
                                                     cAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGTGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCg                                                                                                >  2:2868116/1‑150 (MQ=255)
                                                     cAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGTGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCg                                                                                                >  2:2017875/1‑150 (MQ=255)
                                                       aGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAg                                                                                              <  1:2772155/150‑1 (MQ=255)
                                                         cTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGt                                                                                            <  1:1209643/150‑1 (MQ=255)
                                                          tGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGtt                                                                                           <  2:3091529/150‑1 (MQ=255)
                                                           gggTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTtg                                                                                          <  1:3227776/150‑1 (MQ=255)
                                                            ggTCACGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTtgt                                                                                         <  2:3836370/150‑1 (MQ=255)
                                                             gTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTtgtg                                                                                        >  2:3003907/1‑150 (MQ=255)
                                                              tCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtga                                                                                       <  1:1141877/150‑1 (MQ=255)
                                                              tCAAGGTAAGACTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCTATGCCGGGCGGTGCCTGCCCGCCATCGAGCAGAACAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtga                                                                                       >  2:418054/1‑150 (MQ=255)
                                                                   gTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGAt                                                                                  <  1:239344/150‑1 (MQ=255)
                                                                    tACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATc                                                                                 >  2:4367823/1‑150 (MQ=255)
                                                                         cTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGc                                                                            <  2:2246750/150‑1 (MQ=255)
                                                                          tGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCg                                                                           >  2:2905728/1‑150 (MQ=255)
                                                                             gCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAgg                                                                        <  2:2975348/150‑1 (MQ=255)
                                                                               ggggAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTc                                                                      >  1:1894120/1‑150 (MQ=255)
                                                                                 ggAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGa                                                                    >  1:1486996/1‑150 (MQ=255)
                                                                                 ggAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGa                                                                    >  1:4040146/1‑150 (MQ=255)
                                                                                 ggAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGa                                                                    >  1:254056/1‑150 (MQ=255)
                                                                                  gAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGAc                                                                   >  2:2692749/1‑150 (MQ=255)
                                                                                   cgggCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACg                                                                  <  1:3696360/149‑1 (MQ=255)
                                                                                   aGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACg                                                                  >  1:2785932/1‑150 (MQ=255)
                                                                                       cTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCt                                                              >  1:450935/1‑150 (MQ=255)
                                                                                       cTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCt                                                              <  1:2511562/150‑1 (MQ=255)
                                                                                        tGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTg                                                             >  1:3620476/1‑150 (MQ=255)
                                                                                         ggTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGa                                                            >  2:2064626/1‑150 (MQ=255)
                                                                                           tCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGAcc                                                          >  2:1964103/1‑150 (MQ=255)
                                                                                             tACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCtt                                                        >  1:1448798/1‑150 (MQ=255)
                                                                                                ccGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGa                                                     <  1:3003907/150‑1 (MQ=255)
                                                                                                       aCGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTCTCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCaa                                              >  2:4073534/1‑150 (MQ=255)
                                                                                                         gCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAgg                                            <  1:3119370/150‑1 (MQ=255)
                                                                                                         gCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAgg                                            <  1:546825/150‑1 (MQ=255)
                                                                                                          cccTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGc                                           >  2:4053509/1‑150 (MQ=255)
                                                                                                             tGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCTAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTg                                        >  1:506523/1‑150 (MQ=255)
                                                                                                              ggCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTgg                                       <  1:3568997/150‑1 (MQ=255)
                                                                                                                cAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGcc                                     <  1:1443859/150‑1 (MQ=255)
                                                                                                                  aGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGcccc                                   <  2:1124850/150‑1 (MQ=255)
                                                                                                                    tGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCt                                 <  1:4059982/150‑1 (MQ=255)
                                                                                                                    tGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCt                                 >  2:1389511/1‑150 (MQ=255)
                                                                                                                      aTGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTg                               <  1:1407450/150‑1 (MQ=255)
                                                                                                                       tGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGGTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGa                              >  1:3582450/1‑150 (MQ=255)
                                                                                                                            aaCGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGc                         <  2:3762548/150‑1 (MQ=255)
                                                                                                                             aCGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGcc                        >  2:3366389/1‑150 (MQ=255)
                                                                                                                             aCGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGcc                        <  1:4481878/150‑1 (MQ=255)
                                                                                                                              cGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCg                       >  2:2222960/1‑150 (MQ=255)
                                                                                                                               ggCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCgg                      <  1:598716/150‑1 (MQ=255)
                                                                                                                                gCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCggg                     >  2:817813/1‑150 (MQ=255)
                                                                                                                                 ccAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGgc                    >  1:2220089/1‑150 (MQ=255)
                                                                                                                                 ccAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGgc                    <  2:615937/150‑1 (MQ=255)
                                                                                                                                 ccAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGgc                    >  1:3400417/1‑150 (MQ=255)
                                                                                                                                     gTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGa                <  1:3389476/150‑1 (MQ=255)
                                                                                                                                      tGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGaa               <  1:1314774/150‑1 (MQ=255)
                                                                                                                                       gCTCGATGGCGGGCAGGCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAg              >  1:3334484/1‑150 (MQ=39)
                                                                                                                                         tCGATGGCGGGCAGGCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGACGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGca            >  1:2439903/1‑150 (MQ=255)
                                                                                                                                           gATGGCGGGCAGGCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGcaca          >  2:942615/1‑150 (MQ=255)
                                                                                                                                           gATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGcaca          <  1:2717131/150‑1 (MQ=255)
                                                                                                                                            aTGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACAt         >  1:3782522/1‑150 (MQ=255)
                                                                                                                                              ggcgggcaggCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGt       >  2:529369/3‑150 (MQ=255)
                                                                                                                                              ggcgggcaggCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGt       >  2:3394689/3‑150 (MQ=255)
                                                                                                                                                cgggcaggCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCg     >  2:888873/1‑150 (MQ=255)
                                                                                                                                                 gggcgggGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCgg    >  1:556207/1‑150 (MQ=255)
                                                                                                                                                 gggcgggGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCgg    >  1:2903961/1‑150 (MQ=255)
                                                                                                                                                  ggcgggGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCggt   <  1:243442/150‑1 (MQ=255)
                                                                                                                                                  ggcgggGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCggt   <  1:1723830/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTGCATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTG  >  NC_002947/3003637‑3003933

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: ‑