Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,982,740 C→G 8.7% R161P (CGC→CCC)  PP_2608 ← shikimate dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,982,7400CG 8.7% 203.2 / 10.8 70R161P (CGC→CCC) PP_2608shikimate dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base G (3/3);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

GCGGACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAA  >  NC_002947/2982596‑2982888
                                                                                                                                                |                                                                                                                                                    
gcggACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCAtgcgtgcg                                                                                                                                                 <  2:1605433/150‑1 (MQ=255)
  ggACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGt                                                                                                                                               >  1:203829/1‑150 (MQ=255)
   gACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACGCATGGGTGCGGta                                                                                                                                              <  2:23843/150‑2 (MQ=255)
    aCAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGc                                                                                                                                             >  2:1772749/1‑150 (MQ=255)
       aCGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCTAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTgg                                                                                                                                          <  1:1369730/150‑1 (MQ=255)
                 tGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTc                                                                                                                                      >  1:302227/1‑144 (MQ=255)
                 tGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTc                                                                                                                                      <  2:302227/144‑1 (MQ=255)
                 tGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCagag                                                                                                                                <  2:1642588/150‑1 (MQ=255)
                  gcgcGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAgg                                                                                                                               <  1:1142641/150‑1 (MQ=255)
                  gcgcGGGTACACATGCCGCCTGGGGAACCGTTGCCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAgg                                                                                                                               <  2:1546333/150‑1 (MQ=255)
                   cgcgGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTggg                                                                                                                                         <  2:1319124/139‑1 (MQ=255)
                   cgcgGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTggg                                                                                                                                         >  1:1319124/1‑139 (MQ=255)
                     cgGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGtt                                                                                                                            <  1:1731690/150‑1 (MQ=255)
                      gggTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTa                                                                                                                           <  2:203829/150‑1 (MQ=255)
                            ccATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTgg                                                                                                                     <  1:140974/150‑1 (MQ=255)
                            ccATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTgg                                                                                                                     <  1:377400/150‑1 (MQ=255)
                            ccATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTgg                                                                                                                     <  1:825620/150‑1 (MQ=255)
                              aTGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGACAAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCg                                                                                                                   >  1:20884/1‑150 (MQ=255)
                                   gACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCACAGAAGGTTAAGCTggacaggtc                                                                                                              >  2:107951/1‑143 (MQ=255)
                                      tGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGc                                                                                                           >  2:1619182/1‑150 (MQ=255)
                                                ttGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGt                                                                                                 <  2:1762712/150‑1 (MQ=255)
                                                ttGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGt                                                                                                 >  2:1262178/1‑150 (MQ=255)
                                                     caccaGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATa                                                                                            >  2:829962/1‑150 (MQ=255)
                                                     caccaGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATa                                                                                            <  1:1293770/150‑1 (MQ=255)
                                                       ccaGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAgg                                                                                          >  1:370484/1‑150 (MQ=255)
                                                         aGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAggaa                                                                                        >  1:216053/1‑148 (MQ=255)
                                                         aGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCa                                                                                        >  2:261575/1‑150 (MQ=255)
                                                             cGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATgg                                                                                    >  2:479221/1‑150 (MQ=255)
                                                               aaGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCa                                                                                  >  2:1280060/1‑150 (MQ=255)
                                                                 gTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACt                                                                                <  1:1801866/150‑1 (MQ=255)
                                                                   cTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGc                                                                              <  1:829962/150‑1 (MQ=255)
                                                                    ttCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGcc                                                                             <  1:1467059/150‑1 (MQ=255)
                                                                    ttCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGcc                                                                             >  1:1889343/1‑150 (MQ=255)
                                                                    ttCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACCGcc                                                                             >  2:330369/1‑150 (MQ=255)
                                                                        aGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCAAGGTCTGCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATTCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAaca                                                                         >  1:1788888/1‑150 (MQ=255)
                                                                             aCTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACAccgcc                                                                    <  2:370484/150‑1 (MQ=255)
                                                                               tGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACcgccgc                                                                  >  2:1028966/1‑150 (MQ=255)
                                                                                gAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCa                                                                 >  1:1161772/1‑150 (MQ=255)
                                                                                     gtgtgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTGTGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGa                                                            <  2:1504276/150‑1 (MQ=255)
                                                                                      tgtgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGAt                                                           >  2:77428/1‑150 (MQ=255)
                                                                                      tgtgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGAt                                                           <  1:1848005/150‑1 (MQ=255)
                                                                                      tgtgCTGACCGTAAGCCCGGGGAAGACGGTGCCGAGCAGGTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCATCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGAt                                                           >  1:448866/1‑150 (MQ=255)
                                                                                        tgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATcc                                                         >  2:340636/1‑149 (MQ=255)
                                                                                           ggagcctAAGCCCGGGCACGCCGATGCCGAGCAGTTCGCATACCGCGCCCATGCGGGCGGTGCTGGGGTAGAAGAGGGTTATGCTGGCGATCCCGGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCa                                                      <  2:1349416/144‑1 (MQ=255)
                                                                                             aCCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTCCCAACACCGCCGCAACCGATCACCCGt                                                    >  2:43992/1‑150 (MQ=255)
                                                                                               cGTAAGCCCGGGGAAGACGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAAt                                                                                      >  2:691703/1‑114 (MQ=255)
                                                                                               cGTAAGCCCGGGGAAGACGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAAt                                                                                      <  1:691703/114‑1 (MQ=255)
                                                                                                    gCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGtt                                             >  2:1718417/1‑150 (MQ=255)
                                                                                                     cccGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGttt                                            <  1:1619182/150‑1 (MQ=255)
                                                                                                                    gCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGaa                             >  1:259734/1‑150 (MQ=255)
                                                                                                                    gCCGAGCAGTTCGCATACCGCACCCATCCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGaa                             <  1:479221/150‑1 (MQ=255)
                                                                                                                      cGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCAc                                                        <  1:26605/121‑1 (MQ=255)
                                                                                                                      cGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCAc                                                        >  2:26605/1‑121 (MQ=255)
                                                                                                                      cGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGc                           <  2:665600/150‑1 (MQ=255)
                                                                                                                       gAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGcc                          >  1:569671/1‑150 (MQ=255)
                                                                                                                        aGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGcc                                                                             >  2:1738028/1‑98 (MQ=255)
                                                                                                                        aGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGcc                                                                             <  1:1738028/98‑1 (MQ=255)
                                                                                                                             ttCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCtt                    >  1:482100/1‑150 (MQ=255)
                                                                                                                             ttCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCtt                    >  1:1615062/1‑150 (MQ=255)
                                                                                                                             ttCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCAATGCCAACACCGCCGAAACCGAGCACCATCACCCGTGTGCCCGCCGGCTCGCAGCCCTGccc                    >  2:1755206/1‑148 (MQ=255)
                                                                                                                             ttCGCATACCGCACCAATGTGTGTGGTGCTGGAGTCTTGTATTGTGATGCTGGCGATTCCTGGCTGACCCAGTGCAGAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCtt                    <  1:1772749/150‑1 (MQ=255)
                                                                                                                              cccagcaCCGCACGCATGGGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTtg                   >  1:588668/7‑150 (MQ=255)
                                                                                                                              cccagcaCCGCACGCATGGGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTtg                   >  2:191191/7‑150 (MQ=255)
                                                                                                                                gCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTtgtg                 <  1:1718417/150‑1 (MQ=255)
                                                                                                                                    aCCGCACGCATGGGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCa             >  2:1675473/1‑150 (MQ=255)
                                                                                                                                    aCCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAAATTATCACTAGGGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCa             <  1:1280060/150‑1 (MQ=255)
                                                                                                                                       gCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAgcg          >  2:1694735/1‑150 (MQ=255)
                                                                                                                                           gcATGGGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGa      <  1:191191/149‑1 (MQ=255)
                                                                                                                                           gcATGGGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGa      <  1:1675473/149‑1 (MQ=255)
                                                                                                                                           ccATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGa      >  1:1829426/1‑150 (MQ=255)
                                                                                                                                             aggcgtgcgGTGCTGGGGTCGCAGAGGGTTATGCTGGCGAACCCTGCCTCAGCAAGTGCATAGGCAATAGCAATGCCAACACCGCCGCAACCGACCACCAGTTCCCGTTGGCCCGCCTGCGCGAAGGCATGCTGGTGGGCAGCGCCGAgg    >  1:1716206/3‑150 (MQ=255)
                                                                                                                                             aTGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCCGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAgg    <  2:922526/150‑1 (MQ=255)
                                                                                                                                               gcgtgcgGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGaa  <  1:1755206/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCGGACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAA  >  NC_002947/2982596‑2982888

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