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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,5880GC 6.2% 219.3 / 5.1 65intergenic (‑79/+28)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (36/24);  new base C (2/2);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-01

AAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTC  >  NC_000913/3439343‑3439855
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                           
aaGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGAtt                                                                                                                                                                                                                                                >  2:315709/1‑275 (MQ=255)
                    cAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACGCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCt                                                                                                                                                                                                                            <  2:420569/275‑1 (MQ=255)
                      aggtttACCCCTCCCGTAATTCCGTTGTGCCTTGTGATGCAGCATCCAGTTCCAATGTTTCGCCCCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                                                                                                                                                                                          <  4:418891/270‑1 (MQ=255)
                                                                   ccAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGNACAGCTATTGTCGATAAGTGGGGCTACTACGGGTATAAAAAAACCCGCCTGGGCGGGTTTTTTTACGTTNTTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAANCNTCTCTCTGTACGCCATccccccgtttt                                                                                                                                                                              >  3:37237/1‑264 (MQ=255)
                                                                     aGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCg                                                                                                                                                                            >  2:417936/1‑274 (MQ=255)
                                                                      gttacAATTGTTCTCCGATACTTGCCAGGTTATCAAACTCAACGTTCGCTGTCGCTTCTGCTAGATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCCGGGGATGAACAAAAAGAGTACAGCTAGTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGc                                                                                                                                                                           <  4:274135/270‑1 (MQ=255)
                                                                       tATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                                                                                                                                                                                          >  4:153934/1‑242 (MQ=255)
                                                                           aaTGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCt                                                                                                                                                                                                      <  2:153935/242‑1 (MQ=255)
                                                                            aTGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttttttt                                                                                                                                                                                                                                                               >  4:410443/1‑184 (MQ=255)
                                                                                ttCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGt                                                                                                                                                                                                                                                           <  2:410444/184‑1 (MQ=255)
                                                                                 tCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTGTATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTTGCCTGCCCGGAAGc                                                                                                                                                               >  1:77999/1‑275 (MQ=255)
                                                                                 tCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCCGCTGCACGGAAGc                                                                                                                                                               >  2:240351/1‑275 (MQ=255)
                                                                                 tCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGcc                                                                                                                                                              >  4:243793/1‑275 (MQ=255)
                                                                                  cgccgcTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCca                                                                                                                                                             <  2:243794/275‑1 (MQ=255)
                                                                                   tcccgcaCCTTCCAGTTTATCCAACTCACCTTTCGCTTGGGCTTCTGCTATATGGCCCTCTGACCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCca                                                                                                                                                             <  4:382627/269‑1 (MQ=255)
                                                                                   gccgcTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTAAGCTTACTCTGCAGCCGCTTCTGCTTTCTCTTAACGATCAACCATCTCGATGTAAGCCATCGGCTCGTTGTCGCCTGCACGGAAGCca                                                                                                                                                             >  2:517586/1‑275 (MQ=255)
                                                                                      gcTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCacac                                                                                                                                                          <  4:315708/275‑1 (MQ=255)
                                                                                                aGGTTATCAAACTCACCGTTCGCCTGCGCTTCTGCGATAAGACACGCTGAACACTGGTCAAGTAACCAAATAAATACTCCAGGTGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGTGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATAc                                                                                                                                                <  2:395172/275‑1 (MQ=255)
                                                                                                                          cgcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACAACCGGCACGGCTcgcg                                                                                                                      >  2:182440/1‑275 (MQ=255)
                                                                                                                                       tGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                                                                                                                                                                                            <  7:583304/112‑1 (MQ=255)
                                                                                                                                       tGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                                                                                                                                                                                            >  8:583304/1‑112 (MQ=255)
                                                                                                                                          cgcTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGTTCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACTGAAAGCCCACTCAAGAATACAAGTGTAACCACCGGCACGGCTCGCGCAAAGCGGACCCatat                                                                                                      >  4:237325/1‑272 (MQ=255)
                                                                                                                                            ctctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGcc                                                                                                                                                                          >  4:49985/1‑205 (MQ=255)
                                                                                                                                            ctctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTc                                                                                                     >  2:181167/1‑274 (MQ=255)
                                                                                                                                            ctctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTc                                                                                                     >  2:153107/1‑274 (MQ=255)
                                                                                                                                            ctctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCCCGTTTTCGCCTTCACGGACGCCACACTTCAGAATACGAGTTTCACCACCGGCACGCCTCGCGAACCGCGGGCCCACTTc                                                                                                     >  4:189759/1‑274 (MQ=255)
                                                                                                                                            ctctGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATCACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCAGCCCGGCTCGCGCAACGCGGGACCAATTc                                                                                                     >  1:9432/1‑274 (MQ=255)
                                                                                                                                                gCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCa                                                                                                                                                                      <  2:49986/205‑1 (MQ=255)
                                                                                                                                                      tGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTAc                                                                                                                                                                                                                                                             >  2:545394/1‑112 (MQ=255)
                                                                                                                                                          cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTg                                                                                                                                                                                                                                                         <  4:545393/112‑1 (MQ=255)
                                                                                                                                                           aaGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGt                                                                                                   >  1:108209/1‑261 (MQ=255)
                                                                                                                                                                aCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGAt                                                                                >  2:78837/1‑275 (MQ=255)
                                                                                                                                                                 cacattAAATACTCCAGGTGATGAACAAAAAGAGTACAGCGATTGTAGCTAAGTGGGGATACTACGGGTATAAACAACCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTTTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGAtc                                                                               <  2:295076/270‑1 (MQ=255)
                                                                                                                                                                  cacaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACCCTTCCGCATACGAGTGTACCCACCGGCACGGCTCCCGAAACGCCGGCCCAGTTCGTTAAACAGTTTTGCCCCAAtc                                                                               >  4:281657/1‑274 (MQ=255)
                                                                                                                                                                  cacaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGACTACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAAAAGTTTTGCCACGAtc                                                                               >  1:76405/1‑274 (MQ=255)
                                                                                                                                                                       aaaTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCggg                                                                                                                                                                                                                                                                      <  8:837401/86‑1 (MQ=255)
                                                                                                                                                                       aaaTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCggg                                                                                                                                                                                                                                                                      >  7:837401/1‑86 (MQ=255)
                                                                                                                                                                            cTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCAt                                                                                                                                                                                         >  2:550449/1‑158 (MQ=255)
                                                                                                                                                                            cTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCAc                                                                     >  2:220536/1‑274 (MQ=255)
                                                                                                                                                                            cTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCCGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCAc                                                                     >  1:84670/1‑274 (MQ=255)
                                                                                                                                                                              ccAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttt                                                                                                                                                                                                                                                                   <  7:778335/82‑1 (MQ=255)
                                                                                                                                                                              ccAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttt                                                                                                                                                                                                                                                                   >  8:778335/1‑82 (MQ=255)
                                                                                                                                                                               cAGGGGATGAACAAAAAGAGTACCGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGTTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGCAGCCACACTTCAGAATACGAGTGTAACCACCGGCCCGGCCCGCGACACCCGGGCCCATTTCCGTCAACAGGTTTGCCCAGCACTCGTCATCCCGAGt                                                                 >  1:130247/1‑275 (MQ=255)
                                                                                                                                                                               cAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGTTCAACCATCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGTCACCCACACTTCAGAATACGAGTGTAACCACCGGCCCGGCTCCCTAAACCCGCGCCCATGTCGTTTAACACGTTTTCCACGCACTCGTGATCAccac                                                                  >  4:354037/1‑271 (MQ=255)
                                                                                                                                                                                aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc                                                                                                                                                                                     <  4:550448/158‑1 (MQ=255)
                                                                                                                                                                                  tggaaaacaaaaaaGAGAGCAGGTGTTGTAGGTAAGTGGGGATTCTACGGGTATAAAAACACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGAGTACTCTGNAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTaaa                                                                                               <  3:108209/233‑1 (MQ=255)
                                                                                                                                                                                     aTGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCg                                                                                                                                                                                                                                                                        >  5:772808/1‑70 (MQ=255)
                                                                                                                                                                                     aTGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCg                                                                                                                                                                                                                                                                        <  6:772808/70‑1 (MQ=255)
                                                                                                                                                                                                             gTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGACAGACTACGATTAACaa                                    >  2:283117/1‑274 (MQ=255)
                                                                                                                                                                                                                  tAAGTGGGGATATTAGGGGTTTAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGGCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTAt                              <  4:182439/275‑1 (MQ=255)
                                                                                                                                                                                                                    aggtgtgATCCTACCGGTATAAACAACNNCGCCGGGTCGGGTTTTTTGAGGTTGCTTCGNATTACGCTGGCGCAGCTTCTTCTTTCTCTGAACGATCAACCAGCACGATGTAAGCCATCGGCGAGTTGTCGCCTNCACGGAAGCAACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTAt                              <  3:77999/267‑1 (MQ=255)
                                                                                                                                                                                                                       ggggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTAGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGGGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTc                         <  2:354038/275‑1 (MQ=255)
                                                                                                                                                                                                                        gggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCCCGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTACACCGTTTTGCCACGATCCCGTTATCCCCAGTACCGGGCCANNNCAGACGACGCTTACCaccaccacacacc                         >  3:152071/1‑262 (MQ=255)
                                                                                                                                                                                                                        gggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTGTCACGACTACGGGAGAATGCCAGACTACGATTAGCAAAGCTCTCaagct                        >  1:75684/1‑271 (MQ=255)
                                                                                                                                                                                                                               tacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACCACGATTCGCAACGCTATCAGTATTGGCaa                  >  1:70641/1‑274 (MQ=255)
                                                                                                                                                                                                                                acGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACCCTTCAGAATACGAGTTTAACCACCGGCACGGCCCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCCCGTTATCCCGAGTTCGGGCGAATGCCAGACGACAGTTAGCAACCCTATCAGTCTTGGccct                 >  4:215885/1‑271 (MQ=255)
                                                                                                                                                                                                                                 cGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTACAAACGCTATACCTTTTGGCAagag               >  1:110710/1‑275 (MQ=255)
                                                                                                                                                                                                                                  gggTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAga                                                 <  4:86955/240‑1 (MQ=255)
                                                                                                                                                                                                                                  gggTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAga                                                 >  2:86956/1‑240 (MQ=255)
                                                                                                                                                                                                                                   tgTATAAAAAATCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAAcg                                                                                                                 <  4:166679/174‑1 (MQ=255)
                                                                                                                                                                                                                                   ggTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCCAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGCTCTCGTTATCACGAGTACGGGCGAATGACAGACGACGATTAGCCACGCTATCCGTCTTGACAAGGGt              >  4:186468/1‑274 (MQ=255)
                                                                                                                                                                                                                                    gTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTAt                              <  2:215886/257‑1 (MQ=255)
                                                                                                                                                                                                                                       taAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATc          <  4:181166/274‑1 (MQ=255)
                                                                                                                                                                                                                                             acccggccgGTGCGTGTTTGTTTACGTTGCTTCAGATTACTCTGCAACAGCTTGTGCTTTCTCTGAACGATCAACCAGCTCTATGTAAGCCATCGGCGCGTTGTCGCCTGCACGAAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCt   <  2:435383/270‑1 (MQ=255)
                                                                                                                                                                                                                                             accccGCCGGGGTGGTGTTTTTTACTGTTNNCCAGATTCCTCTGCAGAAGCTGCTGATTCCNCGGAACGACCATCCAGATCGATGTAAGCCATCGGTGGGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGGNTGTAAACACCGGCACGGCTCTCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCt   <  3:75684/273‑1 (MQ=255)
                                                                                                                                                                                                                                              aCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTc  >  2:455082/1‑275 (MQ=255)
                                                                                                                                                                                                                                               cccGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGa                                                                                                                                                                                                                    >  7:542075/1‑64 (MQ=255)
                                                                                                                                                                                                                                               cccGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGa                                                                                                                                                                                                                    <  8:542075/64‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                           
AAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTCGCGAAACGCGGGCCCAGTTCGTTAAACAGTTTTGCCACGATCTCGTTATCACGAGTACGGGCGAATGCCAGACGACGATTAGCAACGCTATCAGTCTTGGCAAGAGTAATCAGCGGCTC  >  NC_000913/3439343‑3439855

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: