Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,883,651 C→G 5.3% A60G (GCC→GGC)  PP_5157 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,883,6510CG 5.3% 303.3 / 7.3 95A60G (GCC→GGC) PP_5157hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (56/33);  new base G (3/2);  total (59/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

CTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCTGC  >  NC_002947/5883506‑5883800
                                                                                                                                                   |                                                                                                                                                     
cTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTg                                                                                                                                                   >  1:1626578/1‑150 (MQ=255)
 tGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACAGAGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTgg                                                                                                                                                  <  1:1676484/150‑1 (MQ=255)
   aTCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGGCCTGGcc                                                                                                                                                >  1:2002556/1‑150 (MQ=255)
     cAGTCCATTGGCTCAGGCAGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGACGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGCGAAGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCaa                                                                                                                                              >  1:1995078/1‑150 (MQ=255)
      aGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGGCCTGGCCaaa                                                                                                                                             <  2:2002556/150‑2 (MQ=255)
        tCCATTGGCTCATGCGGAGCTTAGCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGACTTTGGCCAGGCCCTGGCCAAGGa                                                                                                                                           >  2:97358/1‑150 (MQ=255)
        tCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGa                                                                                                                                           >  2:2117534/1‑150 (MQ=255)
          cATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAAc                                                                                                                                         <  2:705258/150‑1 (MQ=255)
              ggCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAACACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGa                                                                                                                                     <  1:1848739/150‑1 (MQ=255)
                    ggCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTgcg                                                                                                                               >  2:359587/1‑150 (MQ=255)
                     gCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAg                                                                                                                                             >  2:28091/1‑135 (MQ=255)
                     gCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAg                                                                                                                                             <  1:28091/135‑1 (MQ=255)
                         aGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCg                                                                                                                          >  2:1180585/1‑150 (MQ=255)
                          gCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGa                                                                                                                         >  1:1178410/1‑150 (MQ=255)
                                cgacgaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtcgt                                                                                                                   >  2:1952798/1‑150 (MQ=255)
                                 gacgaAATCGCAGACCAGCGCGAGCTGCTCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtcgtc                                                                                                                  >  2:1375387/1‑150 (MQ=255)
                                   cgaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTgc                                                                                                                                <  2:2063992/134‑1 (MQ=255)
                                   cgaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTgc                                                                                                                                >  1:2063992/1‑134 (MQ=255)
                                      aaTCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGc                                                                                                             >  2:2226044/1‑150 (MQ=255)
                                          gCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGCCGACCACCTGGCAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCg                                                                                                         >  1:350836/1‑150 (MQ=255)
                                          gCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCg                                                                                                         >  1:1178051/1‑150 (MQ=255)
                                           ccGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCgg                                                                                                        <  1:1952798/150‑1 (MQ=255)
                                           ccGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCgg                                                                                                        <  1:359587/150‑1 (MQ=255)
                                                   ggCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                                >  2:1464507/1‑150 (MQ=255)
                                                   ggCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTCGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                                >  2:145480/1‑150 (MQ=255)
                                                     cGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGt                                                                                              <  2:1178051/150‑1 (MQ=255)
                                                      gAACTGCGCATTGCCGTACAGGCCGACAGCACGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCATGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTg                                                                                             >  1:1815058/1‑150 (MQ=255)
                                                         cTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTg                                                                                          >  2:1521875/1‑150 (MQ=255)
                                                               aTTGCCGTACAGGCCGCGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGAT‑‑TCGACATCGAATTTAGCCAGGCCCTGGCCAAGGAACTGCATTTGAGCGCCAAATTCGTCGACGCCCCGGCTGCCGAACTGCTAGCgggg                                                                                    >  2:287904/1‑149 (MQ=255)
                                                               aTTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCAGAATTCGTCTAAGCCCCGGCTACCGAAGTGCTGGCGGgt                                                                                    >  2:1937682/1‑150 (MQ=255)
                                                                ttGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGgtg                                                                                   <  1:2226044/150‑1 (MQ=255)
                                                                  gCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACG‑CCAAC‑GACACCTTGCTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCTTGTTCGTCGAAGCCCCGGCTACCGACGTGCTGTCggcgttg                                                                                 <  2:652820/150‑6 (MQ=255)
                                                                      tACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCgaga                                                                             >  1:415334/1‑150 (MQ=255)
                                                                      tACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCCAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCgaga                                                                             <  1:1937682/150‑1 (MQ=255)
                                                                           gCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGa                                                                                                                         >  2:1234399/1‑101 (MQ=255)
                                                                           gCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGa                                                                                                                         <  1:1234399/101‑1 (MQ=255)
                                                                           gCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGc                                                                        >  2:929127/1‑150 (MQ=255)
                                                                              gagaACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAg                                                                     >  2:18270/1‑150 (MQ=255)
                                                                              gagaACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAg                                                                     <  2:170043/150‑1 (MQ=255)
                                                                               agaACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCt                                                                                           >  2:2195369/1‑127 (MQ=255)
                                                                               agaACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCt                                                                                           <  1:2195369/127‑1 (MQ=255)
                                                                                gaACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTa                                                                   <  1:18270/150‑1 (MQ=255)
                                                                                   cTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAg                                                                                                                                             >  2:1690500/1‑73 (MQ=255)
                                                                                   cTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAg                                                                                                                                             <  1:1690500/73‑1 (MQ=255)
                                                                                   cTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGGCGAGTTGGTCGAAGCCCCGGCTACCGAAGTTCTGGCGGGTGTCGAGGGGGGAAAGTACGa                                                                >  2:876512/1‑150 (MQ=255)
                                                                                   cTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGa                                                                >  1:1202139/1‑150 (MQ=255)
                                                                                    tCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCACGGAACTGGATTTGCGCGCCGAGTTCGTCGAATCCACGGCTACCGAAATGATGGTCGGTGTCGAGGTTGACCAGCACGAc                                                               >  2:1082239/1‑150 (MQ=255)
                                                                                    tCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCg                                                                                       >  1:1915282/1‑126 (MQ=255)
                                                                                    tCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCg                                                                                       <  2:1915282/126‑1 (MQ=255)
                                                                                    tCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGAc                                                               >  1:147349/1‑150 (MQ=255)
                                                                                        cATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCTCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGTCGGGTGTCGAGTGTGGCAAGTACGACCTCg                                                           >  2:1654348/1‑150 (MQ=255)
                                                                                           aCGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTAc                                                                  <  1:1304042/140‑1 (MQ=255)
                                                                                           aCGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTAc                                                                  >  2:1304042/1‑140 (MQ=255)
                                                                                            cGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCAtt                                                       <  2:1531412/150‑1 (MQ=255)
                                                                                              ccTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGcccc                                                                                                          <  1:821956/97‑1 (MQ=255)
                                                                                              ccTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGcccc                                                                                                          >  2:821956/1‑97 (MQ=255)
                                                                                               cTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGACCGGGATTTGCGCGCCGAGTTCGTCGAAGCCGCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGAGATCGCATTCAc                                                    >  2:408090/1‑150 (MQ=255)
                                                                                               cTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                                <  1:1259474/106‑1 (MQ=255)
                                                                                               cTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                                >  2:1259474/1‑106 (MQ=255)
                                                                                                ttCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCAcc                                                   >  1:364638/1‑150 (MQ=255)
                                                                                                  caagcaagACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCAcccc                                                 >  1:1774359/1‑150 (MQ=255)
                                                                                                   aagcaagACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCg                                                >  1:265107/1‑150 (MQ=255)
                                                                                                     gcaagACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTc                                              >  1:948028/1‑150 (MQ=255)
                                                                                                       aagACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGa                                            <  1:678311/150‑1 (MQ=255)
                                                                                                        agACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAg                                           >  1:1838659/1‑150 (MQ=255)
                                                                                                         gacgacCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGt                                          <  2:415334/150‑1 (MQ=255)
                                                                                                         gacgacCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGt                                          >  1:887627/1‑150 (MQ=255)
                                                                                                            gacCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGACCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGaa                                       <  1:1654348/150‑1 (MQ=255)
                                                                                                             accaccTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGc                                                                                                             >  2:884903/1‑79 (MQ=255)
                                                                                                             accaccTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGc                                                                                                             <  1:884903/79‑1 (MQ=255)
                                                                                                               caccTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCg                                    >  1:2140514/1‑150 (MQ=255)
                                                                                                                accTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGc                                   <  2:2140514/150‑1 (MQ=255)
                                                                                                                  cTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTc                                 >  1:102120/1‑150 (MQ=255)
                                                                                                                   tGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCGAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTca                                >  2:2257122/1‑150 (MQ=255)
                                                                                                                         gTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTca                                >  1:1906930/1‑144 (MQ=255)
                                                                                                                         gTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTca                                <  2:1906930/144‑1 (MQ=255)
                                                                                                                              tCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACg                       <  1:2117534/150‑1 (MQ=255)
                                                                                                                              tCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACg                       <  2:147349/150‑1 (MQ=255)
                                                                                                                                gACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGTGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGc                     >  2:2226632/1‑150 (MQ=255)
                                                                                                                                gACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGc                     <  2:364638/150‑1 (MQ=255)
                                                                                                                                   aTCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCgg                                                                                      >  2:199514/1‑82 (MQ=255)
                                                                                                                                   aTCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCgg                                                                                      <  1:199514/82‑1 (MQ=255)
                                                                                                                                   aTCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCg                  >  2:137855/1‑150 (MQ=255)
                                                                                                                                   aTCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGCCGGCGGCCCg                  >  2:285258/1‑150 (MQ=255)
                                                                                                                                    tCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTg                                                                                             <  2:1814051/74‑1 (MQ=255)
                                                                                                                                    tCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTg                                                                                             >  1:1814051/1‑74 (MQ=255)
                                                                                                                                    tCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTca                                >  1:1965505/1‑135 (MQ=255)
                                                                                                                                    tCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTca                                <  2:1965505/135‑1 (MQ=255)
                                                                                                                                    tCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGGGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGc                 >  2:168115/1‑150 (MQ=255)
                                                                                                                                    agttccttGGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGc                 >  1:1680858/7‑150 (MQ=255)
                                                                                                                                     cGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCt                <  2:102120/150‑1 (MQ=255)
                                                                                                                                             gCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCAccc                                                  <  2:1516610/108‑1 (MQ=255)
                                                                                                                                             gCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCAccc                                                  >  1:1516610/1‑108 (MQ=255)
                                                                                                                                               cAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGc      >  1:805057/1‑150 (MQ=255)
                                                                                                                                               cAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCAACCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGc      >  1:253400/1‑150 (MQ=255)
                                                                                                                                                   cccTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCATCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCTGc  <  2:1178410/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTT‑‑TCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCTGC  >  NC_002947/5883506‑5883800

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