Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,439,589 G→C 8.9% 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.9% 721.1 / 50.1 214intergenic (‑80/+27)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (100/95);  new base C (8/11);  total (108/106)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-01

Only 100 of 216 total aligned reads displayed.

CTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT‑‑CGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTA  >  NC_000913/3439319‑3439860
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                 
                                                                                                                   cTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCAGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAg                                                                                                                                                            <  2:59593/275‑1 (MQ=255)
                                                                                                                          gTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                                                                                                                                                                                                   >  7:379228/1‑149 (MQ=255)
                                                                                                                             aTCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGc                                                                                                                                                                                                                                                                                <  8:276732/149‑1 (MQ=255)
                                                                                                                               cAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCgg                                                                                                                                                                                                                                                                              <  5:1098989/149‑1 (MQ=255)
                                                                                                                                 aaCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATCGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGt                                                                                                                                                                                                                                                                            >  8:778663/1‑149 (MQ=255)
                                                                                                                                  aCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGtt                                                                                                                                                                                                                                                                           <  8:902350/149‑1 (MQ=255)
                                                                                                                                   cTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttt                                                                                                                                                                                                                                                                          >  8:719607/1‑149 (MQ=255)
                                                                                                                                   cTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttt                                                                                                                                                                                                                                                                          <  6:460829/149‑1 (MQ=255)
                                                                                                                                   cTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttt                                                                                                                                                                                                                                                                          <  5:317864/149‑1 (MQ=255)
                                                                                                                                    tCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttt                                                                                                                                                                                                                                                                         <  5:1078444/149‑1 (MQ=255)
                                                                                                                                    tCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttt                                                                                                                                                                                                                                                                         <  5:652472/149‑1 (MQ=255)
                                                                                                                                    tCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGtttt                                                                                                                                                                                                                                                                         <  7:762041/149‑1 (MQ=255)
                                                                                                                                    tCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGtttt                                                                                                                                                                                                                                                                         <  5:838599/149‑1 (MQ=255)
                                                                                                                                    tCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGATCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttt                                                                                                                                                                                                                                                                         <  5:1125816/149‑1 (MQ=255)
                                                                                                                                     cACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttttt                                                                                                                                                                                                                                                                        <  7:318181/149‑1 (MQ=255)
                                                                                                                                     cACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttttt                                                                                                                                                                                                                                                                        <  5:401683/149‑1 (MQ=255)
                                                                                                                                     cACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttttt                                                                                                                                                                                                                                                                        >  8:1398247/1‑149 (MQ=255)
                                                                                                                                      aCCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                                                                                                                                                                                                       >  5:535431/1‑149 (MQ=255)
                                                                                                                                          ttCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGCAGCCACACTTCAGAATACGAGTGTAACCACCGGCa                                                                                                                                      >  2:223964/1‑274 (MQ=255)
                                                                                                                                               ttGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCt                                                                                                                                                                                                                                                              >  6:774782/1‑149 (MQ=255)
                                                                                                                                                    cTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGa                                                                                                                                                                                                                                                         <  7:981198/149‑1 (MQ=255)
                                                                                                                                                      tCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCCGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT‑‑CGCGa                                                                                                                          >  1:38684/1‑272 (MQ=255)
                                                                                                                                                         gCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAct                                                                                                                                                                                                                                                    >  7:727714/1‑149 (MQ=255)
                                                                                                                                                          cTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTActc                                                                                                                                                                                                                                                   <  7:177303/149‑1 (MQ=255)
                                                                                                                                                          cTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTActc                                                                                                                                                                                                                                                   <  7:1327187/149‑1 (MQ=255)
                                                                                                                                                               tGGCGCTGTGCCCACTGGTCAAGTAACCACAGAAATACTCCAGGGAATAAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT‑‑CGCGAAACGCGGGCCCa                                                                                                              <  2:23607/275‑1 (MQ=255)
                                                                                                                                                                ggCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcagc                                                                                                                                                                                                                                             >  6:1092747/1‑149 (MQ=255)
                                                                                                                                                                    ctctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTGGGTAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT‑‑CGCGAAACGCGGGCCCAGTTc                                                                                                          >  2:93185/1‑274 (MQ=255)
                                                                                                                                                                           cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCtt                                                                                                                                                                                                                                  >  5:212338/1‑149 (MQ=255)
                                                                                                                                                                           cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCtt                                                                                                                                                                                                                                  >  7:1187123/1‑149 (MQ=255)
                                                                                                                                                                           cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCtt                                                                                                                                                                                                                                  >  8:1235083/1‑149 (MQ=255)
                                                                                                                                                                           cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCtt                                                                                                                                                                                                                                  >  5:99768/1‑149 (MQ=255)
                                                                                                                                                                              tGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtct                                                                                                                                                                                                                               <  6:27656/149‑1 (MQ=255)
                                                                                                                                                                                 tCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTg                                                                                                                                                                                                                            >  6:167820/1‑149 (MQ=255)
                                                                                                                                                                                    aGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                                                                                                                                                                                                         <  6:590057/149‑1 (MQ=255)
                                                                                                                                                                                      tAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGa                                                                                                                                                                                                                       <  6:516620/149‑1 (MQ=255)
                                                                                                                                                                                               aaaTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGc                                                                                                                                                                                                              >  7:288433/1‑149 (MQ=255)
                                                                                                                                                                                                 aTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGTTTTTTTatac                                                                                                                                                                                                                                                                  >  8:1262892/1‑92 (MQ=38)
                                                                                                                                                                                                 aTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGTTTTTTTatac                                                                                                                                                                                                                                                                  <  7:1262892/95‑4 (MQ=255)
                                                                                                                                                                                                     tCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGa                                                                                                                                                                                                                                                         <  5:711181/100‑1 (MQ=255)
                                                                                                                                                                                                     tCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATg                                                                                                                                                                                                        >  7:1432292/1‑149 (MQ=255)
                                                                                                                                                                                                       cAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                                                                                                                                                                                                    >  2:128004/1‑151 (MQ=255)
                                                                                                                                                                                                        aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                                                                                                                                                     <  6:212338/149‑1 (MQ=255)
                                                                                                                                                                                                        aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                                                                                                                                                     <  5:574820/149‑1 (MQ=255)
                                                                                                                                                                                                        aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                                                                                                                                                     <  5:1285295/149‑1 (MQ=255)
                                                                                                                                                                                                           ggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGcc                                                                                                                                                                                                  >  8:1160582/1‑149 (MQ=255)
                                                                                                                                                                                                              tGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT‑‑c                                                                                                                              >  1:11011/1‑212 (MQ=255)
                                                                                                                                                                                                                 aCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAc                                                                                                                                                                                                                  >  7:1392809/1‑127 (MQ=255)
                                                                                                                                                                                                                 aCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc                                                                                                                                                                                            >  8:1391810/1‑149 (MQ=255)
                                                                                                                                                                                                                 aCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc                                                                                                                                                                                            >  8:438729/1‑149 (MQ=255)
                                                                                                                                                                                                                   aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcgc                                                                                                                                                                                          >  7:1021446/1‑149 (MQ=255)
                                                                                                                                                                                                                    aaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCt                                                                                                                                                                                                                                   <  7:865825/107‑1 (MQ=255)
                                                                                                                                                                                                                    aaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCt                                                                                                                                                                                                                                   >  8:865825/1‑107 (MQ=255)
                                                                                                                                                                                                                    aaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGcgcg                                                                                                                                                                                         >  6:235938/1‑149 (MQ=255)
                                                                                                                                                                                                                      aaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt                                                                                                                                                                                       <  8:379228/149‑1 (MQ=255)
                                                                                                                                                                                                                            aCAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGc                                                                                                                                                                                                                                    >  6:351430/1‑98 (MQ=255)
                                                                                                                                                                                                                             cAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCt                                                                                                                                                                                <  6:451117/149‑1 (MQ=255)
                                                                                                                                                                                                                              aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTg                                                                                                                                                                               >  8:804117/1‑149 (MQ=255)
                                                                                                                                                                                                                              aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTg                                                                                                                                                                               >  5:1282899/1‑149 (MQ=255)
                                                                                                                                                                                                                              aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTg                                                                                                                                                                               <  5:1455598/149‑1 (MQ=255)
                                                                                                                                                                                                                                 tATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCAc                                                                                                                                                                            <  5:167820/149‑1 (MQ=255)
                                                                                                                                                                                                                                    tttAGGTAAGTTGGGGATATCCCGGGTTAAAAAAACCCGGCCGGGTCGGGTTTTTTTCCGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT‑‑CGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGc                                           <  4:83611/272‑1 (MQ=255)
                                                                                                                                                                                                                                        gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACTATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGcc                                                                                                                                                                     >  7:1364724/1‑149 (MQ=255)
                                                                                                                                                                                                                                         aTAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCca                                                                                                                                                                    <  6:448804/149‑1 (MQ=255)
                                                                                                                                                                                                                                          tAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT‑‑CGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACg                                                                         >  4:158724/1‑237 (MQ=255)
                                                                                                                                                                                                                                             gTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACt                                                                                                                                                                >  7:596301/1‑149 (MQ=255)
                                                                                                                                                                                                                                               tGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT‑‑CGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTa                                                                     <  2:158724/237‑1 (MQ=255)
                                                                                                                                                                                                                                                ggggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTc                                                                                                                                                              >  8:1325902/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 gggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCa                                                                                                                                                             <  8:63946/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 gggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT‑‑CGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTTCCACGAACTCGTTATCACGGATACGGGCGAATTCCAGACGACCATTAGAAACGCtctacgtgt                             >  4:97517/1‑267 (MQ=255)
                                                                                                                                                                                                                                                   gATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTTAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGa                                                                                                                                                           >  5:387140/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   gATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT‑‑CGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCCCGAGTACGGGCGAATGCCAGACGACGATTAGCACAGCTATCAGTCtt                            >  2:323696/1‑274 (MQ=255)
                                                                                                                                                                                                                                                    aTACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCcac                                                                                                                                                                   >  3:72139/1‑140 (MQ=255)
                                                                                                                                                                                                                                                      actacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATa                                                                                                                                                        <  6:99768/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       ctacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATAc                                                                                                                                                       >  7:259826/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         acGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGAAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGa                                                                                                                                                     <  5:222936/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         acGGGTATAAAAAAACCCGCCGGGGCGGGCTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAATCCATCGGCGCGTTTTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCC‑‑ATCGAAACGCGGGCCCAGTTCGTTAAAACGTTTCGCCCCGATCTCGTTTTCCCTGTTTAGGGACAATGCCCGAAGCCGCTTCACAAAGCTATCAGCCTTGGCAag                      >  4:35444/1‑274 (MQ=255)
                                                                                                                                                                                                                                                           gggTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgt                                                                                                                                                   <  8:527636/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              tataAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTa                                                                                                                                                >  6:1280868/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 acaaaaaCCCGCCTGGTCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGCTCAACCAGCTCGTTGTAAGCCATCGGCGCGTTGTGGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCCCGGCT‑‑CGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGATGACGATTAGCaa                                         <  4:46455/245‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaaaCCCGCCGGGGCGGGGTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAacca                                                                                                                                            >  6:1367314/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCa                                                                                                                                                             >  7:400765/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCa                                                                                                                                                             <  8:400765/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccacc                                                                                                                                          >  6:1167465/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccacc                                                                                                                                          >  6:683359/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccacc                                                                                                                                          >  6:423021/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     aaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGAt                                                                                                                                                                                                         <  8:297193/85‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCg                                                                                                                                         <  7:804117/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaaCCCGCCCCGGCGGGTTTTTTTATGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgt                                                                                                                                                   <  6:1221837/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaaCCCGCCCCGGCGGGTTTTTTTATGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgt                                                                                                                                                   >  5:1221837/1‑139 (MQ=255)
                                                                                                                                                                                                                                                                     aaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgt                                                                                                                                                   <  7:706324/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCg                                                                                                                                         >  6:344328/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      aaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCgg                                                                                                                                        <  5:473865/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTaa                                                                                                                                               <  8:1072547/141‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTaa                                                                                                                                               >  7:1072547/1‑141 (MQ=255)
                                                                                                                                                                                                                                                                       aCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT‑‑CGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCt        >  2:176469/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                        cccGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTAc                                                                    >  4:22405/1‑215 (MQ=255)
                                                                                                                                                                                                                                                                         ccgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCAc                                                                                                                                     >  6:1141474/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            gccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATa                                                                                                                                                        <  2:48998/128‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT‑‑CGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGCGCCAGTTATCCCGAGTACGGGCGACAGCCAGACCACGCATAGCAACGGTCTCAGTTTTGGCGAACGTAAACGGCGGAGCAAcca  >  4:130997/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                 
CTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCT‑‑CGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTCAACTA  >  NC_000913/3439319‑3439860

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