Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,020,510 G→T 8.0% intergenic (+21/‑30) ssnA → / → ygfM putative chlorohydrolase/aminohydrolase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,020,5100GT 8.0% 385.5 / 17.7 112intergenic (+21/‑30)ssnA/ygfMputative chlorohydrolase/aminohydrolase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (49/54);  new base T (2/7);  total (51/61)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

Only 100 of 114 total aligned reads displayed.

ACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACAC  >  NC_000913/3020362‑3020651
                                                                                                                                                     |                                                                                                                                             
acaGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑g                                                                                                                                               >  2:376252/1‑149 (MQ=255)
acaGCGTGAGGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATGAATGACAGATGCCCTCTTCCC‑g                                                                                                                                               >  1:923726/1‑149 (MQ=255)
  aGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑Gc                                                                                                                                              <  4:293758/148‑1 (MQ=255)
    cGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAgg                                                                                                                                           <  1:734159/149‑1 (MQ=255)
      tGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGa                                                                                                                                         >  2:44313/1‑149 (MQ=255)
       gATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGaa                                                                                                                                        <  2:1003819/149‑1 (MQ=255)
          ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGGTGGATGCGCTGGCATAAATGACAGATACCCTCTTCCC‑TGCGGGAagag                                                                                                                                     <  2:681718/149‑1 (MQ=255)
           gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAAGAgg                                                                                                                                    <  2:888881/149‑1 (MQ=255)
           gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAAGAgg                                                                                                                                    <  1:846967/149‑1 (MQ=255)
           gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAagag                                                                                                                                     <  3:663675/148‑1 (MQ=255)
           gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAagag                                                                                                                                     <  3:928311/148‑1 (MQ=255)
           gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAgg                                                                                                                                    <  1:45454/149‑1 (MQ=255)
            tCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAggg                                                                                                                                   <  2:618109/149‑1 (MQ=255)
            tCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAgg                                                                                                                                    <  4:588371/148‑1 (MQ=255)
             cAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAAGAggg                                                                                                                                   <  3:143561/148‑1 (MQ=255)
              aaTGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGtt                                                                                                                                 >  1:897764/1‑149 (MQ=255)
              aaTGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGtt                                                                                                                                 >  2:817013/1‑149 (MQ=255)
                 ggTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTaa                                                                                                                               <  3:680721/148‑1 (MQ=255)
                  gtgtGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGC‑GGGAAGAGGGtaagaa                                                                                                                             >  1:15805/1‑144 (MQ=255)
                     tgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAg                                                                                                                          >  2:22216/1‑149 (MQ=255)
                     tgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTAGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAg                                                                                                                          >  1:859546/1‑149 (MQ=255)
                         gctCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGaa                                                                                                                       >  3:687500/3‑148 (MQ=255)
                            cTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAaagaa                                                                                                                   <  2:896256/149‑1 (MQ=255)
                                 aaGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatga                                                                                                              <  1:656570/149‑1 (MQ=255)
                                   gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatgat                                                                                                             <  3:431033/148‑1 (MQ=255)
                                   gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAgatgat                                                                                                             <  4:420187/148‑1 (MQ=255)
                                   gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                            <  1:207240/149‑1 (MQ=255)
                                   gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                            >  1:226476/1‑149 (MQ=255)
                                   gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                            <  1:376252/149‑1 (MQ=255)
                                        tgactgtAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGTGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCA‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAAc                                                                                                        >  4:922561/7‑148 (MQ=255)
                                            tttAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttt                                                                                                   <  1:707253/149‑1 (MQ=255)
                                            tttAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTAGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttt                                                                                                   <  2:859546/149‑1 (MQ=255)
                                             ttAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttt                                                                                                   >  4:573513/1‑148 (MQ=255)
                                              tgacgtCGATGGCGATTCCATTTATGTACAAGCAAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACATATTCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttttt                                                                                                 <  2:1005814/144‑1 (MQ=255)
                                              tAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAtttt                                                                                                  <  3:173855/148‑1 (MQ=255)
                                                aCTTCGATTGCGATTCCATTTATGCACACGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCa                                                                                               >  2:455806/1‑149 (MQ=255)
                                                 cTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAg                                                                                              >  2:140001/1‑149 (MQ=255)
                                                  ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAagaaga                                                                                                                 <  1:521579/129‑1 (MQ=255)
                                                  ttCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAagaaga                                                                                                                 >  2:521579/1‑129 (MQ=255)
                                                   tCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAgg                                                                                             >  3:447348/1‑148 (MQ=255)
                                                   tCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAgg                                                                                             <  4:396620/148‑1 (MQ=255)
                                                    cGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGcc                                                                                           <  1:476097/149‑1 (MQ=255)
                                                        tGCGATTACATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTCTTCAGGCCCGAc                                                                                       <  2:985588/149‑1 (MQ=255)
                                                           gATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGActct                                                                                    >  2:786501/1‑149 (MQ=255)
                                                           gATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGActct                                                                                    >  2:505800/1‑149 (MQ=255)
                                                            aTTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGActct                                                                                    >  3:645822/1‑148 (MQ=255)
                                                             ttacATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGACGTCGGATGGATGCGCTGGCATAAATTACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATCATTCAACACTTTTTCAGCCCCGACTCTg                                                                                   <  4:687500/145‑1 (MQ=255)
                                                               ccATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTc                                                                                 <  3:11073/148‑1 (MQ=255)
                                                                cATTTATGCACCAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATACCCTCTTCCCTGC‑GGGAAGAgggc                                                                                                                                  >  1:561985/1‑97 (MQ=38)
                                                                 aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGACGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGaa                                                                              <  2:923726/149‑1 (MQ=255)
                                                                  tttATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGaa                                                                              >  3:789082/1‑148 (MQ=255)
                                                                   ttATGCACAAGCCAGAAAAGACGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACa                                                                            <  2:1068960/149‑1 (MQ=255)
                                                                    tATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAg                                                                           <  2:472768/149‑1 (MQ=255)
                                                                         acaAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGC‑GGGAAGAgg                                                                                                                                    >  2:1078374/1‑87 (MQ=255)
                                                                         acaAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAAGAgg                                                                                                                                    <  1:1078374/87‑1 (MQ=255)
                                                                          caAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTg                                                                     <  1:130615/149‑1 (MQ=255)
                                                                           aaGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTgg                                                                    >  1:657209/1‑149 (MQ=255)
                                                                           aaGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTgg                                                                    <  2:748968/149‑1 (MQ=255)
                                                                            aGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGa                                                                   >  1:889158/1‑149 (MQ=255)
                                                                             gCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATATATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGaa                                                                  <  2:650279/149‑1 (MQ=255)
                                                                             gCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGaa                                                                  <  2:897764/149‑1 (MQ=255)
                                                                              ccAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGaa                                                                  <  4:789082/148‑1 (MQ=255)
                                                                                  aaaaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGa                                                              <  4:756080/148‑1 (MQ=255)
                                                                                         gCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGTGCTGGAACTGAAGcgccgc                                                      >  1:209102/1‑149 (MQ=255)
                                                                                         gCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGcgccgc                                                      <  1:505800/149‑1 (MQ=255)
                                                                                                 tATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAgg                                               >  4:578242/1‑148 (MQ=255)
                                                                                                     tgGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGa                                           <  3:573513/148‑1 (MQ=255)
                                                                                                             gatggatgCGCTGGCATAAATGACAGATGCCCTCTTCCC‑TGCGGGAag                                                                                                                                       >  4:184552/1‑48 (MQ=255)
                                                                                                              atggatgCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTg                                  >  3:561973/1‑148 (MQ=255)
                                                                                                                  atgCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCg                             >  2:610903/1‑149 (MQ=255)
                                                                                                                      gcTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCgg                         >  2:112431/1‑149 (MQ=255)
                                                                                                                      gcTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCg                          <  4:645822/148‑1 (MQ=255)
                                                                                                                      gcTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATCGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGTTTCGCCgg                         >  2:878511/1‑149 (MQ=255)
                                                                                                                        tGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGACCAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGATCTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCgggg                       <  1:786501/149‑1 (MQ=255)
                                                                                                                        tGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCgggg                       <  1:6816/149‑1 (MQ=255)
                                                                                                                         ggCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCgggg                       >  4:826592/1‑148 (MQ=255)
                                                                                                                         ggCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCgggg                       >  4:727008/1‑148 (MQ=255)
                                                                                                                           cATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgc                     <  4:561973/148‑1 (MQ=255)
                                                                                                                            aTAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACGGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgca                    >  4:226357/1‑148 (MQ=255)
                                                                                                                            aTAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgca                    >  4:236339/1‑148 (MQ=255)
                                                                                                                             tAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcag                   >  4:595729/1‑148 (MQ=255)
                                                                                                                             tAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcag                   >  3:252707/1‑148 (MQ=255)
                                                                                                                                  gACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTgg                                                                    <  1:107497/94‑1 (MQ=255)
                                                                                                                                  gACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTgg                                                                    >  2:107497/1‑94 (MQ=255)
                                                                                                                                    cAGATGCCCTCTTCCC‑GCAGGGTAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGTATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTc            <  3:922561/148‑1 (MQ=255)
                                                                                                                                     aGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTACATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCaa          <  1:455806/149‑1 (MQ=255)
                                                                                                                                      gATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCaa          <  3:595729/148‑1 (MQ=255)
                                                                                                                                       aTGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAAt         >  3:610272/1‑148 (MQ=255)
                                                                                                                                         acccTCTTCCCTGC‑GGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCt      >  1:200544/2‑149 (MQ=255)
                                                                                                                                          cccTCTTCCCTGC‑GGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTa     >  2:773126/1‑149 (MQ=255)
                                                                                                                                          cccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTa     <  2:200544/149‑1 (MQ=255)
                                                                                                                                           ccTCTTCCCTGC‑GGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGaa                                                                  <  3:1051489/87‑1 (MQ=255)
                                                                                                                                           ccTCTTCCCTGC‑GGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGaa                                                                  >  4:1051489/1‑87 (MQ=255)
                                                                                                                                           ccTCTTCCC‑TGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTa     >  4:372822/1‑148 (MQ=255)
                                                                                                                                           ccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCaa          >  1:479323/1‑143 (MQ=255)
                                                                                                                                           ccTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCaa          <  2:479323/143‑1 (MQ=255)
                                                                                                                                            ctctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAGGATGATTGAACAATTTTTCAGGCCCGACTCTTCCGAACAGGCGCTGGACCTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac    >  3:11894/1‑148 (MQ=255)
                                                                                                                                            ctctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac    >  4:34763/1‑148 (MQ=255)
                                                                                                                                            ctctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac    <  3:236339/148‑1 (MQ=255)
                                                                                                                                            ctctTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac    >  4:962581/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
ACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCC‑GCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACAC  >  NC_000913/3020362‑3020651

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑