Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,213,154 G→T 10.2% intergenic (+50/‑21) 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,1540GT10.2% 166.7 / 10.1 59intergenic (+50/‑21)PP_1955/PP_1956cytochrome P450 family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base T (4/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-02

TGCCTGAACTCAAGGTGGTCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGAT  >  NC_002947/2213013‑2213302
                                                                                                                                             |                                                                                                                                                    
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                                                                           >  2:320152/1‑149 (MQ=255)
                                                                    ctctctTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTg                                                                           <  3:155349/149‑1 (MQ=255)
                                                                    ctctctTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTg                                                                           <  4:71539/149‑1 (MQ=255)
                                                                    ctctctTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTg                                                                           <  5:196893/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                                                                                                   >  2:191197/1‑106 (MQ=255)
                                                                                       ggTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCgg                                                                                                   <  1:191197/106‑1 (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)
                                                                                                                                             |                                                                                                                                                    
TGCCTGAACTCAAGGTGGTCCCTGATCAGCCGCTCAACTTTGCTCCGATGGCATTGCTTCCAGTACGCCTCTCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGAT  >  NC_002947/2213013‑2213302

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: