Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,338,494 A→T 9.2% intergenic (+114/+28) PP_3812 → / ← PP_5632 lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,338,4940AT 9.2% 194.2 / 9.2 76intergenic (+114/+28)PP_3812/PP_5632lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base A (26/42);  new base T (5/2);  total (31/45)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01

CGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCG  >  NC_002947/4338349‑4338639
                                                                                                                                                  |                                                                                                                                                   
cGGGCAACGAACTTAAAGTGCATCTGGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAc                                                                                                                                                   >  1:2876/1‑149 (MQ=255)
 gggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAccc                                                                                                                                                 >  6:283740/1‑149 (MQ=255)
   gCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                               <  7:227255/149‑1 (MQ=255)
     aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgcg                                                                                                                                              >  8:229801/1‑148 (MQ=255)
     aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgcg                                                                                                                                              >  1:78896/1‑148 (MQ=255)
       cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                           <  3:273629/149‑1 (MQ=255)
         aaCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAAAAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTAGCGGGTAAACCCGCGAAAgg                                                                                                                                         >  2:281649/1‑149 (MQ=255)
               aaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCgg                                                                                                                                    <  4:228384/148‑1 (MQ=255)
                    cATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATc                                                                                                                              >  1:36277/1‑149 (MQ=255)
                     aTCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTAGCGGGGAAAACAGCGAAAGGGACgggcaata                                                                                                                              >  8:180185/1‑142 (MQ=255)
                     aTCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTAACCGGTAAACCCCCCAAAAGGCCCGCaaaac                                                                                                                              >  8:110566/1‑146 (MQ=255)
                            agtccggCTATCAGCAATAGCCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACACCCCGCCGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCTAAAGGGCCgg                                                                                                                                    <  6:249046/131‑1 (MQ=255)
                                         gCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAccc                                                                                                           <  2:247296/148‑1 (MQ=255)
                                         ccAATAGACCTGCGGGGGGGGTTATTTTGCTTATATTTGCCCAACCCGCCGGCAGGTTTCATTGAACCCGCCCGGGGGGCGGCTGGGCCGGCCCTTCCGGGGGTAAACCCGCTAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCt                                                                                                          <  2:278712/148‑1 (MQ=255)
                                          cAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCtt                                                                                                         <  8:205328/149‑1 (MQ=255)
                                          cAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCtt                                                                                                         <  2:36277/149‑1 (MQ=255)
                                          cAATAGACCCACCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCGCAAACCCtt                                                                                                         >  2:30236/1‑149 (MQ=255)
                                            aTAGACCCGCCTGGGCGGGTCTTTCGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTAGCGGGTAAAACCGCGAAAGGGCCGGAAAATCCAGCTACGGGCTCAAACCCATc                                                                                                        >  8:43670/1‑148 (MQ=255)
                                             tAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCca                                                                                                      <  5:222610/149‑1 (MQ=255)
                                               gACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCGGGTTAACCCGCGAAAGGGCCCGCAAATCCAGCTCAGGGCTCAAACCCTTCcac                                                                                                     >  2:67764/1‑148 (MQ=255)
                                                  ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTAGCCGGGTAACCAGGCCAAGGGCCGGCAAATAATGCTAACGGGGCAAACCCATcaacaaa                                                                                                  >  3:288062/1‑147 (MQ=255)
                                                  ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTAGCGGGGTAACCCGCGGAAGGGCCGGGAAATACGGCTCAGGGGGGAAACCCATGCCCAAtt                                                                                                 >  5:31962/1‑148 (MQ=255)
                                                      cTGGGCGGGTCTTTTGCTTATATTTGCGCAAACTCCCCGCCGGCTTCAATTATTCTCCCAGCGGTGCGGCTGGTCCGGCCCTTTCGCGGTAAAACCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCa                                                                                            <  4:299768/149‑1 (MQ=255)
                                                               tCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAAGATAACCACAACCt                                                                                   >  4:102414/1‑149 (MQ=255)
                                                                  tttGCTTCTTTTTGCGGAAACAGCCCGCAGGATTCATTGAATCTGCCTGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTc                                                                                 <  6:31962/148‑1 (MQ=255)
                                                                   ttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTc                                                                                 <  5:283740/147‑1 (MQ=255)
                                                                       ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAg                                                                                                                                          <  6:196462/87‑1 (MQ=255)
                                                                       ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAg                                                                                                                                          >  5:196462/1‑87 (MQ=255)
                                                                           aTTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCtt                                                                      >  1:22195/1‑149 (MQ=255)
                                                                            tttGCGGAAACCGCCCGCAGGCTTCAATGAATCTGCCAGCGCTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTg                                                                     <  5:265784/149‑1 (MQ=255)
                                                                              tGCGTAAACAGCCCGCAGGCTTCAATGAATCTGCCTGCGGTGGGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTgg                                                                    <  1:281649/148‑1 (MQ=255)
                                                                                  gAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATg                                                                                                 >  5:171117/1‑117 (MQ=255)
                                                                                  gAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATg                                                                                                 <  6:171117/117‑1 (MQ=255)
                                                                                    atcAGCCCGCAGGATTCAATGAATCTGCCTGCGGTGCGGCAGGTCCGGCCCTTTCGGGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGc                                                             <  7:110566/147‑1 (MQ=255)
                                                                                       aGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTg                                                           >  5:79502/1‑148 (MQ=255)
                                                                                       aGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTg                                                           >  7:63783/1‑148 (MQ=255)
                                                                                       aGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGc                                                          <  8:63783/149‑1 (MQ=255)
                                                                                       aGCCCGCAGGCTTAAATTAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGAATCACCTCAGCCTTCGCCACCCCtt                                                                      <  7:192485/138‑1 (MQ=255)
                                                                                          ccGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATg                                                       >  5:160214/1‑149 (MQ=255)
                                                                                            gCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGACCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGcc                                                     <  2:91574/149‑1 (MQ=255)
                                                                                                 tttCAATCAATCTGCCCGGGGTGGGGCTGGGCCGGCCGTTTCGGGGTTAACCCCGCTAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                <  6:53435/148‑1 (MQ=255)
                                                                                                 cgtCAATGAAAGTGCCTGCGGGGCGGCAGGGCCGGCCCTTCCGCGGGTAACCCCGCTAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                <  6:208920/147‑1 (MQ=255)
                                                                                                 cTTCACTGACTCTGCCAGCGGTGCGGCAGGTCCGGCCCTTGCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                <  8:204542/149‑1 (MQ=255)
                                                                                                    cAATGAATCTGCCAGCGGTGCGGCAGTTCCGGCCCTTTGGCGGGTAAACCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTgg                                               <  4:288062/147‑1 (MQ=255)
                                                                                                    cAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTgg                                               <  1:182687/147‑1 (MQ=255)
                                                                                                    cAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGt                                              <  2:67835/148‑1 (MQ=255)
                                                                                                      aTGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATc                                                                                             >  5:256903/1‑100 (MQ=255)
                                                                                                      aTGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATc                                                                                             <  6:256903/100‑1 (MQ=255)
                                                                                                          aTCTGCCTGCGGTGCGGCTGGTCCGGCCCTTTCGCGGTTAACCCGGCTAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAgcg                                        <  6:44095/148‑1 (MQ=255)
                                                                                                          aTCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCa                                                         <  3:198696/131‑1 (MQ=255)
                                                                                                          aTCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCa                                                         >  4:198696/1‑131 (MQ=255)
                                                                                                               ccAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGa                                    <  3:102414/147‑1 (MQ=255)
                                                                                                                    gggtcGGCAGGTCCGGCCCTTTCGCGGGTACACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTaa                             <  3:82505/145‑1 (MQ=255)
                                                                                                                      tGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACa                           <  2:78896/149‑1 (MQ=255)
                                                                                                                         ggCATGTGCGGGCCTTTCCGCGGTAAACCCGCGAAAGGGCCGGGAAATCCTGCTCATGGCTCAAACCCTTCCACAATTATCACCTCATGCTTGCCCACCCCTTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAgg                         <  1:67764/149‑1 (MQ=255)
                                                                                                                         ggCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGc                        <  6:79502/149‑1 (MQ=255)
                                                                                                                           cAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCgg                                <  6:291669/139‑1 (MQ=255)
                                                                                                                           cAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCgg                                >  5:291669/1‑139 (MQ=255)
                                                                                                                            aGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGAc                     >  6:242763/1‑149 (MQ=255)
                                                                                                                                 cGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCt                <  3:28057/149‑1 (MQ=255)
                                                                                                                                 cGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCt                >  4:46376/1‑149 (MQ=255)
                                                                                                                                     ccTTTCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcgc            >  4:160260/1‑149 (MQ=255)
                                                                                                                                     ccTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcgc            <  7:180185/149‑1 (MQ=255)
                                                                                                                                      cTTTCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg           >  5:235379/1‑149 (MQ=255)
                                                                                                                                       tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCctc                                                   <  6:248699/108‑1 (MQ=255)
                                                                                                                                       tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcgc            <  4:51138/147‑1 (MQ=255)
                                                                                                                                       tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcgc            >  3:51138/1‑147 (MQ=255)
                                                                                                                                       tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg           >  2:209400/1‑148 (MQ=255)
                                                                                                                                       tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCGTGGTACAGCGCGGAGCGGTAACAGGCGCCGGCCTGcgcg           >  2:223346/1‑148 (MQ=255)
                                                                                                                                         tCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGt          <  7:270799/147‑1 (MQ=255)
                                                                                                                                         tCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa         <  7:229801/148‑1 (MQ=255)
                                                                                                                                          cgcgGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTAc        <  2:22195/148‑1 (MQ=255)
                                                                                                                                          cccgGGCAACCCCGCGAACGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTAc        <  5:161892/146‑1 (MQ=255)
                                                                                                                                           gcgGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGa      <  3:46376/149‑1 (MQ=255)
                                                                                                                                             tggTAAACCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAAt    <  1:108158/148‑1 (MQ=255)
                                                                                                                                                tAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCg  >  8:43731/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCG  >  NC_002947/4338349‑4338639

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