Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,090,124 C→G 8.9% intergenic (+78/‑68) hisL → / → hisG his operon leader peptide/ATP phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,090,1240CG 8.9% 256.1 / 15.4 79intergenic (+78/‑68)hisL/hisGhis operon leader peptide/ATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (38/33);  new base G (5/2);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-03

CTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGA  >  NC_000913/2089860‑2090382
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                  
ctTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggg                                                                                                                                                                                                                                                          <  2:342965/275‑1 (MQ=255)
        cctcctTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTc                                                                                                                                                                                                                                                              >  4:176888/1‑263 (MQ=255)
            ctTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggg                                                                                                                                                                                                                                                          <  2:176890/263‑1 (MQ=255)
                                  cAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa                                                                                                                                                                                                                        <  1:67671/275‑1 (MQ=255)
                                                 aaaaTATAATAAAAAAGCCTTTGCTGTAANGCGGAAAAGTGGGTTAGGTTAAGAANCATAAGTTGAAAAAACATTCACAGAAACTTTTATGACACGGGTTACATTTAAACACCACCATCATCACCATCATCCTGACNAGTCTTTCAGGCGATGTGTTCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAgg                                                                                                                                                                                                          <  3:90564/275‑1 (MQ=255)
                                                       tataAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACa                                                                                                                                                                                                   >  2:413126/1‑275 (MQ=255)
                                                                       gCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTACCAGGGGTGCATGAACGCATGAGAATGCCCCCGGAAGATCACCTTCCGGGGGCTTTCTGCTGGCGCGGTTg                                                                                                                                                                                                                                     >  3:9619/1‑226 (MQ=255)
                                                                           tCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATaa                                                                                                                                                                                                                                 <  1:9619/226‑1 (MQ=255)
                                                                                  tGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGc                                                                                                                                                                         >  4:103656/1‑275 (MQ=255)
                                                                                   gAAAGTGGTTTAGGTTAAAAAAAATCCTTTGAATAAACATTCCCAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCCCCATCATCCTGACTAGTCTTTCAGGCGATGTGGGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGAGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCa                                                                                                                                                                        <  4:306302/275‑1 (MQ=255)
                                                                                           tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTa                                                                                                                                                                            >  2:306306/1‑263 (MQ=255)
                                                                                           tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATg                                                                                                                                                                >  2:455497/1‑275 (MQ=255)
                                                                                                                         aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCCCGGTTCATAACGGTTCAGACCGGTTTAAAGAGGAATAACAAAATGACAGACACAACTCTTTTACGCATAGCAATGCACAAATCCGGCCGTTTAAGTGTTTACTcc                                                                                                                                  >  4:5778/1‑274 (MQ=255)
                                                                                                                             acaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGtt                                                                                                                                                                                                                                      >  4:309344/1‑171 (MQ=255)
                                                                                                                                 agaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATa                                                                                                                                                                                                                                  <  2:309348/171‑1 (MQ=255)
                                                                                                                                   aCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAAGCGGGCCGTTTAAGGGATGACTCATGCGAAttac                                                                                                                        >  2:356570/1‑273 (MQ=255)
                                                                                                                                     ttttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgg                                                                                                                                                                                                                                                           <  8:24634/142‑1 (MQ=255)
                                                                                                                                     ttttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgg                                                                                                                                                                                                                                                           >  7:24634/1‑142 (MQ=255)
                                                                                                                                         aTGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAacac                                                                                                                                                                                   >  2:453402/1‑210 (MQ=255)
                                                                                                                                             cacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcg                                                                                                                                                                                                                                           <  2:68571/150‑1 (MQ=255)
                                                                                                                                             cacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGt                                                                                                                                                                               <  4:453397/210‑1 (MQ=255)
                                                                                                                                             cacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTa                                                                                                                                                                            >  3:143427/1‑213 (MQ=255)
                                                                                                                                             cacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCtg                                                                                                              >  1:59305/1‑275 (MQ=255)
                                                                                                                                                 cgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCa                                                                                                                                                                        <  1:143427/213‑1 (MQ=255)
                                                                                                                                                  gTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGt                                                                                                                                                                                                                                       >  7:427060/1‑149 (MQ=255)
                                                                                                                                                   ttCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCAt                                                                                                        <  2:203162/275‑1 (MQ=255)
                                                                                                                                                              caccacCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTc                                                                                                                                                                                                                           <  8:427060/149‑1 (MQ=255)
                                                                                                                                                                  accatcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGCCTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTc                                                                                          >  2:267307/1‑274 (MQ=255)
                                                                                                                                                                   ccatcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACAATCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGAGAATTGCTGGCGCGCTGGGGCAATAAAATAAATATTcac                                                                                        >  1:56837/1‑275 (MQ=255)
                                                                                                                                                                      gcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACAcc                                                                                     <  2:5778/274‑1 (MQ=255)
                                                                                                                                                                                  tcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCTCTGTGGCATTAAAATTAATCTTCCCACCCAGCGCCTCAt                                                                          >  1:58759/1‑274 (MQ=255)
                                                                                                                                                                                   ccTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATcc                                                                        >  2:261675/1‑274 (MQ=255)
                                                                                                                                                                                     tGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgg                                                                                                                                                                                                                                                           >  7:388856/1‑94 (MQ=255)
                                                                                                                                                                                     tGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgg                                                                                                                                                                                                                                                           <  8:388856/94‑1 (MQ=255)
                                                                                                                                                                                      gagtatgcTTTCAGGCGATTTTGGCTAGANNACATTCAGATCTTCCATTGGTGCAGTAACGNATTAGAAATACTCCGGAATATCACCTTCCGGGGTCTGTTTTATTGCGAGGTTGATAACGGTTCAGACAGGTTTANAGAGGAATAACAAAAGGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGa                                                                     <  3:59305/268‑1 (MQ=255)
                                                                                                                                                                                           gTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACaaa                                                                                                                                                                                                 >  8:184539/1‑146 (MQ=255)
                                                                                                                                                                                           gTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACaaa                                                                                                                                                                                                 <  7:184539/146‑1 (MQ=255)
                                                                                                                                                                                           gTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCCTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGCATAACAAAATGACAGACAACACTCGTTTACGCAt                                                                                                                                                                       >  2:145463/1‑172 (MQ=255)
                                                                                                                                                                                            tCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATga                                                                                                                                                                                             >  7:285689/1‑149 (MQ=255)
                                                                                                                                                                                                  aGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATgacagaca                                                                                                                                                                                       >  7:402366/1‑149 (MQ=255)
                                                                                                                                                                                                    ggggtgttggCTGGAAGACATTCCAATCTTCCAGGGGGGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATTACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATg                                                       <  4:261671/266‑1 (MQ=255)
                                                                                                                                                                                                     gcagtgtggCTGGAGGACATTCAGATCTGCCAGTGGGGCATGAACGCATGGGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATGGCGCGGTTTATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATg                                                       <  4:455492/266‑1 (MQ=255)
                                                                                                                                                                                                      gATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATGGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGcc                                                     <  2:103657/275‑1 (MQ=255)
                                                                                                                                                                                                             gCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTa                                                                                                                                                                            <  8:285689/149‑1 (MQ=255)
                                                                                                                                                                                                               tGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACg                                                                                                                                                                          >  7:107816/1‑149 (MQ=255)
                                                                                                                                                                                                                  aaGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAaga                                                                                                                                                                                                            <  5:227109/112‑1 (MQ=255)
                                                                                                                                                                                                                  aaGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCCCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCt                                         >  4:360797/1‑275 (MQ=255)
                                                                                                                                                                                                                  aaGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGGTTAAaga                                                                                                                                                                                                            >  6:227109/1‑112 (MQ=255)
                                                                                                                                                                                                                    gACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCggg                                                                                                                                                                                                                                                          <  5:248100/64‑1 (MQ=25)
                                                                                                                                                                                                                    gACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCggg                                                                                                                                                                                                                                                          >  6:248100/1‑64 (MQ=25)
                                                                                                                                                                                                                       aTTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTa                                                                                                                                                                  >  5:51177/1‑149 (MQ=255)
                                                                                                                                                                                                                        ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACGTTCCGGGGGCTTTCTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTAACGCTTAGCTAt                                                                                                                                                                 <  8:402366/149‑1 (MQ=255)
                                                                                                                                                                                                                        ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCAt                                                                                                                                                                       >  5:90153/1‑143 (MQ=255)
                                                                                                                                                                                                                        ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCAt                                                                                                                                                                       <  6:90153/143‑1 (MQ=255)
                                                                                                                                                                                                                        ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTAt                                                                                                                                                                 >  7:453150/1‑149 (MQ=255)
                                                                                                                                                                                                                                 tCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGATATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCCGGCGCGCTGTGGCATTCAACTCAATCTTCACACACAACGCCTGATCGCGATGGCAAAAAACATGCCGATTAATACTCTGCgcgtggatgaaga                          >  4:447217/1‑267 (MQ=255)
                                                                                                                                                                                                                                    aGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTgacgac                         <  2:360802/273‑1 (MQ=255)
                                                                                                                                                                                                                                          gCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTa                                                                                                                                                                            >  2:120457/1‑120 (MQ=255)
                                                                                                                                                                                                                                           cATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCAGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGAATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACGTGCCGATTGAAATTCGGCACATGCGTGCCGCCGACATTccc                >  2:330665/1‑275 (MQ=255)
                                                                                                                                                                                                                                              gAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCa                                                                                                                                                                        <  4:120456/120‑1 (MQ=255)
                                                                                                                                                                                                                                               aaaccATGAAAAAACCCCCGGAAGAACACNNCCAGGGGGGTTTTTTATTTCACGTTTAATANCGGTTAAAAAAGGGTTAAAAAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCNGCCGTTTAAGTGATGACGCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGACATTCCCGGTc            <  3:58759/271‑1 (MQ=255)
                                                                                                                                                                                                                                               aaCGAATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGACATTCCCGGTc            <  2:214603/275‑1 (MQ=255)
                                                                                                                                                                                                                                                aCGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGtga                                                                                                                                         <  8:107816/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   cATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGTTAACGGTTCAGACAGGTTTAACGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGAAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGCCGACGACATTCCCGGTCTGGt        >  2:382448/1‑275 (MQ=255)
                                                                                                                                                                                                                                                   cATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGTTGACTCACGCGAATTGCTGGCGCGCTGTGGCCTTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACa                                                         >  2:16150/1‑226 (MQ=255)
                                                                                                                                                                                                                                                   cATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACa                                                         >  4:103360/1‑226 (MQ=255)
                                                                                                                                                                                                                                                       agaAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGcc                                                     <  2:103361/226‑1 (MQ=255)
                                                                                                                                                                                                                                                          aaGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCAcgc                                                                                                                               <  8:453150/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           aGCACCAGGAAGGTAACCTTCCGGGGAGTTTTTTATTGCGGGGTGGGAAACGGGTAAGACAGTTGTAAAGGTGAATAAACATATGACAGACAACGCCCGTTTACGTATAGCTAGGCAGAAATCCGGCCGTTTAAGTGATGACACCCGCGATGTGCTGGCGCGCTGTAGCATTAGAATGAATCTTCACGCCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGGCGACGACATTCCCGGTCTGGTAatgga  <  4:267303/273‑1 (MQ=255)
                                                                                                                                                                                                                                                             cccccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCtg                                                                                                              <  1:60866/163‑1 (MQ=255)
                                                                                                                                                                                                                                                              ccccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATa                                                                                                                                                                      >  6:106626/1‑106 (MQ=255)
                                                                                                                                                                                                                                                              ccccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATa                                                                                                                                                                      <  5:106626/106‑1 (MQ=39)
                                                                                                                                                                                                                                                              ccccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGgcgcg                                                                                                                 >  1:133777/1‑159 (MQ=255)
                                                                                                                                                                                                                                                               cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAAtt                                                                                                                          >  7:182338/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAAtt                                                                                                                          >  5:99/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                ccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTg                                                                                                                         >  8:468497/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 cGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGc                                                                                                                        >  5:424625/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  ggAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCtgt                                                                                                             <  3:133777/159‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aTCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGgcgc                                                                                                                  >  6:421520/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         cACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGgcgcgc                                                                                                                >  8:426901/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                  
CTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGA  >  NC_000913/2089860‑2090382

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