Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 535,819 A→C 6.6% E85A (GAA→GCA) ‡ rplA → 50S ribosomal protein L1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947535,8190AC 6.6% 239.1 / 6.3 76E85A (GAA→GCA) ‡rplA50S ribosomal protein L1
Reads supporting (aligned to +/- strand):  ref base A (37/34);  new base C (2/3);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTGGGTCCTCG  >  NC_002947/535680‑535966
                                                                                                                                           |                                                                                                                                                   
gtcgAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGTGTCCAGCAGCTTCAgctgctg                                                                                                                                             <  4:274535/148‑2 (MQ=255)
 tcgAGTCCTACGACATCTCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACGAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTTCAgctgctgg                                                                                                                                            <  5:263782/148‑3 (MQ=255)
     gTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGCCCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTg                                                                                                                                       >  2:254451/1‑149 (MQ=255)
      tCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGc                                                                                                                                      <  5:392524/149‑1 (MQ=255)
      tCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGc                                                                                                                                      >  8:45118/1‑149 (MQ=255)
         tACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTTCAgctgc                                                                                                                                               <  3:99776/137‑1 (MQ=255)
         tACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTTCAgctgc                                                                                                                                               >  4:99776/1‑137 (MQ=255)
         tACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCgg                                                                                                                                   >  4:48835/1‑149 (MQ=255)
         tACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCgg                                                                                                                                   >  8:274267/1‑149 (MQ=255)
            gACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGgcgc                                                                                                                                >  7:4153/1‑149 (MQ=255)
               aTCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACAAGGTCGTACGTAGCGCTACCGAGCTGCCCCACGGCACTGGCAAGAACGTTCGCGTTGCTGTCTTCACCCAGGCTCCAGCAGCTGAATCTGCTCTGGCTGCGGCCGCTGa                                                                                                                             >  6:80181/1‑149 (MQ=255)
                tCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGAt                                                                                                                            >  1:381153/1‑149 (MQ=255)
                   ccGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATcgc                                                                                                                         >  7:55321/1‑149 (MQ=255)
                    cGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATcgcg                                                                                                                        >  2:215000/1‑149 (MQ=255)
                      ttAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTa                                                                                                                      <  3:89300/149‑1 (MQ=255)
                      ttAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTa                                                                                                                      <  5:363781/149‑1 (MQ=255)
                          ccTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTa                                                                                                                  <  5:278194/149‑1 (MQ=255)
                           cTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTAt                                                                                                                 >  6:140664/1‑149 (MQ=255)
                              ggTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATg                                                                                                                <  3:48835/147‑1 (MQ=255)
                                tgtTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGAc                                                                                                             >  6:3530/1‑148 (MQ=255)
                                    gACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATc                                                                                                         <  7:45118/148‑1 (MQ=255)
                                    gACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATc                                                                                                         <  3:366524/148‑1 (MQ=255)
                                    gACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCt                                                                                                        >  5:301657/1‑149 (MQ=255)
                                     aCCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTg                                                                                                       <  1:103051/149‑1 (MQ=255)
                                      cccGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGa                                                                                                              >  7:148802/1‑141 (MQ=255)
                                      cccGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGa                                                                                                              <  8:148802/141‑1 (MQ=255)
                                      cccGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGACCTgg                                                                                                      <  5:3530/149‑1 (MQ=255)
                                       ccGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGc                                                                                                     >  6:162356/1‑149 (MQ=255)
                                         gcgTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATcgcg                                                                                                                        >  3:85007/1‑128 (MQ=255)
                                         gcgTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATcgcg                                                                                                                        <  4:85007/128‑1 (MQ=255)
                                           gTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCTAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGcc                                                                                                 <  2:313326/149‑1 (MQ=255)
                                               aTCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCTGTGTCGAGGAGCTGAAGGTGCTCTGGCTGCCGGTGCTGATCGCGTATGTATTGACTATCTGGCTGCCGaaa                                                                                             <  8:14303/149‑1 (MQ=255)
                                                   gACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGaa                                                                                         <  1:200115/149‑1 (MQ=255)
                                                   gACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGa                                                                                          >  4:35087/1‑148 (MQ=255)
                                                      cAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAgg                                                                                      >  1:70886/1‑149 (MQ=255)
                                                         gTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAggcgg                                                                                   <  5:162356/149‑1 (MQ=255)
                                                            gtacgtaGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGgcggcg                                                                                 <  1:254451/148‑1 (MQ=255)
                                                                gtaGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTg                                                                            <  8:55321/149‑1 (MQ=255)
                                                                   gcgcTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGaa                                                                          <  7:217322/148‑1 (MQ=255)
                                                                    cgcTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACt                                                                        <  6:19638/149‑1 (MQ=255)
                                                                     gcTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAATGCGGCGACCTGAACTa                                                                       >  6:287115/1‑149 (MQ=255)
                                                                      cTACCGTGCTGCCACACGGCACTGGCAAGACCTTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTATGTATGGACGATCTGGCTGTCGAAATGAAAGGCGGCGACCTGAACTAt                                                                      <  8:58039/149‑1 (MQ=255)
                                                                        aCCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGaaa                                                                                        <  3:333180/129‑1 (MQ=255)
                                                                        aCCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGaaa                                                                                        >  4:333180/1‑129 (MQ=255)
                                                                         ccGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATg                                                                     <  5:140664/147‑1 (MQ=255)
                                                                            tgctgcCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGAcgtc                                                                <  2:355923/149‑1 (MQ=255)
                                                                               tgcCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGAcgtcgt                                                              >  1:347997/1‑148 (MQ=255)
                                                                                 ccACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTat                                                           >  4:22881/1‑149 (MQ=255)
                                                                                   acacGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTatcg                                                         >  2:42862/1‑149 (MQ=255)
                                                                                   acacGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTatcg                                                         >  4:81115/1‑149 (MQ=255)
                                                                                   acacGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTatcg                                                         <  7:274267/149‑1 (MQ=255)
                                                                                    cacGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTatcgc                                                        <  5:61851/149‑1 (MQ=255)
                                                                                                  gACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCa                                          >  6:43058/1‑149 (MQ=255)
                                                                                                            gTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGt                                <  2:70886/149‑1 (MQ=255)
                                                                                                              tGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAg                              <  6:301657/149‑1 (MQ=255)
                                                                                                               gCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAgg                             >  7:18022/1‑149 (MQ=255)
                                                                                                                cTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGt                            <  1:42862/149‑1 (MQ=255)
                                                                                                                cTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGt                            >  8:224294/1‑149 (MQ=255)
                                                                                                                 tGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGt                            >  8:218473/1‑148 (MQ=255)
                                                                                                                   tCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAg                         >  5:195365/1‑149 (MQ=255)
                                                                                                                     ttCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCt                       >  8:47506/1‑149 (MQ=255)
                                                                                                                      tCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTg                      >  2:36132/1‑149 (MQ=255)
                                                                                                                          ccAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGt                   >  3:93507/1‑148 (MQ=255)
                                                                                                                          ccAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTc                  <  3:81115/149‑1 (MQ=255)
                                                                                                                           cAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGt                   >  4:16931/1‑147 (MQ=255)
                                                                                                                           cAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGt                   >  2:359982/1‑147 (MQ=255)
                                                                                                                            aGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGAGCGTTGTAGGTAAGCTGGGTc                  >  6:269683/1‑147 (MQ=255)
                                                                                                                             gggTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAgg               <  1:36132/149‑1 (MQ=255)
                                                                                                                                tCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAg                              >  2:278931/1‑131 (MQ=255)
                                                                                                                                tCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAg                              <  1:278931/131‑1 (MQ=255)
                                                                                                                                 ccAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCt           <  2:381153/149‑1 (MQ=255)
                                                                                                                                 agagcagcTTCAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCt           >  2:42021/3‑149 (MQ=255)
                                                                                                                                   agcagcTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTgg         <  8:4153/149‑1 (MQ=255)
                                                                                                                                   agcagcTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCCTTGTAGGTCAGCTGGGTCAGGTTCTgg         <  1:215000/149‑1 (MQ=255)
                                                                                                                                       gcTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTGGGTcc     >  5:147690/1‑149 (MQ=255)
                                                                                                                                          gAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTGGGTCCTcg  >  1:40170/1‑149 (MQ=255)
                                                                                                                                           aaGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTGGGTCCTc   <  5:43058/147‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTGGGTCCTCG  >  NC_002947/535680‑535966

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