Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,587,552 A→T 5.1% L290* (TTG→TAG)  treSB ← bifunctional trehalose synthase B/maltokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,587,5520AT 5.1% 1344.2 / 36.2 395L290* (TTG→TAG) treSBbifunctional trehalose synthase B/maltokinase
Reads supporting (aligned to +/- strand):  ref base A (208/166);  new base T (7/14);  total (215/180)
Fisher's exact test for biased strand distribution p-value = 6.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 407 total aligned reads displayed.

GTTCCACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGATCTCG  >  NC_002947/4587404‑4587701
                                                                                                                                                    |                                                                                                                                                     
gtTCCACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCCGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCaa                                                                                                                                                      >  2:4067434/1‑150 (MQ=255)
gtTCCACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCaa                                                                                                                                                      >  1:8380125/1‑150 (MQ=255)
 ttCCACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCATGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAg                                                                                                                                                     <  2:4027596/150‑1 (MQ=255)
   ccACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGc                                                                                                                                                   >  1:3635533/1‑150 (MQ=255)
     acaGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCCGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCa                                                                                                                                                 >  1:1747654/1‑150 (MQ=255)
     acaGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCCGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCa                                                                                                                                                 >  1:790671/1‑150 (MQ=255)
     acaGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAACAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCa                                                                                                                                                 <  2:1688372/150‑1 (MQ=255)
      caGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCAt                                                                                                                                                <  2:4447493/150‑1 (MQ=255)
          tAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCTTGGCCATggaa                                                                                                                                            <  2:256237/150‑4 (MQ=255)
             tCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCGGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATg                                                                                                                                         >  2:6936472/1‑150 (MQ=255)
              cgcgCTCGCGGTCGGTGACCATTTCCAGGGTCGGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCTGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATgc                                                                                                                                        >  2:3737236/1‑150 (MQ=255)
               gcgcTCGCGGTCGGTGACCATTTCCGGGGTCGGTTCATCAGGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATgcg                                                                                                                                       >  1:773605/1‑150 (MQ=255)
                   tCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGc                                                                                                                                   >  2:7420326/1‑150 (MQ=255)
                   tCGCGATCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGc                                                                                                                                   >  2:751781/1‑150 (MQ=255)
                       ggtcggtGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCa                                                                                                                               >  2:3198375/1‑150 (MQ=255)
                            gtGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGa                                                                                                                          <  1:8454695/150‑1 (MQ=255)
                            gtGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGa                                                                                                                          <  1:1027005/150‑1 (MQ=255)
                            gtGACCAGTTCCAGGGTCAGTTCATCATGGTGGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGa                                                                                                                          <  1:5707350/150‑1 (MQ=255)
                             tGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGaa                                                                                                                         <  1:4394844/150‑1 (MQ=255)
                                ccATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGt                                                                                                                      >  2:1344381/1‑150 (MQ=255)
                                ccATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGt                                                                                                                      <  1:6897476/150‑1 (MQ=255)
                                 cATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTg                                                                                                                     >  2:3344251/1‑150 (MQ=255)
                                 cATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTg                                                                                                                     >  1:8310690/1‑150 (MQ=255)
                                 cATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTg                                                                                                                     <  2:1141584/150‑1 (MQ=255)
                                    ttCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGaa                                                                                                                  >  1:2512104/1‑150 (MQ=255)
                                     tCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGaaa                                                                                                                 <  2:3514104/150‑1 (MQ=255)
                                     tCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGaaa                                                                                                                 <  2:3092876/150‑1 (MQ=255)
                                       cAGGGTCAGTTCATTATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGc                                                                                                               >  2:7956150/1‑150 (MQ=255)
                                       cAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGc                                                                                                               >  1:8742725/1‑150 (MQ=255)
                                        aGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGcc                                                                                                              >  2:8263309/1‑150 (MQ=255)
                                         gggTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCa                                                                                                             <  2:4374195/150‑1 (MQ=255)
                                         gggTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCa                                                                                                             <  1:5912942/150‑1 (MQ=255)
                                            tCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCatat                                                                                                          <  1:2196902/150‑1 (MQ=255)
                                                ttCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCa                                                                                                      <  2:1449701/150‑1 (MQ=255)
                                                ttCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCa                                                                                                      >  2:8302519/1‑150 (MQ=255)
                                                ttCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCa                                                                                                      >  1:857982/1‑150 (MQ=255)
                                                   atcatGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATtc                                                                                                   >  1:7419546/1‑150 (MQ=255)
                                                            gcgcAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTcgtcgccgtc                                                                                          <  1:3820302/150‑1 (MQ=255)
                                                              gcAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTcgtcgccgtcgc                                                                                        >  2:2308857/1‑150 (MQ=255)
                                                              gcAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTcgtcgccgtcgc                                                                                        >  2:7106549/1‑150 (MQ=255)
                                                              gcAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTcgtcgccgtcgc                                                                                        >  1:1242266/1‑150 (MQ=255)
                                                                    aGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGc                                                                                  <  2:8733759/150‑1 (MQ=255)
                                                                        cGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGaa                                                                              <  2:7780701/150‑1 (MQ=255)
                                                                        cGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGaa                                                                              <  2:7226069/150‑1 (MQ=255)
                                                                          cccACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCgggaggcg                                                                            >  1:2241335/1‑143 (MQ=255)
                                                                            cACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTAc                                                                          <  1:395824/150‑1 (MQ=255)
                                                                            cACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTAc                                                                          <  1:3020801/150‑1 (MQ=255)
                                                                            cACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACAGGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTAc                                                                          <  1:4842607/150‑1 (MQ=255)
                                                                             aCTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACg                                                                         >  2:1192244/1‑150 (MQ=255)
                                                                              cTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACgg                                                                        <  2:6520183/150‑1 (MQ=255)
                                                                               tGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACggg                                                                       >  2:3629794/1‑150 (MQ=255)
                                                                                 gCAGTTGGCCGGGATTTCCGGGGTGTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGgcg                                                                     >  2:2663313/1‑150 (MQ=255)
                                                                                    gTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCgt                                                                  >  2:2208746/1‑150 (MQ=255)
                                                                                            ggATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGACGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAg                                                          >  2:6356647/1‑150 (MQ=255)
                                                                                             gATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAgg                                                         >  1:4288183/1‑150 (MQ=255)
                                                                                              aTTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGc                                                        >  1:1345429/1‑150 (MQ=255)
                                                                                                ttCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCa                                                      >  2:6125387/1‑150 (MQ=255)
                                                                                                  ccGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACt                                                    >  1:5482939/1‑150 (MQ=255)
                                                                                                  ccGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACt                                                    >  1:4675616/1‑150 (MQ=255)
                                                                                                    ggggTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTgg                                                  >  1:7597926/1‑150 (MQ=255)
                                                                                                      ggTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGtt                                                <  2:4759942/150‑1 (MQ=255)
                                                                                                      ggTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGtt                                                <  2:3902248/150‑1 (MQ=255)
                                                                                                          tGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCt                                            >  2:8116919/1‑150 (MQ=255)
                                                                                                          tGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATACGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCt                                            >  2:8453569/1‑150 (MQ=255)
                                                                                                            gcgcAGGATGTCGGTGAGAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTc                                          <  1:706514/150‑1 (MQ=255)
                                                                                                             cgcAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCg                                         <  1:8269123/150‑1 (MQ=255)
                                                                                                              gcAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCgg                                        >  1:3726371/1‑150 (MQ=255)
                                                                                                              gcAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCgg                                        >  2:5193056/1‑150 (MQ=255)
                                                                                                               cAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGc                                       <  2:5699605/150‑1 (MQ=255)
                                                                                                               cAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGc                                       >  1:4961438/1‑150 (MQ=255)
                                                                                                                    tGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCcagca                                  <  1:3440951/150‑1 (MQ=255)
                                                                                                                         gTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATg                             <  1:4587447/150‑1 (MQ=255)
                                                                                                                              aGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATc                        >  1:937412/1‑150 (MQ=255)
                                                                                                                              aGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATc                        >  1:5255506/1‑150 (MQ=255)
                                                                                                                               gggAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCg                       >  1:6210352/1‑150 (MQ=255)
                                                                                                                                 gAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCggg                     >  2:4652020/1‑150 (MQ=255)
                                                                                                                                 gAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCggg                     <  1:6647702/150‑1 (MQ=255)
                                                                                                                                 gAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCggg                     <  1:414227/150‑1 (MQ=255)
                                                                                                                                   aGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGt                   <  1:6335322/150‑1 (MQ=255)
                                                                                                                                        tCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTg              >  1:3282238/1‑150 (MQ=255)
                                                                                                                                           tacATGGCCTTGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCg           >  1:2649738/3‑150 (MQ=255)
                                                                                                                                           tCCATGGCCAAGGCCATGTACCTGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCg           <  2:534382/150‑1 (MQ=255)
                                                                                                                                            ccATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGt          >  2:2973333/1‑150 (MQ=255)
                                                                                                                                            acATGGCCTTGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGt          >  2:4121661/2‑150 (MQ=255)
                                                                                                                                            acATGGCCTTGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGt          >  1:2804149/2‑150 (MQ=255)
                                                                                                                                            acATGGCCTTGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGt          >  2:1062402/2‑150 (MQ=255)
                                                                                                                                             cATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTc         >  1:7766517/1‑150 (MQ=255)
                                                                                                                                             cATGGCCAAGGCCATGTACATGCGCGGCACCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTc         <  2:1747654/150‑1 (MQ=255)
                                                                                                                                              aTGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCg        >  1:8123113/1‑150 (MQ=255)
                                                                                                                                              aTGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCg        >  1:7195104/1‑150 (MQ=255)
                                                                                                                                              aTGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCg        >  1:2751689/1‑150 (MQ=255)
                                                                                                                                              aTGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCg        >  2:6736287/1‑150 (MQ=255)
                                                                                                                                               tGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGa       >  1:1588988/1‑150 (MQ=255)
                                                                                                                                                  ccAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGAtct    >  1:6557975/1‑150 (MQ=255)
                                                                                                                                                  ccAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGAtct    >  1:208923/1‑150 (MQ=255)
                                                                                                                                                   cAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGAtctc   >  1:7523080/1‑150 (MQ=255)
                                                                                                                                                    aaGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCgtgt                                                                <  2:7379374/88‑1 (MQ=255)
                                                                                                                                                    aaGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGATCTCg  >  2:5994450/1‑150 (MQ=255)
                                                                                                                                                    aaGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGATCTCg  >  2:5557359/1‑150 (MQ=255)
                                                                                                                                                    aaGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGATCTCg  >  2:4974550/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTCCACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTCGATCTCG  >  NC_002947/4587404‑4587701

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