Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,424,258 A→C 6.4% E321A (GAG→GCG) ‡ PP_5467 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,424,2580AC 6.4% 285.3 / 6.1 94E321A (GAG→GCG) ‡PP_5467hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (44/44);  new base C (2/4);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CTCAAGAGTGCCCTGAGCGCGATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATCGGGCAGAAGCTGGTC  >  NC_002947/1424110‑1424406
                                                                                                                                                     |                                                                                                                                                    
cTCAAGAGTGCCCTGAGCGCGATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAagag                                                                                                                                                     <  1:1970777/150‑1 (MQ=255)
 tCAAGAGTGCCCTGAGCGCGATGCTGATTGGCGCGTGCAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGCAGTTTAGCGGCTTCAAATCGCTGAAAGAgg                                                                                                                                                    >  1:1812697/1‑150 (MQ=255)
 tCAAGAGTGCCCTGAGCGCGATGCTGATTGGCGCGTACATCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAgg                                                                                                                                                    >  1:1861711/1‑150 (MQ=255)
   aaGAGTGCCCTGAGCGCGATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCt                                                                                                                                                  >  1:1693553/1‑150 (MQ=255)
         gCCCTGAGCGCGATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGc                                                                                                                                            <  2:1917981/150‑1 (MQ=255)
           ccTGAGCGCGATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGTTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGCCTCTTTCAGCGa                                                                                                                                          <  1:869102/150‑1 (MQ=255)
            cTGAGCGCGATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCATTGAGGGGCTTTGCCCGAAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGaa                                                                                                                                         <  1:893452/150‑1 (MQ=255)
             tGAGCGCGATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAAt                                                                                                                                        >  1:12099/1‑150 (MQ=255)
             tGAGCGCGATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAAt                                                                                                                                        >  2:259157/1‑150 (MQ=255)
                  gcgATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACcag                                                                                                                                   <  2:1226499/150‑1 (MQ=255)
                  gcgATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACcag                                                                                                                                   <  1:1183815/150‑1 (MQ=255)
                   tgATTCTGATGGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACcagc                                                                                                                                  <  2:868704/149‑1 (MQ=255)
                   cgATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACcagc                                                                                                                                  <  1:91066/150‑1 (MQ=255)
                   cgATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACcagc                                                                                                                                  <  2:1224465/150‑1 (MQ=255)
                    gATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACcagca                                                                                                                                 <  2:2404436/150‑1 (MQ=255)
                         tGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACcag                                                                                                                                   >  2:850626/1‑143 (MQ=255)
                         tGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACcag                                                                                                                                   <  1:850626/143‑1 (MQ=255)
                         tGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACca                                                                                                                                    >  1:1475306/1‑142 (MQ=255)
                         tGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACca                                                                                                                                    <  2:1475306/142‑1 (MQ=255)
                          gATTGGCGCGTACCGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCttt                                                                                                                                                >  2:2058182/1‑129 (MQ=255)
                          gATTGGCGCGTACCGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCttt                                                                                                                                                <  1:2058182/129‑1 (MQ=255)
                           aTTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTc                                                                                                                          <  2:2166143/150‑1 (MQ=255)
                                cgcgTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCaa                                                                                                                     <  1:970133/150‑1 (MQ=255)
                                 gcgTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAg                                                                                                                    <  1:1602418/150‑1 (MQ=255)
                                  cgTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAgg                                                                                                                   >  1:711419/1‑150 (MQ=255)
                                      cAGCCGGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCaaa                                                                                                               <  2:551662/150‑1 (MQ=255)
                                         cccGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGtt                                                                                                            >  1:690608/1‑150 (MQ=255)
                                          ccGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTc                                                                                                           >  2:1567158/1‑150 (MQ=255)
                                                acCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGcc                                                                                                     >  1:1029518/1‑150 (MQ=255)
                                                 ccGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGCCTCTTTCAGCGa                                                                                                                                          >  2:1331062/1‑112 (MQ=255)
                                                 acGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGCCTCTTTCAGCGa                                                                                                                                          <  1:1331062/111‑1 (MQ=255)
                                                  cGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCtt                                                                                                   >  1:945839/1‑150 (MQ=255)
                                                    gCCCGGCGGTGTTGGTTGAAGCGGGGAGCGGCGTTGCCCGCAAGGAAACGGTCGAATGCTCGTGGTGTATCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCt                                                                                                 <  1:2090182/150‑1 (MQ=255)
                                                     ccAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTg                                                                                                >  2:1500697/1‑150 (MQ=255)
                                                        ggCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCa                                                                                             >  1:2329008/1‑150 (MQ=255)
                                                         gCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGCCTCTTTCAGc                                                                                                                                            <  1:908318/102‑1 (MQ=255)
                                                          cGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCagg                                                                                           <  1:259157/150‑1 (MQ=255)
                                                             tgtTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCaggagg                                                                                        >  2:2203353/1‑150 (MQ=255)
                                                               ttGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGtt                                                                                      <  2:690608/150‑1 (MQ=255)
                                                               ttGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGtt                                                                                      <  1:73430/150‑1 (MQ=255)
                                                                  gCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCttt                                                                                                                                                <  2:138221/89‑1 (MQ=255)
                                                                  gCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCttt                                                                                                                                                >  1:138221/1‑89 (MQ=255)
                                                                  gCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGACTTCTGCCAGGAGGTTggg                                                                                   >  2:1869810/1‑150 (MQ=255)
                                                                    tGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCt                                                                                 >  1:73524/1‑150 (MQ=255)
                                                                     gAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCaa                                                                                                                     <  2:2174501/113‑1 (MQ=255)
                                                                     gAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCaa                                                                                                                     >  1:2174501/1‑113 (MQ=255)
                                                                            tGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACcagca                                                                                                                                 >  2:2186365/1‑94 (MQ=255)
                                                                            tGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACcagca                                                                                                                                 <  1:2186365/94‑1 (MQ=255)
                                                                            tGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTcgtcgt                                                                         <  1:1415652/150‑1 (MQ=255)
                                                                             gAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTg                                                                        >  2:266516/1‑150 (MQ=255)
                                                                             gAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTGCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCATTCTGCCAGGAGGTTGGGCTGTCGTCGTg                                                                        <  1:194063/150‑1 (MQ=255)
                                                                              aGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTgg                                                                       >  1:853474/1‑150 (MQ=255)
                                                                              aGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTgg                                                                       <  1:374177/150‑1 (MQ=255)
                                                                               gCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGa                                                                      <  1:1940670/150‑1 (MQ=255)
                                                                                cGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTcgtcgt                                                                         >  1:1796255/1‑146 (MQ=255)
                                                                                cGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTcgtcgt                                                                         <  2:1796255/146‑1 (MQ=255)
                                                                                 ggCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGGCTTCTGCCAGGAGTTTGGACTGTCGTCGTGGAGc                                                                    >  2:633568/1‑150 (MQ=255)
                                                                                  gCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCa                                                                   >  1:1203690/1‑150 (MQ=255)
                                                                                        cccGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTcgtcgt                                                                         <  2:1600115/138‑1 (MQ=255)
                                                                                        cccGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTcgtcgt                                                                         >  1:1600115/1‑138 (MQ=255)
                                                                                         ccGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCaa                                                                                                                     >  2:873641/1‑93 (MQ=255)
                                                                                         ccGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCaa                                                                                                                     <  1:873641/93‑1 (MQ=255)
                                                                                                   aaGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCCCATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCgccaag                                                  <  1:633568/150‑1 (MQ=255)
                                                                                                   aaGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCa                                                                                                        >  1:2136513/1‑96 (MQ=255)
                                                                                                   aaGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCa                                                                                                        <  2:2136513/96‑1 (MQ=255)
                                                                                                      gTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCt                                                                                                 <  1:1767918/100‑1 (MQ=255)
                                                                                                      gTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCt                                                                                                 >  2:1767918/1‑100 (MQ=255)
                                                                                                       tCGAGAACTGCTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCgccaaggcca                                              >  2:590182/1‑150 (MQ=255)
                                                                                                         gagaACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCgccaaggccaag                                            >  1:1146518/1‑150 (MQ=255)
                                                                                                                actCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACcaca                                     >  1:681619/1‑150 (MQ=255)
                                                                                                                  tCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACcagca                                                                                                                                 <  1:2056445/56‑1 (MQ=255)
                                                                                                                  tCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACcagca                                                                                                                                 >  2:2056445/1‑56 (MQ=255)
                                                                                                                    ggTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCCGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGa                                 >  2:592299/1‑150 (MQ=255)
                                                                                                                    ggTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCaaa                                                                                                               >  1:149106/1‑72 (MQ=255)
                                                                                                                    ggTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCaaa                                                                                                               <  2:149106/72‑1 (MQ=255)
                                                                                                                      tGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACt                               >  2:1626862/1‑150 (MQ=255)
                                                                                                                           cTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTg                                                                        >  1:1577854/1‑105 (MQ=255)
                                                                                                                           cTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTg                                                                        <  2:1577854/105‑1 (MQ=255)
                                                                                                                                  ggCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCTAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAGCCGCAGCGACTCGGAATGCGcc                    <  2:1234019/150‑1 (MQ=255)
                                                                                                                                   gCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGc                                                                                                      >  2:1587570/1‑67 (MQ=255)
                                                                                                                                   gCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGc                                                                                                      <  1:1587570/67‑1 (MQ=255)
                                                                                                                                     ttCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGt                                                                 >  1:1909620/1‑102 (MQ=255)
                                                                                                                                     ttCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGt                                                                 <  2:1909620/102‑1 (MQ=255)
                                                                                                                                     ttCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATc                 >  1:1611750/1‑150 (MQ=255)
                                                                                                                                     tggtattCGCTGAAAGCCTCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATc                 >  2:2395773/7‑150 (MQ=255)
                                                                                                                                      tCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATCg                <  1:2304610/150‑1 (MQ=255)
                                                                                                                                          aTCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATCGGGCa            <  1:2203353/150‑1 (MQ=255)
                                                                                                                                           tCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGTGGTTACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATCGGGCAg           <  1:1500697/150‑1 (MQ=255)
                                                                                                                                             gCTGAAAGCCTCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATCg                <  1:1670458/143‑1 (MQ=255)
                                                                                                                                               tGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATCGGGCAGAAGc       <  1:592299/150‑1 (MQ=255)
                                                                                                                                                gAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATCGGGCAGAAGCt      <  2:853474/150‑1 (MQ=255)
                                                                                                                                                 aaaGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCCGACTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATCGGGCAGAAGCTg     >  2:2392495/1‑150 (MQ=255)
                                                                                                                                                 aaaGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGaa                          <  1:698417/129‑1 (MQ=255)
                                                                                                                                                 aaaGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGaa                          >  2:698417/1‑129 (MQ=255)
                                                                                                                                                 aaaGAAGCTTTTAGCAAATAGCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGCAGGTTGGGCTGTCGACGTGGAGCAGTTACTTTGGCGCAAAGGCCTAGAACCACAGCGACTCGGAATACGCCATCAGGCaataactg     >  2:457445/1‑143 (MQ=255)
                                                                                                                                                   agagGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATCGGGCAGAAGCTGGt   <  1:1567158/150‑1 (MQ=255)
                                                                                                                                                   agagGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATCGGGCAGAAGCTGGt   <  2:1693553/150‑1 (MQ=255)
                                                                                                                                                    gagGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATCGGGCAGAAGCTGGTc  <  1:590182/150‑1 (MQ=255)
                                                                                                                                                     agGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATgc                       >  1:1372957/1‑128 (MQ=255)
                                                                                                                                                     agGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATgc                       <  2:1372957/128‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTCAAGAGTGCCCTGAGCGCGATGCTGATTGGCGCGTACAGCCCGTACACCGGCCAGGCGGTGTTGGCTGAAGCGGTGAGCGGCTTTGCCCGCAAGGCCAAGGTCGAGAACTACTCGGTGAGCTACAGCGGCTTCAAATCGCTGAAAGAGGCTTTCAGCGAATACCAGCAGCAAGTCCCCAAGGCAAAGTTCGCAGCCTTCTGCCAGGAGGTTGGGCTGTCGTCGTGGAGCAGTGACTTTGGCGCCAAGGCCAAGAACCACAGCGACTCGGAATGCGCCATCGGGCAGAAGCTGGTC  >  NC_002947/1424110‑1424406

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