Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,439,588 G→C 6.5% intergenic (‑79/+28) yhdN ← / ← rplQ DUF1992 family protein/50S ribosomal subunit protein L17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,5880GC 6.5% 1073.2 / 53.2 293intergenic (‑79/+28)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (132/142);  new base C (14/5);  total (146/147)
Fisher's exact test for biased strand distribution p-value = 3.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

Only 100 of 312 total aligned reads displayed.

AGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCG  >  NC_000913/3439439‑3439737
                                                                                                                                                     |                                                                                                                                                     
aGGTTATCAAACTCACCTTTCGCTTGCGCTTCAGCTATATGGCGCACTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTAAGGGTATAAAAAAACacgccg                                                                                                                                                       >  1:1591293/1‑150 (MQ=255)
  gTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCggg                                                                                                                                                     >  2:3622284/1‑150 (MQ=255)
  gTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCggg                                                                                                                                                     >  1:214205/1‑150 (MQ=255)
      tCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCgg                                                                                                                                                 >  2:2600493/1‑150 (MQ=255)
         aaCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtt                                                                                                                                              >  2:3433632/1‑150 (MQ=255)
          aCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttt                                                                                                                                             <  2:959426/150‑1 (MQ=255)
          aCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttt                                                                                                                                             <  1:3948014/150‑1 (MQ=255)
             cACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGtttttt                                                                                                                                          <  1:2868468/150‑1 (MQ=255)
             cACCTTTCGCTTGCGCTTCTGCTATATGGCGCACTGCCCACTGGTAAATTAACCACATAAATACTCCTGGGTATGAACAAAAAGAGTACAGCTTTTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                                                                          <  1:229388/150‑1 (MQ=255)
             cACCTCTCGCTTGCGCTTCGGCGATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAACTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                                                                          <  2:2784232/150‑1 (MQ=255)
              aCCTTTCGCTTGCGCTTCTGCTATATGGCGCGCTGCCCACAGGTCAAGTAACCACAAAAAAACAACCGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACCGGTATAAAAAAACCCGCCGGGGCGGGttttttt                                                                                                                                         >  1:705328/1‑150 (MQ=255)
                cTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTAc                                                                                                                                       >  2:1239672/1‑150 (MQ=255)
                cTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGACGAACAAAAAGATTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTAc                                                                                                                                       <  2:760138/150‑1 (MQ=255)
                  ttCGCTTGCGATTCTGCTATATGGCGCTCTGCCCACTGGACAATTAACCACATAAATACTCCCGGGGATGAACAAAAAGAGTACAGCTATTGTAGAAACGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGt                                                                                                                                     <  1:652534/150‑1 (MQ=255)
                       ttGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCtt                                                                                                                                >  2:861742/1‑150 (MQ=255)
                       ttGCGCTTCTGCTATATGGCGCTCGGCCCACTGTTAAAGTAACCACATAAATACTCCCGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCtt                                                                                                                                <  1:3513074/150‑1 (MQ=255)
                        tGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTc                                                                                                                               >  1:1671944/1‑150 (MQ=255)
                           gcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGa                                                                                                                            <  2:2577911/150‑1 (MQ=255)
                              tCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTa                                                                                                                         <  2:516277/150‑1 (MQ=255)
                                tGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAct                                                                                                                       >  2:3211005/1‑150 (MQ=255)
                                  cTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACCGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTActct                                                                                                                     >  2:3091818/1‑150 (MQ=255)
                                   tataTGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTg                                                                                                                    <  2:3383289/150‑1 (MQ=255)
                                      atGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcag                                                                                                                 <  2:1617881/150‑1 (MQ=255)
                                          cgcTCTTCCAACTGGTCAAGTAACCACCTAAATACTCCAGGGGATGAACAAAAAGATTACAGCTATTTTAGATAAGTGGGGATACGACGTGTATAAAAAAACCCGCCGGGGCTGGTTTTTTTACGTTGCTTCAGATTACTCTGCagcagc                                                                                                             <  1:1960089/150‑1 (MQ=255)
                                          cgcTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCagcagc                                                                                                             >  2:2991209/1‑150 (MQ=255)
                                          cgcTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTATTACGTTGCTTCAGATTACTCTGCagcagc                                                                                                             <  2:2455270/150‑1 (MQ=255)
                                              ctGCCCACTGGTCAAGTAACCACATAACTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCt                                                                                                         <  2:3876102/150‑1 (MQ=255)
                                              ctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGACCAAAAAGAGTACAGCTATTGTAGATAAGAGGGGATACTACGGGTATACAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCttat                                                                                                         <  2:3040397/150‑3 (MQ=255)
                                                   cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCCGCTTCTGCttt                                                                                                    >  2:3709678/1‑150 (MQ=255)
                                                     cTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGAACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtct                                                                                                  <  1:1864443/150‑1 (MQ=255)
                                                        gTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTg                                                                                               <  1:1497219/150‑1 (MQ=255)
                                                        gTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTg                                                                                               >  1:1931298/1‑150 (MQ=255)
                                                         tCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGa                                                                                              >  2:1133302/1‑150 (MQ=255)
                                                         tCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGa                                                                                              >  2:3308482/1‑150 (MQ=255)
                                                          cAAGTAACCACATAAATACTCCCGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                             <  1:2905546/150‑1 (MQ=255)
                                                            aGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACg                                                                                           >  2:3503917/1‑150 (MQ=255)
                                                            aGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACg                                                                                           >  2:393136/1‑150 (MQ=255)
                                                             gTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGa                                                                                          <  1:1733158/150‑1 (MQ=255)
                                                               aaCCACATAAAGACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCCTTCTCTGAACGATc                                                                                        <  2:2677235/150‑1 (MQ=255)
                                                                  cacaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAc                                                                                     >  2:480546/1‑150 (MQ=255)
                                                                    caTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCa                                                                                   >  1:1958178/1‑150 (MQ=255)
                                                                        aaTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTc                                                                               <  1:1516844/150‑1 (MQ=255)
                                                                            cTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATg                                                                           >  1:76669/1‑150 (MQ=255)
                                                                             tCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGt                                                                          >  1:37187/1‑150 (MQ=255)
                                                                             tCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGt                                                                          >  1:2393953/1‑150 (MQ=255)
                                                                                tggtgATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATAATACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAAGATCAACCAGCTCGATGTAAg                                                                       <  1:3634753/146‑1 (MQ=255)
                                                                                aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                                                                       >  1:102087/1‑150 (MQ=255)
                                                                                 ggggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGc                                                                      <  1:3022287/150‑1 (MQ=255)
                                                                                 ggggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGc                                                                      <  1:379626/150‑1 (MQ=255)
                                                                                   ggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCa                                                                    >  1:1038576/1‑150 (MQ=255)
                                                                                   ggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCa                                                                    >  1:118380/1‑150 (MQ=255)
                                                                                   ggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCa                                                                    <  1:3130680/150‑1 (MQ=255)
                                                                                    gATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCAt                                                                   >  1:2442408/1‑150 (MQ=255)
                                                                                     aTGAACAAAAAGAGTACCGCTATTGTAGATCAGTGGGGATACTACGGGTATAAACAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATc                                                                  <  2:2295842/150‑1 (MQ=255)
                                                                                     aTGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATc                                                                  <  1:975407/150‑1 (MQ=255)
                                                                                     aTGAACAAAAAGACTACAGCTATTGTAGATAAGTGGGGATACTACGGCTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATc                                                                  <  2:715012/150‑1 (MQ=255)
                                                                                        aaCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc                                                               <  1:1130923/150‑1 (MQ=255)
                                                                                        aaCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc                                                               >  2:2011850/1‑150 (MQ=255)
                                                                                        aaCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc                                                               <  1:3486476/150‑1 (MQ=255)
                                                                                           aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGcgcg                                                            >  2:583123/1‑150 (MQ=255)
                                                                                            aaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATc                                                                                        <  2:3129904/121‑1 (MQ=255)
                                                                                            aaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGAGCAACCAGCTCGATGTAAGCCATCGGCGCGt                                                           >  2:2889201/1‑150 (MQ=255)
                                                                                             aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt                                                          >  2:1704891/1‑150 (MQ=255)
                                                                                             aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt                                                          >  2:1705095/1‑150 (MQ=255)
                                                                                                gagTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTc                                                       <  1:847381/150‑1 (MQ=255)
                                                                                                    aCAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCt                                                   <  1:937123/150‑1 (MQ=255)
                                                                                                     cAGCTATTGTAGATAAGTGGGGACACTACGGGTATCAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTg                                                  <  1:2001477/150‑1 (MQ=255)
                                                                                                        cTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCAc                                               >  2:3816070/1‑150 (MQ=255)
                                                                                                         tATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACg                                              >  1:3366475/1‑150 (MQ=255)
                                                                                                            tGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGaa                                           <  1:3769441/150‑1 (MQ=255)
                                                                                                                gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCca                                       <  1:2958235/150‑1 (MQ=255)
                                                                                                                  tAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCcaca                                     <  2:1834049/150‑1 (MQ=255)
                                                                                                                   aaGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCacac                                    >  1:1555252/1‑150 (MQ=255)
                                                                                                                       ggggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCa                                <  1:3430902/150‑1 (MQ=255)
                                                                                                                          gATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGaa                             <  1:2986256/150‑1 (MQ=255)
                                                                                                                          gATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTATTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGaa                             <  1:265716/150‑1 (MQ=255)
                                                                                                                              ctacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACg                         <  1:3485272/150‑1 (MQ=255)
                                                                                                                               tacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGa                        <  1:50035/150‑1 (MQ=255)
                                                                                                                                acGGGTATAAATAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGGTTTCTCTGAACGATCAACCAGCTCGATGTAAGCGATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAg                       >  1:2514633/1‑150 (MQ=255)
                                                                                                                                acGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACTTTGCTACATATTACTCTGCAGCAGCTTCTGATTTCTCTTAACGATCAACCAGCTCGATTTAAGCATTCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAg                       <  2:1776065/150‑1 (MQ=255)
                                                                                                                                 cGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgt                      <  1:3251563/150‑1 (MQ=255)
                                                                                                                                   ggTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgtgt                    <  1:3924795/150‑1 (MQ=255)
                                                                                                                                     tataAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTaa                  <  1:3054902/150‑1 (MQ=255)
                                                                                                                                      ataAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAac                 <  2:1029347/150‑1 (MQ=255)
                                                                                                                                      ataAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAac                 <  2:505958/150‑1 (MQ=255)
                                                                                                                                      ataAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGGTGTCGCCTGCACGGAAGCCCCACTTCAGAATACGAGTGTAac                 >  1:1293051/1‑150 (MQ=255)
                                                                                                                                       taAAAAAACCCGCCGGGGCGGGTTTTTTTAAGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAacc                >  1:657210/1‑150 (MQ=255)
                                                                                                                                       taAAAAAACCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCCTCGGCGCGGTGTCGCCTGCACGGACGCCCCACTTCAGAATACGAGTGTAacc                >  2:1899425/1‑150 (MQ=255)
                                                                                                                                        aaaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAacca               >  2:1395340/1‑150 (MQ=255)
                                                                                                                                         aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac              <  1:2723845/150‑1 (MQ=255)
                                                                                                                                         aaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTAGATGTAAGCCATCGGCGCGTTGTCGCCTGCACCGAAGCCACACATCAGAATAAGAGTGTAaccac              >  1:3587581/1‑150 (MQ=255)
                                                                                                                                         aaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTATGCTTTCTCTGAACGATCAACCAGCGCGATGTAAGCCCTCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac              >  1:2073077/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCGTCTGCTTTCTCTGAACGATAAACCAGCTCGATGTTAGCCATAGGCGCGGTGTCGCCAGCACGGAAGCCACACGTCAGAATACGAGTGTAaccacc             >  1:2223790/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCGTCAGATTACTCTGCAGCACCGTCTGCTTTCTCGGGACGATCAACCAGCTCGCTGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccacc             >  1:2536106/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCCGCCCCGGCGGGTTTTTTTAAGTTTCTTCAGATTACTATGCAGCCGCGTCTGCTTTCGCTGAACGATCAACCAGCTCGAGGTAAGCCATCGGCGCGGTGTCGCCTGCACGGAAGCCCCACTTCAGAATACGAGTGTAaacacc             >  1:540239/1‑150 (MQ=255)
                                                                                                                                            aaaCCCGCCGGGGCGGGTTTTTTTACGTTGCGTCAGATTACTCTGCAGCCTCTTATTCTGTATAAGAACGAACCACCCGCAAGATGGAAGCCATCGGCGCGTTGTCGCCAGCACGGAAGCCACACTTCAGAATACGAGTGTAACCCCCgg           >  1:574984/1‑150 (MQ=255)
                                                                                                                                                ccgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACg       >  1:2802944/1‑150 (MQ=255)
                                                                                                                                                ccgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACg       >  2:3420113/1‑150 (MQ=255)
                                                                                                                                                 cgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACgg      <  1:1796759/150‑1 (MQ=255)
                                                                                                                                                     ggggCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTcg  <  1:2608177/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCG  >  NC_000913/3439439‑3439737

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