Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,090,126 C→T 8.7% intergenic (+80/‑66) 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,1260CT 8.7% 237.7 / 11.6 81intergenic (+80/‑66)hisL/hisGhis operon leader peptide/ATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (40/34);  new base T (5/2);  total (45/36)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-02

CTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTG  >  NC_000913/2089860‑2090399
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                 
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                                                                                                                                                                                                                                                                           >  6:248100/1‑64 (MQ=25)
                                                                                                                                                                                                                    gACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCggg                                                                                                                                                                                                                                                                           <  5:248100/64‑1 (MQ=25)
                                                                                                                                                                                                                       aTTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTa                                                                                                                                                                                   >  5:51177/1‑149 (MQ=255)
                                                                                                                                                                                                                        ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACGTTCCGGGGGCTTTCTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTAACGCTTAGCTAt                                                                                                                                                                                  <  8:402366/149‑1 (MQ=255)
                                                                                                                                                                                                                        ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCAt                                                                                                                                                                                        <  6:90153/143‑1 (MQ=255)
                                                                                                                                                                                                                        ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCAt                                                                                                                                                                                        >  5:90153/1‑143 (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)
                                                                                                                                                                                                                                                                          aCCTTCCGCGGGCTTTTTTATCACGCGGTTTATAACGGCTCAGACAGGTTTAAAGAGGAATAACACAATGACAAACATCACCCGTTTACGCATAGCTATGCAGTAATCCGGCCGTTTAAGTGATGACTCACGCGAATTTCTGGCGCACTGTGGCATTAAAATTAATCTTCTCCCCAAGCACCTGATCGAGTTGGCAGACACAATTCACATTGATTTTCTGCGCGCGCGTTCCGACCACATTCCCGTCCTGGCAACGGAGCGCTTGGttaccctgg  >  4:384429/1‑267 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                 
CTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTG  >  NC_000913/2089860‑2090399

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