Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,254,520 T→A 7.1% I38I (ATA→ATT) ‡ PP_4630 ← MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,254,5200TA 7.1% 248.1 / 5.9 83I38I (ATA→ATT) ‡PP_4630MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (38/39);  new base A (2/4);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CCTGCCCGGAGAGCGCATTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTCCCTGCTCGTTCATCTGCACGAAG  >  NC_002947/5254372‑5254667
                                                                                                                                                    |                                                                                                                                                     
ccTGCCCGGAGAGCGCATTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGt                                                                                                                                                       <  2:317311/149‑1 (MQ=255)
 cTGCCCGGAGAGCGCATTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTa                                                                                                                                                      >  2:306250/1‑149 (MQ=255)
      cGGAGAGCGCATTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGc                                                                                                                                                 <  5:367206/149‑1 (MQ=255)
        gagaGCGCATTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGcc                                                                                                                                                >  3:377023/1‑148 (MQ=255)
          gagCGCATTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGATACAGCCCAt                                                                                                                                            <  6:268801/149‑1 (MQ=255)
           agCGCATTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGATACAGCCCAtc                                                                                                                                           <  3:372719/149‑2 (MQ=255)
            gagcATTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCAGATACAGCCCAtcc                                                                                                                                          <  5:61705/147‑3 (MQ=255)
              gcATTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACg                                                                                                                                       <  4:352855/149‑1 (MQ=255)
               cATTGTTGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGATACAGCCCAtccgta                                                                                                                                       <  7:87315/149‑6 (MQ=255)
                aTTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGc                                                                                                                                     >  1:201835/1‑149 (MQ=255)
                aTTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGc                                                                                                                                     <  8:74874/149‑1 (MQ=255)
                aTTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGc                                                                                                                                     <  8:70111/149‑1 (MQ=255)
                        gagcgagAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAg                                                                                                                             >  7:92566/1‑149 (MQ=255)
                               aaTGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGaa                                                                                                                      <  6:375738/149‑1 (MQ=255)
                                 tgctgcTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAgg                                                                                                                    <  3:230886/149‑1 (MQ=255)
                                 atgcgctCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAgg                                                                                                                    <  8:258484/145‑1 (MQ=255)
                                  gctgctCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAggg                                                                                                                   >  2:375678/1‑149 (MQ=255)
                                   ctgctCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGATGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAggg                                                                                                                   >  5:72590/1‑148 (MQ=255)
                                   ctgctCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAggg                                                                                                                   <  7:265111/148‑1 (MQ=255)
                                     gctCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGCCCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCa                                                                                                                >  5:363038/1‑149 (MQ=255)
                                      ctCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCa                                                                                                                <  8:102496/148‑1 (MQ=255)
                                      ctCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCAc                                                                                                               >  4:35898/1‑149 (MQ=255)
                                         aCCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGc                                                                                                            >  2:145767/1‑149 (MQ=255)
                                           cTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGt                                                                                                          <  3:35898/149‑1 (MQ=255)
                                           cTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGt                                                                                                          >  8:209802/1‑149 (MQ=255)
                                             tGCTGACCGCCACGCAGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTg                                                                                                        >  6:205720/1‑149 (MQ=255)
                                              gCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAg                                                                                                                             <  5:25416/127‑1 (MQ=255)
                                              gCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAg                                                                                                                             >  6:25416/1‑127 (MQ=255)
                                              gCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAg                                                                                                                     <  5:359398/135‑1 (MQ=255)
                                              gCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAg                                                                                                                     >  6:359398/1‑135 (MQ=255)
                                              gCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGtt                                                                                                         <  8:155204/147‑1 (MQ=255)
                                               cTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGa                                                                                                       <  3:334161/148‑1 (MQ=255)
                                                tGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGACTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGaca                                                                                                     >  5:69062/1‑149 (MQ=255)
                                                  aCCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACAcc                                                                                                   >  7:163351/1‑149 (MQ=255)
                                                   ccGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCg                                                                                                  <  4:241988/149‑1 (MQ=255)
                                                    cGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTAATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGACCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCgg                                                                                                 >  7:391754/1‑149 (MQ=255)
                                                         cgccgcGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATACACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCaa                                                                                             <  1:90703/149‑1 (MQ=255)
                                                         cgccgcGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAg                                                                                            >  2:203759/1‑149 (MQ=255)
                                                              cgcTCCAGCCACCCGAGGATACGGCTGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTAGGGATGGGCTGTATCAGCCCATAATGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGa                                                                                       <  1:321748/149‑1 (MQ=255)
                                                              cgcTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGa                                                                                       <  4:377023/149‑1 (MQ=255)
                                                              cgcTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGa                                                                                       <  7:57891/149‑1 (MQ=255)
                                                              cgcTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGa                                                                                       <  1:168196/149‑1 (MQ=255)
                                                                 tCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGacac                                                                                    >  4:270704/1‑149 (MQ=255)
                                                                 tCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGacac                                                                                    >  1:414310/1‑149 (MQ=255)
                                                                 tCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGacac                                                                                    >  8:112510/1‑149 (MQ=255)
                                                                  ccAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCaa                                                                                             >  4:419611/1‑139 (MQ=255)
                                                                  ccAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCaa                                                                                             <  3:419611/139‑1 (MQ=255)
                                                                    aGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACtt                                                                                  >  8:339480/1‑148 (MQ=255)
                                                                     gCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCa                                                                                >  5:76126/1‑149 (MQ=255)
                                                                          ccGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAAt                                                                           >  4:399867/1‑149 (MQ=255)
                                                                          ccGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAAt                                                                           <  8:129260/149‑1 (MQ=255)
                                                                            gAGGATCCGGCTGATTTCCTCGATAGCAGCCAGGGATTCGAGGCGGTGACCACTGTCGGTACGTATCGGCTGTATGAGCCCATAACGGCCTTCGCTGGCCCGAAGGGTCACGGCGTTGAGACCTGTCAAGGCTGACTCTTCACGAATCg                                                                         <  4:203715/149‑1 (MQ=255)
                                                                            gAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGACTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCg                                                                         <  6:69062/149‑1 (MQ=255)
                                                                             aGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCgg                                                                        <  1:203759/149‑1 (MQ=255)
                                                                             aGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCgg                                                                        <  7:209802/149‑1 (MQ=255)
                                                                             aGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCgg                                                                        <  1:281862/149‑1 (MQ=255)
                                                                                aTACGGCGGATTTCATCGATATCAGCCAGGGACTACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGaaa                                                                     <  8:198791/149‑1 (MQ=255)
                                                                                aTACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGaaa                                                                     <  7:339480/149‑1 (MQ=255)
                                                                                    ggcggATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGc                                                                 <  8:92566/149‑1 (MQ=255)
                                                                                          tttCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCctgc                                                           <  6:26472/149‑1 (MQ=255)
                                                                                          tttCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCctg                                                            <  6:72590/148‑1 (MQ=255)
                                                                                            tCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCAATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCctgct                                                          >  3:28520/1‑148 (MQ=255)
                                                                                              aTCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTgctgct                                                       >  5:391958/1‑149 (MQ=255)
                                                                                              aTCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTgctgct                                                       <  2:244606/149‑1 (MQ=255)
                                                                                                  atatCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGa                                                    >  6:35007/1‑148 (MQ=255)
                                                                                                   tatCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGaaa                                                  <  4:202393/149‑1 (MQ=255)
                                                                                                       aGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCg                                              >  1:298346/1‑149 (MQ=255)
                                                                                                       aGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAACCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTc                                               <  1:306250/148‑1 (MQ=255)
                                                                                                         ccAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCa                                            <  3:270704/149‑1 (MQ=255)
                                                                                                           aGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATc                                          >  6:299854/1‑149 (MQ=255)
                                                                                                                  aCAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCctgctg                                                         >  1:408355/1‑127 (MQ=255)
                                                                                                                  aCAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCctgctg                                                         <  2:408355/127‑1 (MQ=255)
                                                                                                                    aGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCgg                                 >  3:25121/1‑149 (MQ=255)
                                                                                                                        ggTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTa                             >  4:130798/1‑149 (MQ=255)
                                                                                                                            aCCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTc                         >  2:93061/1‑149 (MQ=255)
                                                                                                                              cACTGTCGGTACGGATGGGCTGTATCAGCCCAGTAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTcc                        >  1:272712/1‑149 (MQ=255)
                                                                                                                                   tCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTCCCTGc                    <  5:252991/147‑1 (MQ=255)
                                                                                                                                      gTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGTGTTGACACCCGTCAAGCGTCACACTTAACGAATCGGAAACAGCTCCTCCTGCTCCAAACTAGCAACCATCACCGGTTTACGTCACTGCTCGt                >  5:244141/1‑148 (MQ=255)
                                                                                                                                        aCGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTCCCTGCTCGTTCa             >  3:93064/1‑149 (MQ=255)
                                                                                                                                           taTGGGCTGATACAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTCCCTGCTCGTTCATCt          >  1:371373/2‑149 (MQ=255)
                                                                                                                                           gATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTCCCTGCTCGTTCATCt          <  5:299854/149‑1 (MQ=255)
                                                                                                                                            aTGGGCTGATACAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTCCCTGCTCGTTCATCTg         >  2:69349/1‑149 (MQ=255)
                                                                                                                                               ggCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTCCCTGCTCGTTCATCTGCAc      >  4:19397/1‑149 (MQ=255)
                                                                                                                                                gCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTCCCTGCTCGTTCATCTGCAc      >  5:217121/1‑148 (MQ=255)
                                                                                                                                                    tATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTCCCTGCTCGTTCATCTGCACGAAg  <  5:205720/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCTGCCCGGAGAGCGCATTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTCCCTGCTCGTTCATCTGCACGAAG  >  NC_002947/5254372‑5254667

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: