Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,090,124 C→G 8.2% 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.2% 816.9 / 49.8 239intergenic (+78/‑68)hisL/hisGhis operon leader peptide/ATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (110/108);  new base G (9/11);  total (119/120)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

Only 100 of 243 total aligned reads displayed.

AGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAG  >  NC_000913/2089854‑2090393
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                             
            ccccTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAGACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCttt                                                                                                                                                                                                                                                               >  2:559499/1‑275 (MQ=255)
                          gCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATa                                                                                                                                                                                                                                             <  1:153402/279‑1 (MQ=255)
                                        cAAATCCAATGACAAAAAATATATAAAAAAGCCCTTGCTTTCTAACGTTGAAGTGGGTTAGGTTAAAAGAAATAAGTTGAATAAACATTCACAGAGACTTTTATGACCAGCGTTCAAATTAGACACCACCATCATCACCATCATTCTGGCTATTCGTTCAAGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa                                                                                                                                                                                                                                   <  4:559495/275‑1 (MQ=255)
                                           aTCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATACACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGACAGCCCCCGTAAGATCACCTCACGGGGGCTTTTTTATTCCTCGGTTGATAACTGTTCaaccag                                                                                                                                                                                                                                >  2:684531/1‑270 (MQ=255)
                                                          atatatAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa                                                                                                                                                                                                                                   <  2:227876/257‑1 (MQ=255)
                                                             tataAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGGGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa                                                                                                                                                                                                                                   <  1:130337/254‑1 (MQ=255)
                                                              ataAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGCTTAAAGAGGAATACCAAAATg                                                                                                                                                                                                         >  1:90235/1‑279 (MQ=255)
                                                                            gCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGacaaca                                                                                                                                                                                               <  2:600216/275‑1 (MQ=255)
                                                                                 atAACGTGAAAAGTGTTTAGGTTAAAGGATCTCATGTGAATGAACAGTGAACGGTACTATTAATACACGCGTTCAATTTAAACACTACCATCATCAACATCATCCTGACTTGTCCTGCCGGCGAGGTGTGGTGGAAGACCTTCAGATCTTCCAGTGGTGCATGAACGCATGGTAAAGCCCCCGGAAGATCACCTTCCGGGGGGTTTTTTATTGCGCGGTTGATAACTGTTCAGACAAGTTTAGAGAGGGATAACAAAATGACAGACAACACTCGt                                                                                                                                                                                          <  4:611270/274‑1 (MQ=255)
                                                                                                               aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAAACAACCCGCGTTTCCGAATAGCATTTGAGAAATCCGGCCGt                                                                                                                                                            >  4:583195/1‑275 (MQ=255)
                                                                                                                            aCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGTGGAATAACAAAATGACAGACAACACTCGTTTCCGCATAGCTATGCAGAAATCCGGCCGTTTAATGGATGACt                                                                                                                                               >  2:154554/1‑275 (MQ=255)
                                                                                                                                    agagACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCggg                                                                                                                                                                                                                                                                     <  7:668713/149‑1 (MQ=255)
                                                                                                                                    agagACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggg                                                                                                                                                                                                                                                                     >  5:169069/1‑149 (MQ=255)
                                                                                                                                     gagaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCgggg                                                                                                                                                                                                                                                                    <  7:461976/149‑1 (MQ=255)
                                                                                                                                     gagaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCgggg                                                                                                                                                                                                                                                                    <  6:730136/149‑1 (MQ=255)
                                                                                                                                     gagaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCgggg                                                                                                                                                                                                                                                                    <  5:205766/149‑1 (MQ=255)
                                                                                                                                        aCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGCCAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGc                                                                                                                                   >  4:649213/1‑275 (MQ=255)
                                                                                                                                          ttttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCttt                                                                                                                                                                                                                                                               >  7:252810/1‑149 (MQ=255)
                                                                                                                                            ttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCttttt                                                                                                                                                                                                                                                             <  5:736940/149‑1 (MQ=255)
                                                                                                                                               ttacacTCGTTTAATTTAAACACCACCATCATCACCTTCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAACACATTCAGCTCCACCAGTGGTGCATGAACGCATGAGAAAGCCCCAGGAAGATCAACTTCCGGGGGCTTTTTTAt                                                                                                                                                                                                                                                          <  6:169069/147‑1 (MQ=255)
                                                                                                                                                gACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAtt                                                                                                                                                                                                                                                         <  8:578978/149‑1 (MQ=255)
                                                                                                                                                gACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAtt                                                                                                                                                                                                                                                         <  8:252810/149‑1 (MQ=255)
                                                                                                                                                    cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcgc                                                                                                                                                                                                                                                     <  8:631250/149‑1 (MQ=255)
                                                                                                                                                    cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCAt                                                                                                                   >  1:268088/1‑279 (MQ=255)
                                                                                                                                                      cgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAt                                                                                                                                                                                                                                                          >  4:708556/1‑142 (MQ=255)
                                                                                                                                                        ttCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGtt                                                                                                                                                                                                                                                 <  5:1612145/149‑1 (MQ=255)
                                                                                                                                                        ttCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGtt                                                                                                                                                                                                                                                 <  6:207694/149‑1 (MQ=255)
                                                                                                                                                          cAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcg                                                                                                                                                                                                                                                      <  2:708561/142‑1 (MQ=255)
                                                                                                                                                          cAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGa                                                                                                                                                                                                                                               <  6:155960/149‑1 (MQ=255)
                                                                                                                                                          cAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACg                                                                                                                                                                                     <  2:419216/207‑1 (MQ=255)
                                                                                                                                                            aTTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATa                                                                                                                                                                                                                                             >  5:573161/1‑149 (MQ=255)
                                                                                                                                                             tttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCggg                                                                                                                                                                                                                                                                     >  8:807810/1‑124 (MQ=255)
                                                                                                                                                             tttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAAAAAGCCCCCGGAAGGTGATCTTCCgg                                                                                                                                                                                                                                                                      <  7:1203299/123‑1 (MQ=255)
                                                                                                                                                             tttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAAAAAGCCCCCGGAAGGTGATCTTCCgg                                                                                                                                                                                                                                                                      >  8:1203299/1‑123 (MQ=255)
                                                                                                                                                               tAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACg                                                                                                                                                                                                                                          <  6:67587/149‑1 (MQ=255)
                                                                                                                                                                  acacCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGtt                                                                                                                                                                                                                                       <  7:803370/149‑1 (MQ=255)
                                                                                                                                                                   caccacCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTc                                                                                                                                                                                                                                      <  5:536467/149‑1 (MQ=255)
                                                                                                                                                                    accaccATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCa                                                                                                                                                                                                                                     >  8:755944/1‑149 (MQ=255)
                                                                                                                                                                     ccaccaTCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAg                                                                                                                                                                                                                                    >  6:426915/1‑149 (MQ=255)
                                                                                                                                                                      caccatcatcaccatcatCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATAAAAAAGCCCCCGGAAGGTGATCTTCCggg                                                                                                                                                                                                                                                                     <  5:1353773/115‑1 (MQ=38)
                                                                                                                                                                      caccatcatcaccatcatCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATAAAAAAGCCCCCGGAAGGTGATCTTCCggg                                                                                                                                                                                                                                                                     >  6:1353773/1‑115 (MQ=38)
                                                                                                                                                                       accatcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGtt                                                                                                                                                                                                                                       >  2:280914/1‑144 (MQ=255)
                                                                                                                                                                         catcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACgg                                                                                                                                                                                                                                         >  8:526737/1‑140 (MQ=255)
                                                                                                                                                                         catcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACgg                                                                                                                                                                                                                                         <  7:526737/140‑1 (MQ=255)
                                                                                                                                                                         catcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAg                                                                                                                                                                                                                                <  8:1124184/149‑1 (MQ=255)
                                                                                                                                                                         catcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTcaca                                                                                                  <  2:391852/275‑1 (MQ=255)
                                                                                                                                                                          atcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAgg                                                                                                                                                                                                                               <  6:1414278/149‑1 (MQ=255)
                                                                                                                                                                           tcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa                                                                                                                                                                                                                                   <  4:280911/144‑1 (MQ=255)
                                                                                                                                                                           tcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGt                                                                                                                                                                                                                              <  7:339275/149‑1 (MQ=255)
                                                                                                                                                                                  catcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAg                                                                                                                                                                                                                                >  5:384416/1‑140 (MQ=255)
                                                                                                                                                                                  catcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAaga                                                                                                                                                                                                                       >  5:1186058/1‑149 (MQ=255)
                                                                                                                                                                                   atcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGAATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAagac                                                                                                                                                                                                                      >  5:1230857/1‑148 (MQ=255)
                                                                                                                                                                                    tcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggggg                                                                                                                                                                                                                                                                   <  2:441354/103‑1 (MQ=255)
                                                                                                                                                                                    tcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGAc                                                                                                                                                                                                                                  <  4:435164/136‑1 (MQ=255)
                                                                                                                                                                                         ctaactactCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCa                                                                                                                                                                                                                                     <  7:1424380/125‑1 (MQ=255)
                                                                                                                                                                                         cTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCGTTTTTATTGCGCGGTTGATAACGGTTCa                                                                                                                                                                                                                                     >  8:1424380/1‑128 (MQ=255)
                                                                                                                                                                                            aCTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGATAGGTTTAAAGAGGAATAACAAAATGACGGACAACACTCGTTTACGCATAGCTATGCAGAACTCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGACCCCGa                                                                                >  2:494617/1‑274 (MQ=255)
                                                                                                                                                                                              tAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACaaaa                                                                                                                                                                                                           <  6:989796/149‑1 (MQ=255)
                                                                                                                                                                                               aGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAAt                                                                                                                                                                                                          <  6:149006/149‑1 (MQ=255)
                                                                                                                                                                                               aGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAAt                                                                                                                                                                                                          >  6:560513/1‑149 (MQ=255)
                                                                                                                                                                                                   tttCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATgaca                                                                                                                                                                                                      <  7:496318/149‑1 (MQ=255)
                                                                                                                                                                                                        ggCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGacaa                                                                                                                                                                                                 <  8:1319860/149‑1 (MQ=255)
                                                                                                                                                                                                        ggCGATGTGTGCTGGAAGACAATCAGATCGTCCAGTGGTGCATGAACGCATGAGAGAGCCCCCGGAAGATCACCTTCCGAGGGCTATTTTATTGCGCGGTTGATATCGGTTCGGACAGGTTTAGAGAGGAATAACAAAATGACAGacaa                                                                                                                                                                                                 >  8:1048476/1‑149 (MQ=255)
                                                                                                                                                                                                             cgtgtgCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGCAAGCCCCCGGCAGAGCACCTTCCGGGCGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACATGTATAAAGAGGAATAACAAAATGACAGACAACACTc                                                                                                                                                                                            <  6:1230857/148‑1 (MQ=255)
                                                                                                                                                                                                                   cTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTaa                                                                                                                >  4:399170/1‑219 (MQ=255)
                                                                                                                                                                                                                      gAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGt                                                                                                                                                                                                                              <  5:1497984/106‑1 (MQ=255)
                                                                                                                                                                                                                       aaGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCTGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa                                                                                                                                                                                                                                   <  6:1472039/100‑1 (MQ=255)
                                                                                                                                                                                                                       aaGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACGTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa                                                                                                                                                                                                                                   >  5:1472039/1‑100 (MQ=255)
                                                                                                                                                                                                                            aTTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTa                                                                                                                                                                             >  7:1630152/1‑149 (MQ=255)
                                                                                                                                                                                                                             ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTAt                                                                                                                                                                            >  8:428670/1‑149 (MQ=255)
                                                                                                                                                                                                                              tCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATg                                                                                                                                                                           >  6:1053572/1‑149 (MQ=255)
                                                                                                                                                                                                                                aGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCa                                                                                                                                                                         <  8:1661638/149‑1 (MQ=255)
                                                                                                                                                                                                                                 gATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGacaa                                                                                                                                                                                                 >  6:1317608/1‑124 (MQ=255)
                                                                                                                                                                                                                                 gATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAg                                                                                                                                                                        >  5:651348/1‑149 (MQ=255)
                                                                                                                                                                                                                                  aTCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGa                                                                                                                                                                       >  6:283346/1‑149 (MQ=255)
                                                                                                                                                                                                                                        cAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGt                                                                                                                                                                                          >  7:1246175/1‑124 (MQ=255)
                                                                                                                                                                                                                                                cATGAAGGCTGGGGAAAGCCCCCTGATGGTGCNTNTTCGGGGGGTTTTTTATAGCGCTGTTGATAACGGGTCAGACAGGTNTAAAGAGGAATAACAAAAGGACAGACAACACGCGTTTACGCATAGCAATGCAAAAATCCGGGCGTTTAANTGATGACTCACGCGAATTGCTGGCACGCTGTGGAATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGCCAGAAAACATGCCGATGGATATTCTGCGCGTGCGTGACGACGACATTCCCGGTc                       <  3:90235/279‑1 (MQ=255)
                                                                                                                                                                                                                                                cATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGAAATTccc                           >  2:331232/1‑275 (MQ=255)
                                                                                                                                                                                                                                                cATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAAAATGCCGATTGATATTCTGCGCGGGCGTGACGACGACATTccc                           >  2:357032/1‑275 (MQ=255)
                                                                                                                                                                                                                                                 aTGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTaa                                                                                                                                                        <  6:891466/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTaa                                                                                                                                                        <  8:1618486/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   gAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGACATTCCCGGt                        <  4:195047/275‑1 (MQ=255)
                                                                                                                                                                                                                                                     aCGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGtga                                                                                                                                                    <  6:1282682/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        cATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGtgatga                                                                                                                                                 <  5:1074820/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         aTGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGcc                                                                                                                                                              >  8:875926/1‑135 (MQ=255)
                                                                                                                                                                                                                                                         aTGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGcc                                                                                                                                                              <  7:875926/135‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ccccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAAt                                                                                                                                      >  7:1370790/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    cccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAAtt                                                                                                                                     >  6:1541350/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     ccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACgg                                                                                                                                                                                                                                         >  5:1409220/1‑48 (MQ=255)
                                                                                                                                                                                                                                                                      cGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCAt                                                                                                                                                                                  >  7:1328288/1‑102 (MQ=255)
                                                                                                                                                                                                                                                                      cGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGtga                                                                                                                                                    <  6:102807/132‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGtga                                                                                                                                                    >  5:102807/1‑132 (MQ=255)
                                                                                                                                                                                                                                                                       ggAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCt                                                                                                                                  <  5:283346/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ggAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTCACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCt                                                                                                                                  <  7:755944/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGtt                                                                                                                                                                                                                                                 >  8:1359656/1‑37 (MQ=255)
                                                                                                                                                                                                                                                                        gAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTa   <  2:642463/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aTCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGgcgc                                                                                                                             >  5:972460/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            aTCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGgcgc                                                                                                                             >  7:206924/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                             tCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAagag                                                                                                                                                                                                                      >  6:419056/1‑59 (MQ=255)
                                                                                                                                                                                                                                                                             tCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGACATTCCCGGTc                       <  4:366077/250‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                             
AGTCCTCTTAGCCCCCTCCTTTCCCCGCTCATTCATTAAACAAATCCATTGCCATAAAATATATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAG  >  NC_000913/2089854‑2090393

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