Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,364,797 G→T 6.9% L158M (CTG→ATG)  PP_2961 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,364,7970GT 6.9% 251.1 / 4.5 81L158M (CTG→ATG) PP_2961LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (33/42);  new base T (2/4);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCACGATGGGTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCTGTAGAAAAGCCCCGATCAGGG  >  NC_002947/3364648‑3364943
                                                                                                                                                     |                                                                                                                                                   
gCCACGATGGGTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAg                                                                                                                                                     <  1:740925/150‑1 (MQ=255)
  cACGATGGGTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGTCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAggg                                                                                                                                                   <  1:523948/150‑1 (MQ=255)
  cACGATGGGTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAATTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAggg                                                                                                                                                   <  1:2062419/150‑1 (MQ=255)
   aCGATGGGTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAgggg                                                                                                                                                  >  1:252084/1‑150 (MQ=255)
   aCGATGGGTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAgggg                                                                                                                                                  >  2:1127337/1‑150 (MQ=255)
    cGATGGGTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCATCATAGACTCGCCGTGTGCTGGTGCCCAGGGGc                                                                                                                                                 <  1:1773837/150‑1 (MQ=255)
     gATGGGTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCa                                                                                                                                                >  2:1823090/1‑150 (MQ=255)
     gATGGGTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCa                                                                                                                                                <  2:2060124/150‑1 (MQ=255)
      aTGGGTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCTAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCAc                                                                                                                                               <  2:851335/150‑1 (MQ=255)
       tGGGTGTAGTGCCAGCGGTCGTACTGCTCAAGCCAGGTAGAGGCGATGCACTTTTGGCCTGCCAGTTCGGCGGGTTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACg                                                                                                                                              <  1:532789/150‑1 (MQ=255)
       tGGGTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACg                                                                                                                                              <  2:644707/150‑1 (MQ=255)
        gggTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCCTGGGCAcca                                                                                                                                             <  2:345516/150‑3 (MQ=255)
         ggTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCCTGGGCAccag                                                                                                                                            <  2:494237/150‑4 (MQ=255)
               gTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAg                                                                                                                                      >  1:569510/1‑150 (MQ=255)
                 gCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGcc                                                                                                                                    >  2:2068073/1‑150 (MQ=255)
                         tCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGctc                                                                                                                            <  2:407404/150‑1 (MQ=255)
                             aaTGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAggg                                                                                                                                                   >  1:1343075/1‑123 (MQ=255)
                             aaTGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAggg                                                                                                                                                   <  2:1343075/123‑1 (MQ=255)
                                 cTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCaccaa                                                                                                                    >  1:1600094/1‑150 (MQ=255)
                                  tCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCac                                                                                                                       <  2:1166079/146‑1 (MQ=255)
                                  tCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCac                                                                                                                       >  1:1166079/1‑146 (MQ=255)
                                    aagcaagGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCaccaaacc                                                                                                                 >  1:1645755/1‑150 (MQ=255)
                                    aagcaagGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCaccaaacc                                                                                                                 >  2:733662/1‑150 (MQ=255)
                                      gcaagGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCaccaaaccaa                                                                                                               >  1:434976/1‑150 (MQ=255)
                                           gTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAgg                                                                                                                                                    >  2:1088210/1‑108 (MQ=255)
                                           gTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAgg                                                                                                                                                    <  1:1088210/108‑1 (MQ=255)
                                                  cGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCa                                                                                                   <  2:1600094/150‑1 (MQ=255)
                                                   gATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGCTTAGCCGTGTGCTGCATTCCGTGGTTACTGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAg                                                                                                  <  1:1487232/150‑1 (MQ=255)
                                                        aGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGTCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGa                                                                                             >  2:2424957/1‑150 (MQ=255)
                                                         gTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGAt                                                                                            <  1:1823090/150‑1 (MQ=255)
                                                             tgGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGttt                                                                                        >  1:825902/1‑150 (MQ=255)
                                                             tgGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGttt                                                                                        <  1:245962/150‑1 (MQ=255)
                                                             tgGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGttt                                                                                        <  2:1946067/150‑1 (MQ=255)
                                                             tgGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGAGCCGCTACACCATGCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGttt                                                                                        <  2:1797493/150‑1 (MQ=255)
                                                               gCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTc                                                                                      >  1:1530090/1‑150 (MQ=255)
                                                                 cTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCCTGGGCAccag                                                                                                                                            >  1:98010/1‑91 (MQ=255)
                                                                 cTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCCTGGGCAccag                                                                                                                                            <  2:98010/94‑4 (MQ=255)
                                                                 cTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATg                                                                                           >  2:1812327/1‑143 (MQ=255)
                                                                 cTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATg                                                                                           <  1:1812327/143‑1 (MQ=255)
                                                                 cTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCa                                                                                    >  2:536029/1‑150 (MQ=255)
                                                                  tGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTc                                                                                      <  1:2176350/147‑1 (MQ=255)
                                                                  tGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTc                                                                                      >  2:2176350/1‑147 (MQ=255)
                                                                       gTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGGTGGCGAGGTCGATGTGTTCGATGAGCa                                                                              >  1:1541991/1‑150 (MQ=255)
                                                                        ttCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGCTGGCCAGGTCGATGTTTTCCATGAGCAg                                                                            <  1:1775475/150‑1 (MQ=255)
                                                                        ttCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCat                                                                             <  1:1127337/150‑2 (MQ=255)
                                                                        ttCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAg                                                                            <  1:2150834/150‑1 (MQ=255)
                                                                              gggCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAgg                                                                                                 >  2:1987664/1‑124 (MQ=255)
                                                                              gggCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAgg                                                                                                 <  1:1987664/124‑1 (MQ=255)
                                                                              gggCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTc                                                                      <  1:582403/150‑1 (MQ=255)
                                                                                 cTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGtt                                                                   <  2:1530090/150‑1 (MQ=255)
                                                                                  tCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTg                                                                  <  1:2174799/150‑1 (MQ=255)
                                                                                         gctggcCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTg                                                                  <  1:1674354/143‑1 (MQ=255)
                                                                                         gctggcCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTg                                                                  >  2:1674354/1‑143 (MQ=255)
                                                                                          ctggcCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGtt                                                          <  2:1109329/150‑1 (MQ=255)
                                                                                           tggcCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTc                                                         <  1:262797/150‑1 (MQ=255)
                                                                                            ggcCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCa                                                        <  2:569510/150‑1 (MQ=255)
                                                                                            ggcCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCa                                                        <  1:940688/150‑1 (MQ=255)
                                                                                               cATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGgc                                                     >  1:621533/1‑150 (MQ=255)
                                                                                                aTGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGccacca                                                                                                                                       <  1:2365834/68‑1 (MQ=255)
                                                                                                aTGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGccacca                                                                                                                                       >  2:2365834/1‑68 (MQ=255)
                                                                                                aTGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGACACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCAAGGGCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGgcg                                                    >  2:314437/1‑150 (MQ=255)
                                                                                                              aGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGca                                      <  1:1172694/150‑1 (MQ=255)
                                                                                                                     tGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCgtg                               >  1:1071655/1‑148 (MQ=255)
                                                                                                                     tGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCggg                               <  1:99683/150‑1 (MQ=255)
                                                                                                                       cagcaTAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGt                             <  2:841685/150‑1 (MQ=255)
                                                                                                                       cagcaTAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGt                             <  2:1645755/150‑1 (MQ=255)
                                                                                                                          caTAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGgc                          >  1:839243/1‑150 (MQ=255)
                                                                                                                            tAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGgcgc                        >  1:1699543/1‑150 (MQ=255)
                                                                                                                             aGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCt                       <  2:252084/150‑1 (MQ=255)
                                                                                                                                 tCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCTGTAg                   <  1:2424957/150‑1 (MQ=255)
                                                                                                                                   gCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCa                                                  >  2:1295371/1‑117 (MQ=255)
                                                                                                                                   gCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCa                                                  <  1:1295371/117‑1 (MQ=255)
                                                                                                                                   gCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCTGTAGaa                 >  2:415302/1‑150 (MQ=255)
                                                                                                                                     cGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGCCGCTGTAGaaaa               >  2:379168/1‑150 (MQ=255)
                                                                                                                                         tgCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGtt                                                          >  1:454688/1‑103 (MQ=255)
                                                                                                                                         tgCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGtt                                                          <  2:454688/103‑1 (MQ=255)
                                                                                                                                         tgCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCTGTAGAAAAGccc           <  2:839243/150‑1 (MQ=255)
                                                                                                                                         tgCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCTGTAGAAAAGccc           <  2:421848/150‑1 (MQ=255)
                                                                                                                                           cTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCTGTAGAAAAGCCCCg         <  1:230200/150‑1 (MQ=255)
                                                                                                                                            tGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTaa                                                               >  2:1611457/1‑95 (MQ=255)
                                                                                                                                            tGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTaa                                                               <  1:1611457/95‑1 (MQ=255)
                                                                                                                                            tGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCt                       <  1:2450412/135‑1 (MQ=255)
                                                                                                                                            tGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCt                       >  2:2450412/1‑135 (MQ=255)
                                                                                                                                            ccgTGCCCCTGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCTGTAGAAAAGCCCCGa        >  2:1672067/3‑150 (MQ=255)
                                                                                                                                              gTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCTGTAGaaa                >  1:1267700/1‑140 (MQ=255)
                                                                                                                                              gTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCTGTAGaaa                <  2:1267700/140‑1 (MQ=255)
                                                                                                                                                gCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCTGTAGAAAAGCCCCGATCAg    >  2:459696/1‑150 (MQ=255)
                                                                                                                                                  ccAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGGTAATAAGTTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGGGTTGTAGAAAAGCCCCGATCAggg  >  1:2368563/1‑150 (MQ=255)
                                                                                                                                                  ccAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGTGGTGGCGCTGTAGAAAAGCCCCGATCAggg  <  2:1071655/150‑1 (MQ=255)
                                                                                                                                                  ccAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCTGTAGAAAAGCCCCGATCAggg  >  1:262075/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCACGATGGGTGAAGTGCCAGCGGTCGTAATGCTCAAGCAAGGTGAAGGCGATGCAGTTGTGGCCTGCCAGTTCGGCGGGCTCGCGTGGCTGGCCATGGCGATCGAGATAGGCCGGTGCAGCATAGACTCGCCGTGTGCTGGTGCCCAGGGGCACGGCCACCAGCCCCTCGCTCTTCACCAAACCAATGCGGATGGCCAGGTCGATGTTTTCCATGAGCAGGTCTTCGTTGTAATCGTTCAGGCGCAGGTCAATCTGCACATCGGGGTGGCGCTGTAGAAAAGCCCCGATCAGGG  >  NC_002947/3364648‑3364943

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