Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,439,588 G→C 8.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 8.5% 1348.5 / 95.1 388intergenic (‑79/+28)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (197/158);  new base C (22/11);  total (219/169)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 417 total aligned reads displayed.

AGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCG  >  NC_000913/3439439‑3439737
                                                                                                                                                     |                                                                                                                                                     
 ggTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                                                                      >  2:3577965/1‑150 (MQ=255)
  gTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCggg                                                                                                                                                     >  2:1355737/1‑150 (MQ=255)
  gTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCggg                                                                                                                                                     >  1:3836029/1‑150 (MQ=255)
    tATCAAACCCACCTTTCGCTTGCGCTCCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGc                                                                                                                                                   >  2:2960621/1‑150 (MQ=255)
     aTCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCg                                                                                                                                                  >  2:5450415/1‑150 (MQ=255)
      tCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCgg                                                                                                                                                 >  1:6732471/1‑150 (MQ=255)
        aaaCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGt                                                                                                                                               >  2:4870391/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)
           cTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttt                                                                                                                                            >  1:4452654/1‑150 (MQ=255)
            tCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttttt                                                                                                                                           <  1:1536761/150‑1 (MQ=255)
            tCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACCCCAGGGTATGAACAAAAAGAGTACAGCTATTGTATATATGTGGGGATACAACGGGTATAAAAAAACCCGCCGGGGCGGGttttt                                                                                                                                           <  1:386712/150‑1 (MQ=255)
            tCACCTTTCGCTTGCGCTTATGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACAAAAATAATACAGGGGATGAACAAAAAGAGTAAAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttttt                                                                                                                                           >  1:5315162/1‑150 (MQ=255)
             cACCTTGCGCTGGCGCTTCTGCTATATGGCGCTCGGCCCACTGGTCAAGTAACCACATAACTTCTCCAGGGGATGAACAAAAAGAGTACAGCTTTTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                                                                          <  2:5842321/150‑1 (MQ=255)
               ccTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTATATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGTTTTTTTa                                                                                                                                        <  1:6633771/150‑1 (MQ=255)
                 tttCGCTTGCGCTTCTGCTATCTGGCGCTCTCCCCACTGGCCAATTACCCACATAAATACCCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATACGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACg                                                                                                                                      <  1:5656371/150‑1 (MQ=255)
                       gtGCGCTTCTGCTTTATGGCGCTCGGCCCACTGGTCAAGTAACCACATAAATACTCCCGGGTATGAACAAAATGAGAACAGCTTTTGTAGATAAGTGGGGATACTCGGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCtt                                                                                                                                <  1:2857680/149‑1 (MQ=255)
                           gcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTCACCACATAAATACTCCCGCGGATGAACAAACAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTATTTTACGTTGCTTCAGa                                                                                                                            >  1:3621858/1‑150 (MQ=255)
                            cTTCTGCTCTATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGAt                                                                                                                           <  2:5304919/150‑1 (MQ=255)
                                   tataTGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTg                                                                                                                    <  1:381461/150‑1 (MQ=255)
                                        ggCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcagca                                                                                                               <  1:6086082/150‑1 (MQ=255)
                                          cgcgCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAAAAAAAAGAGTACAGCTATGGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCagcagc                                                                                                             <  1:5443636/150‑1 (MQ=255)
                                             tctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTc                                                                                                          >  2:3154102/1‑150 (MQ=255)
                                             tctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTc                                                                                                          >  2:463553/1‑150 (MQ=255)
                                             tctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGt                                                                                                                                               >  2:755448/1‑113 (MQ=255)
                                                   cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                                                                    >  1:6182833/1‑150 (MQ=255)
                                                      tGGTCACGTAACCACAAAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTTTAGATAAGTGGGGATACTACGGTTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtctc                                                                                                 <  1:3018495/150‑1 (MQ=255)
                                                        gTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTg                                                                                               >  1:1758657/1‑150 (MQ=255)
                                                        gTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCGATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTg                                                                                               >  2:1503388/1‑150 (MQ=255)
                                                          cAAGTAACCACCTAAATCCTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                             <  2:3191499/150‑1 (MQ=255)
                                                          cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                             <  2:6752790/150‑1 (MQ=255)
                                                          cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                             <  2:3954640/150‑1 (MQ=255)
                                                             gTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGa                                                                                          >  1:1357016/1‑150 (MQ=255)
                                                              tAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACGCTGCAGCAGCTTCTGCTTTCTCTGAACGAt                                                                                         >  2:5236399/1‑150 (MQ=255)
                                                              tAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACACGCCGGGGCGGGTTTTTTTACTTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGAt                                                                                         >  1:2526609/1‑150 (MQ=255)
                                                               aaCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATc                                                                                        >  1:3551423/1‑150 (MQ=255)
                                                               aaCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATc                                                                                        <  2:4452654/150‑1 (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)
                                                                   acaTATATACTCCAGGGGATGAACAAAAAGAGTACTGCTATTGTAGATAAGTGGGGATACTTCGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                                                                    <  2:615908/150‑1 (MQ=255)
                                                                   acaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                                                                    >  1:18679/1‑150 (MQ=255)
                                                                         aTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCg                                                                              >  1:5188889/1‑150 (MQ=255)
                                                                          tACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGCTCAACCAGCTCGa                                                                             >  1:1320094/1‑150 (MQ=255)
                                                                          tACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGa                                                                             <  1:4279859/150‑1 (MQ=255)
                                                                           aCTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGAt                                                                            <  1:4672902/150‑1 (MQ=255)
                                                                           aCTCCAGGGGATGAACAAAAAGAGTACAGATATTGTAGATAAGTGGGGATACTACGGGAATAAAAAAACCCGCCGGGGCGGGATTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACTATCAACCAGCTCGAt                                                                            <  1:1262087/150‑1 (MQ=255)
                                                                                aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCg                                                                              <  2:5938494/143‑1 (MQ=255)
                                                                                aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                                                                       >  2:5003252/1‑150 (MQ=255)
                                                                                aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                                                                       >  1:2150247/1‑150 (MQ=255)
                                                                                 ggggATGAACAAAAAGAGTAGAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGc                                                                      <  2:6527518/150‑1 (MQ=255)
                                                                                  gggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTCTAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTACGcc                                                                     >  1:5366475/1‑150 (MQ=255)
                                                                                           aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGcgcg                                                            <  1:742474/150‑1 (MQ=255)
                                                                                             aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt                                                          >  1:4728883/1‑150 (MQ=255)
                                                                                              aaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTg                                                         <  2:601137/150‑1 (MQ=255)
                                                                                                  gcaccgCTATTGTAGATACGGGGGGATACCCCGGGTAAAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGc                                                     <  2:2871847/145‑1 (MQ=255)
                                                                                                   tACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGcc                                                    <  2:2210852/150‑1 (MQ=255)
                                                                                                      aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGc                                                 >  1:4573164/1‑150 (MQ=255)
                                                                                                       gCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCa                                                <  1:3148592/150‑1 (MQ=255)
                                                                                                            tGTTGATATGTGGGTATACTAGGGTTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCCTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGaa                                           <  2:26112/150‑1 (MQ=255)
                                                                                                              tAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGc                                         <  1:3212645/150‑1 (MQ=255)
                                                                                                              tAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGc                                         >  1:2958994/1‑150 (MQ=255)
                                                                                                                gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCca                                       <  1:4415040/150‑1 (MQ=255)
                                                                                                                gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCca                                       >  1:4778511/1‑150 (MQ=255)
                                                                                                                gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCca                                       >  2:1035265/1‑150 (MQ=255)
                                                                                                                gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCca                                       <  2:6592892/150‑1 (MQ=255)
                                                                                                                      tGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTc                                 >  1:4037404/1‑150 (MQ=255)
                                                                                                                      tGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTc                                 <  1:81415/150‑1 (MQ=255)
                                                                                                                       ggggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCa                                <  1:1035265/150‑1 (MQ=255)
                                                                                                                       ggggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCa                                >  2:4773739/1‑150 (MQ=255)
                                                                                                                        gggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAg                               >  2:1916387/1‑150 (MQ=255)
                                                                                                                        gggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAg                               >  2:5798551/1‑150 (MQ=255)
                                                                                                                           aTACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAAt                            <  1:4116770/150‑1 (MQ=255)
                                                                                                                              ctacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACg                         <  2:4093485/150‑1 (MQ=255)
                                                                                                                               tacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGa                        >  1:1822185/1‑150 (MQ=255)
                                                                                                                                  gggTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgtg                     >  1:6390919/1‑150 (MQ=255)
                                                                                                                                   aacgtaAAAAAACCCGCCCCGGCGGGTTTTTTTAAGTTTCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCCGCTGGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCCAACTTCAGAATACGAgtgt                    >  1:857350/5‑150 (MQ=255)
                                                                                                                                     tataAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTaa                  >  2:82989/1‑150 (MQ=255)
                                                                                                                                      ataAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAac                 <  2:2405849/150‑1 (MQ=255)
                                                                                                                                       taAAAAAACCCGCCGGGGCGGTTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAacc                <  2:35267/150‑1 (MQ=255)
                                                                                                                                        aaaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCAGCGGCGCGGTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccc               >  2:3425644/1‑149 (MQ=255)
                                                                                                                                         aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAg                               <  1:4801816/133‑1 (MQ=255)
                                                                                                                                         aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTaa                  >  2:329416/1‑146 (MQ=255)
                                                                                                                                         aaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac              >  2:2463228/1‑150 (MQ=255)
                                                                                                                                           aaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTaa                  <  1:48856/144‑1 (MQ=255)
                                                                                                                                           aaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCg            >  2:389198/1‑150 (MQ=255)
                                                                                                                                             aaCCCGCCCCGGCGGGTTTTTTTACGTTTCTTCAGATTACTCTGCAGCAGCGTCGGCTTTCTCTGGACGATCAACCAGCTCGATGTAAGCCCACGGCGCGTTGTCGCCAGCACGGAAGCCACACATCAGAATACGAGTGTAACCACCCGc          >  1:4214461/1‑150 (MQ=255)
                                                                                                                                             aaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACg                         >  2:6620572/1‑135 (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)
                                                                                                                                                 cgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACgg      >  2:6220554/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)
                                                                                                                                                    cgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTc   <  1:5798551/150‑1 (MQ=255)
                                                                                                                                                     ggggCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGa                                                                             >  2:4716404/1‑75 (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: