Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,857,758 A→C 6.4% E84A (GAG→GCG) ‡ PP_5134 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,857,7580AC 6.4% 197.3 / 2.6 70E84A (GAG→GCG) ‡PP_5134hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/28);  new base C (3/2);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTCCCACCGCCCCAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATGCGCGCGAGCCGGTGCTGGCC  >  NC_002947/5857610‑5857906
                                                                                                                                                     |                                                                                                                                                    
gtgtCCCACCGCCCCAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGACCCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAg                                                                                                                                                     >  1:403651/1‑150 (MQ=255)
gtgtCCCACCGCCCCAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAg                                                                                                                                                     <  1:779270/150‑1 (MQ=255)
   tCCCACCGCCCCAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAgcgc                                                                                                                                                  <  2:499366/150‑1 (MQ=255)
     ccACCGCCCCAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCt                                                                                                                                                >  2:1145773/1‑150 (MQ=255)
      cACCGCCCCAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGCTCGCCTg                                                                                                                                               >  1:939545/1‑150 (MQ=255)
      cACCGCCCCAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTg                                                                                                                                               >  1:1381691/1‑150 (MQ=255)
             ccAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGCTCGCCTgggcgctc                                                                                                                                        <  2:939545/150‑7 (MQ=255)
             ccAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGcc                                                                                                                                        <  1:633358/150‑1 (MQ=255)
             ccAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGcc                                                                                                                                        <  2:1318424/150‑1 (MQ=255)
                    tactCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGGCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCg                                                                                                                                 >  1:372766/1‑150 (MQ=255)
                    tactCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCg                                                                                                                                 >  2:521194/1‑150 (MQ=255)
                      ctCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGc                                                                                                                               <  2:1381691/150‑1 (MQ=255)
                       tCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGACCAGGCGAGCGCCTGGACGGCCCTGCCCGGCt                                                                                                                              <  1:160694/150‑1 (MQ=255)
                           cATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAg                                                                                                                          >  2:8981/1‑150 (MQ=255)
                            aTGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGc                                                                                                                         >  1:1311382/1‑150 (MQ=255)
                                 gAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGtttt                                                                                                                    <  1:521194/150‑1 (MQ=255)
                                 gAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGtttt                                                                                                                    >  2:752629/1‑150 (MQ=255)
                                       gCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAg                                                                                                              <  1:1204258/150‑1 (MQ=255)
                                         aCGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAggg                                                                                                            >  2:862500/1‑150 (MQ=255)
                                           ggCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGccc                                                                                                                                  <  2:598361/126‑1 (MQ=255)
                                           ggCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGccc                                                                                                                                  >  1:598361/1‑126 (MQ=255)
                                             cTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGCCCCGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGc                                                                                                        <  2:403651/150‑1 (MQ=255)
                                                tCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGaa                                                                                                     >  2:64114/1‑150 (MQ=255)
                                                         ggTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGCAgcgcgcgcc                                                                                            >  2:931801/1‑147 (MQ=255)
                                                               ccTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGTACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTc                                                                                      <  2:921570/150‑1 (MQ=255)
                                                               ccTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTc                                                                                      <  2:1311382/150‑1 (MQ=255)
                                                                 tttGCTGGAGCTGTACGACACGCCTGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGa                                                                                    >  1:1413721/1‑150 (MQ=255)
                                                                 tttGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACTGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGa                                                                                    >  2:1046557/1‑150 (MQ=255)
                                                                        gAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGtt                                                                                                                      <  2:1064420/109‑1 (MQ=255)
                                                                        gAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGtt                                                                                                                      >  1:1064420/1‑109 (MQ=255)
                                                                                 gACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGcc                                                                                                                                        <  1:712340/82‑1 (MQ=255)
                                                                                 gACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGcc                                                                                                                                        >  2:712340/1‑82 (MQ=255)
                                                                                 gACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTa                                                                    >  2:608036/1‑150 (MQ=255)
                                                                                   cacGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTAtt                                                                  >  1:156700/1‑150 (MQ=255)
                                                                                    acGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTg                                                                 <  1:931801/150‑1 (MQ=255)
                                                                                     cGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTgc                                                                <  1:608036/150‑1 (MQ=255)
                                                                                      gCCAGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTgcg                                                               >  1:415956/1‑150 (MQ=255)
                                                                                              tGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGtttt                                                                                                                    <  1:865766/89‑1 (MQ=255)
                                                                                              tGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGtttt                                                                                                                    >  2:865766/1‑89 (MQ=255)
                                                                                              tGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAActgc                                                       >  1:1281819/1‑150 (MQ=255)
                                                                                                  aGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGc                                                   <  2:947381/150‑1 (MQ=255)
                                                                                                      gcgaTCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGt                                               >  1:1424358/4‑150 (MQ=255)
                                                                                                       ccATCGCCCTGCTCGACTACCTTGAACGCCTCGAGCGCCGAAGCGAGCCCCTGGACCTGCCCGCCCCGCTGGAGCGGTTTCTGTACGGCGGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACCTCTGGCCAGTc                                              >  1:219305/1‑150 (MQ=255)
                                                                                                        cATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCa                                                 >  1:347429/1‑146 (MQ=255)
                                                                                                        cATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCa                                                 <  2:347429/146‑1 (MQ=255)
                                                                                                         aTCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAgg                                                                                                             >  1:624485/1‑85 (MQ=255)
                                                                                                         aTCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAgg                                                                                                             <  2:624485/85‑1 (MQ=255)
                                                                                                          tCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATg                                           <  1:65124/150‑1 (MQ=255)
                                                                                                                   tCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGt                                  <  2:1281819/150‑1 (MQ=255)
                                                                                                                      aCTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTAtg                               >  2:1316773/1‑150 (MQ=255)
                                                                                                                        tACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTAtgtg                             >  2:887207/1‑150 (MQ=255)
                                                                                                                            ttGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCg                         >  1:545674/1‑150 (MQ=255)
                                                                                                                              gAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGAt                       <  1:8981/150‑1 (MQ=255)
                                                                                                                               ctcGCCTGGGCGCTCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCa                                                                         <  1:919598/97‑1 (MQ=255)
                                                                                                                               ctcGCCTGGGCGCTCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCa                                                                         >  2:919598/3‑99 (MQ=255)
                                                                                                                                aCGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATgc                     >  2:396189/1‑150 (MQ=255)
                                                                                                                                  ggagggAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCAAGTCATGCCGGCCTGTATGTGATCGATgcgc                   >  1:410267/5‑150 (MQ=255)
                                                                                                                                   ccTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATgcgcg                  <  1:862500/150‑1 (MQ=255)
                                                                                                                                    cTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCaggag                                                                             >  2:1417205/1‑90 (MQ=255)
                                                                                                                                    cTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCaggag                                                                             <  1:1417205/90‑1 (MQ=255)
                                                                                                                                    cTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATgcgcgc                 >  2:740834/1‑150 (MQ=255)
                                                                                                                                    cTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATgcgcgc                 >  2:1058451/1‑150 (MQ=255)
                                                                                                                                         gccgtccAGGCGCTCGCCTGGGCGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATGCGCGCGAGc             >  2:933574/6‑150 (MQ=255)
                                                                                                                                           cgtccAGGCGCTCGCCTGGGCGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGa                            >  1:331574/4‑133 (MQ=255)
                                                                                                                                           cgtccAGGCGCTCGCCTGGGCGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGa                            <  2:331574/130‑1 (MQ=255)
                                                                                                                                                aGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATGCGCGCGAGCCGGTGCt      >  2:329929/1‑150 (MQ=255)
                                                                                                                                                 ggCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATGCGCGCGAGCCGGTGCTg     <  2:156700/150‑1 (MQ=255)
                                                                                                                                                 ggCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTACTGGCCAGTCATGCCGGCCTGTATGTGATCGATGCGCGCGAGCCGGTGCTg     <  1:752629/150‑1 (MQ=255)
                                                                                                                                                    gAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATGCGCGCGAGCCGGTGCTGGcc  >  2:957240/1‑150 (MQ=255)
                                                                                                                                                    gAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATGCGCGCGAGCCGGTGCTGGcc  >  1:1173202/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTGTCCCACCGCCCCAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATGCGCGCGAGCCGGTGCTGGCC  >  NC_002947/5857610‑5857906

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