Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,254,522 T→A 7.2% I38L (ATA→TTA) ‡ 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,5220TA 7.2% 250.4 / 7.3 83I38L (ATA→TTA) ‡PP_4630MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (37/40);  new base A (2/4);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

CGGAGAGCGCATTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTCCCTGCTCGTTCATCTGCACGAAG  >  NC_002947/5254378‑5254667
                                                                                                                                                |                                                                                                                                                   
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                                                                                                                                     <  8:74874/149‑1 (MQ=255)
          aTTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGc                                                                                                                                     <  8:70111/149‑1 (MQ=255)
          aTTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGc                                                                                                                                     >  1:201835/1‑149 (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                                                                                                                             >  6:25416/1‑127 (MQ=255)
                                        gCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAg                                                                                                                             <  5:25416/127‑1 (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                                                                                       <  7:57891/149‑1 (MQ=255)
                                                        cgcTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGa                                                                                       <  1:168196/149‑1 (MQ=255)
                                                        cgcTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGa                                                                                       <  4:377023/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                                                                        <  7:209802/149‑1 (MQ=255)
                                                                       aGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCgg                                                                        <  1:281862/149‑1 (MQ=255)
                                                                       aGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCgg                                                                        <  1:203759/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)
                                                                                                                                                |                                                                                                                                                   
CGGAGAGCGCATTGTCGCGAGCGAGAATGCTGCTCACCTTGCTGACCGCCACGCCGCGCTCCAGCCACCCGAGGATACGGCGGATTTCATCGATATCAGCCAGGGAATACAGGCGGTGACCACTGTCGGTACGGATGGGCTGTATCAGCCCATAACGGCGTTCCCAGGCGCGAAGGGTCACGGCGTTGACACCGGTCAAGCGTGACACTTCACGAATCGGAAACAGCTCCTGCTGCTCGAAACTCGCATCGATCACCGGTTTACGTCCCTGCTCGTTCATCTGCACGAAG  >  NC_002947/5254378‑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: