Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,439,589 G→C 8.4% intergenic (‑80/+27) 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,5890GC 8.4% 1413.7 / 96.4 404intergenic (‑80/+27)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (206/164);  new base C (23/11);  total (229/175)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 419 total aligned reads displayed.

GGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGC  >  NC_000913/3439440‑3439738
                                                                                                                                                     |                                                                                                                                                     
 gTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCggg                                                                                                                                                      >  2:1355737/1‑150 (MQ=255)
   tATCAAACCCACCTTTCGCTTGCGCTCCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGc                                                                                                                                                    >  2:2960621/1‑150 (MQ=255)
     tCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCgg                                                                                                                                                  >  1:6732471/1‑150 (MQ=255)
         aCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttt                                                                                                                                              <  2:301091/150‑1 (MQ=255)
         aCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttt                                                                                                                                              <  2:1814340/150‑1 (MQ=255)
         aCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttt                                                                                                                                              <  2:455662/150‑1 (MQ=255)
           tCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttttt                                                                                                                                            <  1:1536761/150‑1 (MQ=255)
            cACCTTTCGCTTGGGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCTGGGGATTAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGTATACTACGGTTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                                                                           <  1:4293690/150‑1 (MQ=255)
              ccTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTATATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGTTTTTTTa                                                                                                                                         <  1:6633771/150‑1 (MQ=255)
                tttCGCTTGCGCTTCTGCTATCTGGCGCTCTCCCCACTGGCCAATTACCCACATAAATACCCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATACGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACg                                                                                                                                       <  1:5656371/150‑1 (MQ=255)
                tttCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACg                                                                                                                                       >  2:4063531/1‑150 (MQ=255)
                      gtGCGCTTCTGCTTTATGGCGCTCGGCCCACTGGTCAAGTAACCACATAAATACTCCCGGGTATGAACAAAATGAGAACAGCTTTTGTAGATAAGTGGGGATACTCGGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCtt                                                                                                                                 <  1:2857680/149‑1 (MQ=255)
                        gcgcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTct                                                                                                                               >  1:3149244/1‑149 (MQ=255)
                          gcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTCACCACATAAATACTCCCGCGGATGAACAAACAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTATTTTACGTTGCTTCAGa                                                                                                                             >  1:3621858/1‑150 (MQ=255)
                               tGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAct                                                                                                                        <  2:3409818/150‑1 (MQ=255)
                                  tataTGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTg                                                                                                                     <  1:381461/150‑1 (MQ=255)
                                    tatGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATGACTCTgca                                                                                                                   >  2:3496997/1‑150 (MQ=255)
                                       ggCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcagca                                                                                                                <  1:6086082/150‑1 (MQ=255)
                                         cgcTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCagcagc                                                                                                              <  1:5443681/150‑1 (MQ=255)
                                           ctctGCCCCCTGGTCCATGAACCACATAAATCCTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTTGATAAGGGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGTTTTTTTTACGTTGCGTCAGATTACTCTGCAGCAGCtt                                                                                                            <  1:4989270/150‑1 (MQ=255)
                                            tctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTc                                                                                                           >  2:463553/1‑150 (MQ=255)
                                            tctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTc                                                                                                           >  2:3154102/1‑150 (MQ=255)
                                            tctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGt                                                                                                                                                >  2:755448/1‑113 (MQ=255)
                                                cccACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCt                                                                                                       >  1:4805890/1‑150 (MQ=255)
                                                  cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                                                                     >  1:6182833/1‑150 (MQ=255)
                                                     tGGTCACGTAACCACAAAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTTTAGATAAGTGGGGATACTACGGTTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtctc                                                                                                  <  1:3018495/150‑1 (MQ=255)
                                                     tGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATGGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtctc                                                                                                  <  1:5746863/150‑1 (MQ=255)
                                                       gTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTg                                                                                                >  1:1758657/1‑150 (MQ=255)
                                                         cAAGTAACCACCTAAATCCTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                              <  2:3191499/150‑1 (MQ=255)
                                                         cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                              <  2:3954640/150‑1 (MQ=255)
                                                         cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                              <  1:3864410/150‑1 (MQ=255)
                                                         cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                              <  2:6752790/150‑1 (MQ=255)
                                                            gTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGa                                                                                           >  1:1357016/1‑150 (MQ=255)
                                                              aaCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATc                                                                                         <  2:5835106/150‑1 (MQ=255)
                                                              aaCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATc                                                                                         >  1:3551423/1‑150 (MQ=255)
                                                               aCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCa                                                                                        >  2:5249704/1‑150 (MQ=255)
                                                                ccACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCaa                                                                                       >  1:6117778/1‑150 (MQ=255)
                                                                ccACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCaa                                                                                       >  1:6564470/1‑150 (MQ=255)
                                                                 cacaTAAATACTCCAGGGGATGAACAAAAATAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAc                                                                                      <  1:3308624/150‑1 (MQ=255)
                                                                 cacaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAc                                                                                      <  1:6172211/150‑1 (MQ=255)
                                                                  acaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                                                                     >  1:18679/1‑150 (MQ=255)
                                                                    aTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAg                                                                                   >  1:4441794/1‑150 (MQ=255)
                                                                        aTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCg                                                                               >  1:5188889/1‑150 (MQ=255)
                                                                         tACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGCTCAACCAGCTCGa                                                                              >  1:1320094/1‑150 (MQ=255)
                                                                          aCTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGAt                                                                             <  1:4672902/150‑1 (MQ=255)
                                                                          aCTCCAGGGGATGAACAAAAAGAGTACAGATATTGTAGATAAGTGGGGATACTACGGGAATAAAAAAACCCGCCGGGGCGGGATTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACTATCAACCAGCTCGAt                                                                             <  1:1262087/150‑1 (MQ=255)
                                                                            tCCAGGGGATGAACAAAAAGAGTACAGCTATTGGAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTGTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCGTTCTCTGAACGATCAACCAGCTCGAGGt                                                                           >  2:5985657/1‑150 (MQ=255)
                                                                              cAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                         <  2:5936366/150‑1 (MQ=255)
                                                                              cAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                         <  1:4787833/150‑1 (MQ=255)
                                                                               aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                                                                        >  1:2150247/1‑150 (MQ=255)
                                                                                 gggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGcc                                                                      >  1:5366536/1‑150 (MQ=255)
                                                                                 gggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGcc                                                                      >  2:5163226/1‑150 (MQ=255)
                                                                                          aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGcgcg                                                             <  1:742474/150‑1 (MQ=255)
                                                                                          aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGcgcg                                                             >  1:3870737/1‑150 (MQ=255)
                                                                                            aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt                                                           <  2:5432922/150‑1 (MQ=255)
                                                                                                 gcaccgCTATTGTAGATACGGGGGGATACCCCGGGTAAAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGc                                                      <  2:2871847/145‑1 (MQ=255)
                                                                                                  tACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGcc                                                     <  2:2210852/150‑1 (MQ=255)
                                                                                                   aCAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCt                                                    <  2:4998412/150‑1 (MQ=255)
                                                                                                     aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGc                                                  >  2:6192309/1‑150 (MQ=255)
                                                                                                     aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGc                                                  >  1:4495319/1‑150 (MQ=255)
                                                                                                     aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGc                                                  >  1:4573164/1‑150 (MQ=255)
                                                                                                             tAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGc                                          >  1:2958994/1‑150 (MQ=255)
                                                                                                               gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCca                                        <  2:6592892/150‑1 (MQ=255)
                                                                                                               gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCca                                        >  2:1035265/1‑150 (MQ=255)
                                                                                                                     tGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTc                                  >  1:4037404/1‑150 (MQ=255)
                                                                                                                     tGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTc                                  <  1:81415/150‑1 (MQ=255)
                                                                                                                      ggggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCa                                 >  2:4773739/1‑150 (MQ=255)
                                                                                                                      ggggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCa                                 <  1:1035265/150‑1 (MQ=255)
                                                                                                                       gggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAg                                >  2:1916387/1‑150 (MQ=255)
                                                                                                                          aTACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAAt                             <  1:4116770/150‑1 (MQ=255)
                                                                                                                              tacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGa                         <  1:3224362/150‑1 (MQ=255)
                                                                                                                              tacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGa                         >  1:1822185/1‑150 (MQ=255)
                                                                                                                                 gggTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgtg                      >  1:6390919/1‑150 (MQ=255)
                                                                                                                                  aacgtaAAAAAACCCGCCCCGGCGGGTTTTTTTATGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgtgt                     >  2:1521066/5‑150 (MQ=255)
                                                                                                                                  aacgtaAAAAAACCCGCCCCGGCGGGTTTTTTTAAGTTTCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCCGCTGGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCCAACTTCAGAATACGAgtgt                     >  1:857350/5‑150 (MQ=255)
                                                                                                                                    tataAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTaa                   >  2:82989/1‑150 (MQ=255)
                                                                                                                                    cgtaAAAAAACCCGCCCCGGCGGGTTTTTTTACGGTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTaa                   >  2:5632461/3‑150 (MQ=255)
                                                                                                                                     ataAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAac                  <  2:2405849/150‑1 (MQ=255)
                                                                                                                                       aaaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCGTCTGCTTGCTCTGAACGATCAACCCGCGCGCTGTAAGCCATCGGCGCGGTGTCGCCAGCACGGCAGCCCCCCATCAGAATACGAGTGTAacca                >  1:5327070/1‑150 (MQ=255)
                                                                                                                                        aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAg                                >  2:4801816/1‑133 (MQ=255)
                                                                                                                                        aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTaa                   >  2:329416/1‑146 (MQ=255)
                                                                                                                                        aaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac               >  2:6112675/1‑150 (MQ=255)
                                                                                                                                        aaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac               >  2:2463228/1‑150 (MQ=255)
                                                                                                                                         aaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccacc              >  1:253068/1‑150 (MQ=255)
                                                                                                                                          aaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCg             >  2:6517237/1‑150 (MQ=255)
                                                                                                                                          aaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCg             >  2:389198/1‑150 (MQ=255)
                                                                                                                                           aaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAc                                                                                      <  1:4738176/76‑1 (MQ=255)
                                                                                                                                            aaCCCGCCCCGGCGGGTTTTTTTACGTTTCTTCAGATTACTCTGCAGCAGCGTCGGCTTTCTCTGGACGATCAACCAGCTCGATGTAAGCCCACGGCGCGTTGTCGCCAGCACGGAAGCCACACATCAGAATACGAGTGTAACCACCCGc           >  1:4214461/1‑150 (MQ=255)
                                                                                                                                            aaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACg                          >  2:6620572/1‑135 (MQ=255)
                                                                                                                                             aCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCa          <  2:3566954/150‑1 (MQ=255)
                                                                                                                                             aCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCa          <  2:1319176/150‑1 (MQ=255)
                                                                                                                                               ccgccgGGGCGGGTTTTTTTACGTTTCTTAAGATTACTCTGCATCAGCTTCTTCTTTATCTGAACGATCAACCAGCACGATTTAAGCCAACGGCGCGTTGTCGCCAGCACGGAAGCCACACATCAGAATACGAGTGTAACCACCCGCACg        >  1:2431683/1‑150 (MQ=255)
                                                                                                                                               ccgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACg        <  1:5561049/150‑1 (MQ=255)
                                                                                                                                               ccgcccCGGCGGGTTTTTTTATGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACg        >  2:3682888/1‑150 (MQ=255)
                                                                                                                                                cgccgGGGCGGGTTTTTTTACGTTGCTTAAGATTAATCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCAACGGCGCGTTGTCGCCTGCACCGAAGCCACACATCAGAATACGAGTGTAACCACCCGCACgg       >  1:5854967/1‑150 (MQ=255)
                                                                                                                                                cgccgGGGCGGGTTTTTTTAAGTTTGTTCAGATTACTCTTCAGCAGGTTATGCTTTCTGTGTAAGATAAACCAGCTCGATTTAAGCCATAGGCGCGTTTTAGCCTTCACGGAATCCACACATAAGAATAAGAGTGTAAACAACGGCACgg       >  1:2235961/1‑150 (MQ=255)
                                                                                                                                                 gccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGc                                                      <  1:2622695/102‑1 (MQ=255)
                                                                                                                                                 gccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGc                                                      >  2:2622695/1‑102 (MQ=255)
                                                                                                                                                    ggggCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGa                                                                              >  2:4716404/1‑75 (MQ=255)
                                                                                                                                                     gggCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTcgc  >  2:4357786/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGC  >  NC_000913/3439440‑3439738

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: