Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,213,155 A→T 10.0% intergenic (+51/‑20) PP_1955 → / → PP_1956 cytochrome P450 family protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,213,1550AT10.0% 165.1 / 8.2 60intergenic (+51/‑20)PP_1955/PP_1956cytochrome P450 family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base T (4/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

TGCCTGAACTCAAGGTGGTCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATT  >  NC_002947/2213013‑2213303
                                                                                                                                              |                                                                                                                                                    
tGCCTGAACTCAAGGTGGTCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCgcgaagcgaa                                                                                                                                                >  1:338329/1‑149 (MQ=255)
  ccTGAACTCAAGGTGGTCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCgcgaagcgaa                                                                                                                                                >  5:172255/1‑147 (MQ=255)
  ccTGAACTCAAGGTGGTCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCgcgaagcgaa                                                                                                                                                <  5:348942/147‑1 (MQ=255)
  ccTGAACTCAAGGTGGTCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCgcgaagcgaa                                                                                                                                                <  6:172255/147‑1 (MQ=255)
  ccTGAACTCAAGGTGGTCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCgcgaagcgaa                                                                                                                                                <  7:101126/147‑1 (MQ=255)
   cTGAACTCAAGGTGGTCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATgg                                                                                                                                             >  1:135231/1‑149 (MQ=255)
      aaCTCAAGGTGGTCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCTTCGCGAATgtcgg                                                                                                                                          <  6:145538/149‑5 (MQ=255)
       aCTCAAGGTGGTCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCTTCGCGAATgtcgg                                                                                                                                          <  3:231119/148‑5 (MQ=255)
        cTCAAGGTGGTCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAgg                                                                                                                                        >  4:188169/1‑149 (MQ=255)
              gtggtCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTa                                                                                                                                  <  3:57206/149‑1 (MQ=255)
               tggtCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTAt                                                                                                                                 <  3:42795/149‑1 (MQ=255)
                    ccTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGaat                                                                                                                             >  7:308485/1‑148 (MQ=255)
                       taTCAGCCGCTCAACTTTGCTCCGAGGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGa                                                                                                                                          <  4:311528/131‑1 (MQ=255)
                        aTCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGaataat                                                                                                                          <  2:238477/147‑1 (MQ=255)
                        aTCACCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGaataat                                                                                                                          <  3:228265/147‑1 (MQ=255)
                           aGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTTCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGt                                                                                                                      <  3:151109/148‑1 (MQ=255)
                           aGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTa                                                                                                                     <  3:144562/149‑1 (MQ=255)
                           aGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTa                                                                                                                     <  8:13157/149‑1 (MQ=255)
                            gCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTac                                                                                                                    >  4:350792/1‑149 (MQ=255)
                              cGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGTATGGTGAGGTTTTTATGAATAATAAGTacg                                                                                                                   >  8:105529/1‑148 (MQ=255)
                                  cAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTacgcacgc                                                                                                              >  3:296454/1‑149 (MQ=255)
                                   atcTTTGCTCCGTTGGCATTGATTCAAGTACGCCTCTCTTTGCAAGTAATCTTGTAAAAACCGTAGCAACCGGCTCCAATTATAACAATATGGTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATTAGTacgcac                                                                                                                <  4:313319/144‑1 (MQ=255)
                                   aaCTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGt                                                                                                                                            >  8:144685/1‑118 (MQ=255)
                                   aaCTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGt                                                                                                                                            <  7:144685/118‑1 (MQ=255)
                                   aaCTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCt                                                                                                             <  3:350792/149‑1 (MQ=255)
                                         gCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTAtctc                                                                                                        <  8:308485/148‑1 (MQ=255)
                                           tCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTcgcg                                                                                                     >  7:18492/1‑149 (MQ=255)
                                                    aTTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGcc                                                                                            >  6:196893/1‑149 (MQ=255)
                                                    aTTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGcc                                                                                            >  8:163184/1‑149 (MQ=255)
                                                             aGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCAtgtg                                                                                   >  4:155349/1‑149 (MQ=255)
                                                             aGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCAtgtg                                                                                   >  5:291163/1‑149 (MQ=255)
                                                              gTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTgg                                                                                  >  8:178410/1‑149 (MQ=255)
                                                                    ctctctTTGCAAGTAATCTGGTAAACACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCTAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTg                                                                            <  7:118458/149‑1 (MQ=255)
                                                                    ctctctTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTg                                                                            <  4:71539/149‑1 (MQ=255)
                                                                    ctctctTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTg                                                                            >  2:320152/1‑149 (MQ=255)
                                                                    ctctctTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTg                                                                            <  5:196893/149‑1 (MQ=255)
                                                                    ctctctTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTg                                                                            <  3:155349/149‑1 (MQ=255)
                                                                           tGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATa                                                                      <  2:77163/148‑1 (MQ=255)
                                                                             cAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACca                                                                   <  4:296454/149‑1 (MQ=255)
                                                                                      tGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCTCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTtgctg                                                          >  4:146914/1‑149 (MQ=255)
                                                                                      tGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGtttt                                                                              <  1:35963/129‑1 (MQ=255)
                                                                                      tGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGtttt                                                                              >  2:35963/1‑129 (MQ=255)
                                                                                      tGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGaaa                                                                        <  3:254304/135‑1 (MQ=255)
                                                                                       ggTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCgg                                                                                                    <  1:191197/106‑1 (MQ=255)
                                                                                       ggTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCgg                                                                                                    >  2:191197/1‑106 (MQ=255)
                                                                                              aCCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCGTGTGGTTTTTGGAAATACCACATTTGCTGCCGaaaa                                                   >  5:314753/1‑148 (MQ=255)
                                                                                              aCCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAat                                                  <  3:188169/149‑1 (MQ=255)
                                                                                                   aGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGTATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCaa                                             <  7:105529/149‑1 (MQ=255)
                                                                                                        ccGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACAttt                                                              >  8:172029/1‑127 (MQ=255)
                                                                                                        ccGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACAttt                                                              <  7:172029/127‑1 (MQ=255)
                                                                                                               cAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCGTGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGc                                   <  6:314753/147‑1 (MQ=255)
                                                                                                                aataatAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGccc                                 >  2:156852/1‑148 (MQ=255)
                                                                                                                 ataataACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCaa                               >  7:70654/1‑149 (MQ=255)
                                                                                                                     taaCAATATGCTGCCGACATTCGCGAAGCGAATGGTGATGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCACTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGACAGACGTCCAAATTgg                           >  2:276036/1‑149 (MQ=255)
                                                                                                                     taaCAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTgg                           >  8:153589/1‑149 (MQ=255)
                                                                                                                       aCAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGATGCAGAAAATATCAAGCCAGCCGGCCCAATTgtag                         >  1:255958/1‑146 (MQ=255)
                                                                                                                            atGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACt                    >  6:229143/1‑149 (MQ=255)
                                                                                                                                  tcaccATTCGCTTCGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAATCCAGCCGGCCCAATTGGCGGCACTGATAtt              >  5:294031/5‑149 (MQ=255)
                                                                                                                                  tcaccATTCGCTTCGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATAtt              >  7:288507/5‑149 (MQ=255)
                                                                                                                                  ccGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATAtt              >  7:283963/1‑149 (MQ=255)
                                                                                                                                    accATTCGCTTCGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTc             >  2:87659/3‑148 (MQ=255)
                                                                                                                                     ccATTCGCTTCGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGa           >  2:53587/2‑149 (MQ=255)
                                                                                                                                      cATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGa           >  3:281055/1‑148 (MQ=255)
                                                                                                                                      cATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGAc          <  2:338329/149‑1 (MQ=255)
                                                                                                                                          cgcgAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAg       <  7:163184/148‑1 (MQ=255)
                                                                                                                                             gaagcgaaTGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGAt   >  5:324483/1‑149 (MQ=255)
                                                                                                                                              aagcgaaTGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGAtt  <  7:153589/149‑1 (MQ=255)
                                                                                                                                              aagcgaaTGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGAt   >  4:300484/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TGCCTGAACTCAAGGTGGTCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATT  >  NC_002947/2213013‑2213303

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: