Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,191,239 T→C 5.4% K136R (AAG→AGG)  yehC ← putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,191,2390TC 5.4% 227.8 / 2.4 75K136R (AAG→AGG) yehCputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (31/39);  new base C (2/2);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGATTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGC  >  NC_000913/2191094‑2191382
                                                                                                                                                 |                                                                                                                                               
aTTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTAc                                                                                                                                              <  7:95840/149‑1 (MQ=255)
     cTGAGTCATTTTTTATCACCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACcttct                                                                                                                                         <  7:308915/149‑1 (MQ=255)
        aGTCATTTTTTATCACCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGc                                                                                                                                      >  8:225274/1‑149 (MQ=255)
            aTTTTTTATCCCCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGctct                                                                                                                                   <  7:225274/148‑1 (MQ=255)
            aTTTTTTATCACCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTg                                                                                                                                  >  5:182824/1‑149 (MQ=255)
             ttttttATCACCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTgg                                                                                                                                 <  3:233607/149‑1 (MQ=255)
             ttttttATCACCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTg                                                                                                                                  >  3:117321/1‑148 (MQ=255)
                   aTCACCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTAtt                                                                                                                           >  7:131899/1‑149 (MQ=255)
                    tCACCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTAtt                                                                                                                           <  2:51836/148‑1 (MQ=255)
                      aCCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTtgg                                                                                                                        >  2:81952/1‑149 (MQ=255)
                      aCCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTtg                                                                                                                         <  8:104365/148‑1 (MQ=255)
                      aCCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTAACTTCTTGCTCTGGACTATTtgg                                                                                                                        >  5:21478/1‑149 (MQ=255)
                       ccAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTtggt                                                                                                                       >  1:16232/1‑149 (MQ=255)
                       ccAAACCATTGCCACTGCGATTTACCAGAAATTTTTTAAATGTCGCTTTATTTACGGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTtggt                                                                                                                       <  5:18684/149‑1 (MQ=255)
                        caacaccttGCCACTGAGATTTACCAGCACTTTTTTAAATGTCGCCTTATTTCCCGGAGCAATCCCCGCTGGCCGGGAAAACAACTTAATTCTGTTTTGCAGCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTtggtg                                                                                                                      <  2:210833/142‑1 (MQ=255)
                         aaaCCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTtggtgg                                                                                                                     <  6:151890/149‑1 (MQ=255)
                           aCCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGaa                                                                                                                   <  1:283601/149‑1 (MQ=255)
                               ttGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATc                                                                                                               <  7:229091/149‑1 (MQ=255)
                               ttGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATc                                                                                                               <  5:296984/149‑1 (MQ=255)
                               ttGCCACTGAGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCAGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGAATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATc                                                                                                               >  6:310774/1‑149 (MQ=255)
                                  ccACTGAGAATTACCCTCAATATTTTAAATGTCTCTTTAATTACCTGAGCAATACCCGCTGGACGGTAAAAAAACTTAATTCTGATTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTATGGACTATTTGGTGGAATATCCAg                                                                                                            <  1:260542/149‑1 (MQ=255)
                                   cACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAg                                                                                                            >  3:48079/1‑148 (MQ=255)
                                    aCTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGaa                                                                                                          >  2:9550/1‑149 (MQ=255)
                                      tGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACa                                                                                                        >  3:48542/1‑149 (MQ=255)
                                      tGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACa                                                                                                        >  4:176837/1‑149 (MQ=255)
                                        cGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACAtt                                                                                                      <  1:152148/149‑1 (MQ=255)
                                         gATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTc                                                                                                     <  8:131899/149‑1 (MQ=255)
                                         gATTTACCAGCAATTTTTTAAATGGCGCTTTATTTACCTGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTGTTACCTTCTTGCTCTGGACTATTTTGTGGAATATACAGAACATTc                                                                                                     <  2:189036/149‑1 (MQ=255)
                                          aTTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGAAAACTTCAGTGCATTCTTACATTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACAtt                                                                                                      >  8:301826/1‑147 (MQ=255)
                                           tttACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAg                                                                                                   <  1:81952/149‑1 (MQ=255)
                                            ttACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAAATTAATTCTGTTTTGCATCGAAAACTTAAGTGCATTCATACCTTCTTGCTCTCGAATATTTCGTGGAATATCCAGAACATTCAg                                                                                                   >  3:163406/1‑148 (MQ=255)
                                               ccAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATa                                                                                                >  7:311501/1‑148 (MQ=255)
                                                cAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATaa                                                                                               >  4:141996/1‑148 (MQ=255)
                                                cAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATaa                                                                                               >  5:299087/1‑148 (MQ=255)
                                                 aGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATaaa                                                                                              >  4:281818/1‑148 (MQ=255)
                                                   cAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATaaaaaa                                                                                           <  6:196782/149‑1 (MQ=255)
                                                    aaTTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCCTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAAt                                                                                          >  6:257161/1‑149 (MQ=255)
                                                    aaTTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCCTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAAt                                                                                          >  6:252786/1‑149 (MQ=255)
                                                       tttttAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATgat                                                                                       >  8:185672/1‑147 (MQ=255)
                                                         tttAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTAGctc                                                                                                                                    <  1:155140/102‑1 (MQ=255)
                                                         tttAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTAGctc                                                                                                                                    >  2:155140/1‑102 (MQ=255)
                                                           tAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGctttc                                                                                    <  1:64489/148‑1 (MQ=255)
                                                             aaTGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAATCAACTTAATTCTGTTTTGCATCGGAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGctttctt                                                                                  <  6:247201/148‑1 (MQ=255)
                                                             aaTGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGctttctt                                                                                  <  2:52278/148‑1 (MQ=255)
                                                             aaTGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCCTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGctttctt                                                                                  <  5:252786/148‑1 (MQ=255)
                                                             aaTGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCCTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGctttctt                                                                                  <  5:257161/148‑1 (MQ=255)
                                                               tGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTttat                                                                               <  3:281818/149‑1 (MQ=255)
                                                                  cGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAg                                                                            >  8:252913/1‑149 (MQ=255)
                                                                   gCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGt                                                                           >  4:231612/1‑149 (MQ=255)
                                                                       tatttaCCTGATCAATACCCGCTGGCCGTTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTggg                                                                        <  4:163406/148‑1 (MQ=255)
                                                                         tttaCCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTgg                                                                         <  5:117343/145‑1 (MQ=255)
                                                                         tttaCCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTgg                                                                         >  6:117343/1‑145 (MQ=255)
                                                                            aCCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGtt                                                                   >  8:323002/1‑148 (MQ=255)
                                                                            aCCGGAGCAAAACCCGCTCGCCGGTAAAACAACTTAATTCTGTTTTGAATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGCCTATTTGGTGGAATATCCAGCACATTCAGATAAAAAATGCTTTCTTTAGTAGTGCGCAGtt                                                                   >  4:267740/1‑148 (MQ=255)
                                                                             ccGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCCTCGCCAACTTCAGTGCATTCTTACCTTCTTGCTCTGGCCTATTTGGTCGAATATCCAGAACATTCAGATAAAAAATGATTTCTTTATTAGTGGGCAGttaa                                                                 >  6:100713/1‑147 (MQ=255)
                                                                             ccGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTa                                                                 <  3:181979/149‑1 (MQ=255)
                                                                                gAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTAAGATAAAAAAAGCTTTCTTTATTAGTGGGCAGTTTATTc                                                              >  8:292090/1‑149 (MQ=255)
                                                                                 aGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCg                                                             <  2:16232/149‑1 (MQ=255)
                                                                                 aGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCg                                                             <  7:252913/149‑1 (MQ=255)
                                                                                       aCCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATaa                                                       >  4:57205/1‑149 (MQ=255)
                                                                                           gCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAAtttt                                                   >  2:255275/1‑149 (MQ=255)
                                                                                               gCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGAt                                                <  4:117321/148‑1 (MQ=255)
                                                                                               gCCGGTAAAAAAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTAACTTCTTGCTCTGGACTATTTGGTTGAATATCCATAACATTCATATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGAt                                                <  5:310774/148‑1 (MQ=255)
                                                                                                       aacaacTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTTGAATATACATAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGATTTTTACTtg                                       <  7:301826/149‑1 (MQ=255)
                                                                                                               aaTTCTGTTTTGCATCGCAATCTTCAGTGCATTCTTACCTTCTTGCGCTGGACTATTTGGTGGCATATCCAGAACATTCAGAGAAAAAATGCTTTCTTTATTAGTGGGCTGTTTATTCGGCAGAATTTTGATTTTTACTTGCTGCCCg                                >  2:245662/1‑148 (MQ=255)
                                                                                                               aaTTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGATTTTTACTTGCTGccc                                 <  6:299087/147‑1 (MQ=255)
                                                                                                                aTTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGATTTTTACTTGCTGCCCgg                               >  6:37888/1‑148 (MQ=255)
                                                                                                                 ttCTGTTTTGCATAGCAAACTTCAGTGCATTCTTACCTTCTTTCTCTTGTCTATTTGGTTGAATATCCAGAACATTCAGATAAAAAATTCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGATTTTTACTTGCTGCCCGGaa                             <  7:185672/149‑1 (MQ=255)
                                                                                                                      ttttGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGATTTTTACTTGCTGCCCGGAATTTg                         <  3:231612/148‑1 (MQ=255)
                                                                                                                      ttttGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGATTTTTACTTGCTGCCCGGAATTTg                         <  7:155440/148‑1 (MQ=255)
                                                                                                                             tCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGATTTTTACTTGCTGCCCGGAATTTGCCCCTAtt                 >  2:297447/1‑149 (MQ=255)
                                                                                                                                       cAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGATTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCAc         <  3:57205/147‑1 (MQ=255)
                                                                                                                                          tGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATCTCCAGAACATTCAGATACAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGATTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACtggtg    <  4:48079/149‑1 (MQ=255)
                                                                                                                                            cATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGATTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACtggtgg   <  3:141996/148‑1 (MQ=255)
                                                                                                                                            cATTCTTAACTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTAATTAGTGCGCAGTTTATTCGGCATAATTTTGATTATTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGc  <  7:292090/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ATTAGCTGAGTCATTTTTTATCACCAAACCATTGCCACTGCGATTTACCAGCAATTTTTTAAATGTCGCTTTATTTACCGGAGCAATACCCGCTGGCCGGTAAAACAACTTAATTCTGTTTTGCATCGCAAACTTCAGTGCATTCTTACCTTCTTGCTCTGGACTATTTGGTGGAATATCCAGAACATTCAGATAAAAAATGCTTTCTTTATTAGTGGGCAGTTTATTCGGCATAATTTTGATTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGC  >  NC_000913/2191094‑2191382

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