Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 155,248 A→G 5.0% I60T (ATA→ACA)  htrE ← putative outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913155,2480AG 5.0% 231.0 / 2.6 79I60T (ATA→ACA) htrEputative outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base A (35/40);  new base G (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

CGGGAGAATTAATATGAAACTGCAATAAATTCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTC  >  NC_000913/155102‑155391
                                                                                                                                                  |                                                                                                                                               
cGGGAGAATTAATATGAAACTGCAATAAATTCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGtata                                                                                                                                               <  1:1012974/149‑1 (MQ=255)
        ttAATATGAAACTGCAATAAATTCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGGTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGtt                                                                                                                                       <  1:541224/149‑1 (MQ=255)
             atGAAACTGCAATAAATTCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCtt                                                                                                                                  <  2:287085/149‑1 (MQ=255)
               gAAACTGCAATAAATTCTTTAATGTTATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTTGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTAGTAAACACCCGGTATAGCGGGGTTACCTTCa                                                                                                                                <  2:631263/149‑1 (MQ=255)
               gAAACTGCAATAAATTCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCa                                                                                                                                >  1:283165/1‑149 (MQ=255)
                  aCTGCAATAAATTCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGTTATAGCGGGGTTACATTCACTg                                                                                                                             >  1:989408/1‑149 (MQ=255)
                  aCTGCAATAAATTCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTg                                                                                                                             <  2:1021284/149‑1 (MQ=255)
                      cAATAAATTCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTAcc                                                                                                                         >  2:119739/1‑149 (MQ=255)
                          aaaTTCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCt                                                                                                                     <  2:469585/149‑1 (MQ=255)
                           aaTTCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTg                                                                                                                    <  1:207434/149‑1 (MQ=255)
                           aaTTCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACAGCCGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTg                                                                                                                    >  1:1070034/1‑149 (MQ=255)
                            aTTCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTga                                                                                                                   >  2:37086/1‑149 (MQ=255)
                              tCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTgaga                                                                                                                 >  2:110817/1‑149 (MQ=255)
                                tttAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCATTTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACAGTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCCCGGTCCCGGCTGCGATc                                                                                                               >  2:705098/1‑149 (MQ=255)
                                tttAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATc                                                                                                               >  1:792257/1‑149 (MQ=255)
                                 ttAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGT‑AACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCaa                                                                                                             <  1:559248/149‑1 (MQ=255)
                                    aTGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAata                                                                                                           >  1:6809/1‑149 (MQ=255)
                                     tgtgATACAAGCCTGGGCGTTGTGTTTTCCTTCAATTGCGACAAATGTAATACTTTGGGTAATGATTGGTTGGTCGTTTACATAAACAATGCCGTCGTATCCACCCGGTATAGCGGGGTTACCTGAACTGTACAGGCTGAGATGAttat                                                                                                          >  2:661694/1‑149 (MQ=255)
                                     tgtgATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAatat                                                                                                          >  2:823569/1‑149 (MQ=255)
                                     tgtgATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAatat                                                                                                          <  1:431385/149‑1 (MQ=255)
                                      gtgATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGAAGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATAtt                                                                                                         >  1:453424/1‑149 (MQ=255)
                                       tgATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAACTGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTa                                                                                                        <  1:1124204/149‑1 (MQ=255)
                                           aCAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATg                                                                                                    >  2:507370/1‑149 (MQ=255)
                                           aCAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTACATg                                                                                                    >  1:537083/1‑149 (MQ=255)
                                                  tGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTg                                                                                             <  2:103493/149‑1 (MQ=255)
                                                   gggCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGc                                                                                            >  2:349413/1‑149 (MQ=255)
                                                   gggCGTTCTTTTTTCCTTAAAATGCGACAAATCTAATACTTTGGTTAATGCTTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGTTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGc                                                                                            <  2:136042/149‑1 (MQ=255)
                                                     gCGTTCTTTTATCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGccc                                                                                          <  2:699897/149‑1 (MQ=255)
                                                        ttCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCAt                                                                                       <  2:291967/149‑1 (MQ=255)
                                                         tCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGTGATCAATATTAGATGCATTCTGCCCCATc                                                                                      <  1:641161/149‑1 (MQ=255)
                                                            tttttCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAgg                                                                                   <  1:415349/149‑1 (MQ=255)
                                                                  cTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTa                                                                             <  1:936205/149‑1 (MQ=255)
                                                                    tCAATTGCGACAAATGTAATTCTTTGGTTATTGATTGGTTGGTCGTTTTCATAAACTCTGACGTCGTAAACACCCGGTATTGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATg                                                                           <  1:149506/149‑1 (MQ=255)
                                                                         tGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCAGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCat                                                                      <  1:340147/149‑1 (MQ=255)
                                                                               aaaTGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCaa                                                                <  1:539716/149‑1 (MQ=255)
                                                                                aaTGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGTGTTACCTTCAGTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAAc                                                               <  1:679909/149‑1 (MQ=255)
                                                                                aaTGTAATACTTTGGTTAATGATTGGGTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATAATCAAc                                                               <  2:1021203/149‑1 (MQ=255)
                                                                                    tAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGt                                                           >  1:859610/1‑149 (MQ=255)
                                                                                    tAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCACATCAGGAAGGTATGGTCTTATTCAACGCGt                                                           >  2:897850/1‑149 (MQ=255)
                                                                                       tACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCa                                                        <  2:797675/149‑1 (MQ=255)
                                                                                        acattGGTTAATGATTGGTGTGTCGTTTACATTAACTCTGACGGCGTTTACAGCCGGTATAGCGGTGTTACCTTCACTGTACCGGCTGAGAGCATTATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAg                                                       <  2:940700/146‑1 (MQ=255)
                                                                                         cTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGc                                                      <  2:246514/149‑1 (MQ=255)
                                                                                         cTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGc                                                      <  1:466079/149‑1 (MQ=255)
                                                                                          tttGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTGTAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCa                                                     >  1:720525/1‑149 (MQ=255)
                                                                                          tttGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCa                                                     >  2:209628/1‑149 (MQ=255)
                                                                                          tttGGTTAATGATTGGTTGGTCGGTTACATAAACACTGACGTCGTAAACACCCGGTGTAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCa                                                     >  1:690716/1‑149 (MQ=255)
                                                                                           ttGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCAc                                                    >  1:684726/1‑149 (MQ=255)
                                                                                             ggTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTg                                                  <  2:792257/149‑1 (MQ=255)
                                                                                              tttAATGATTGGTTGGTCGTTTAAATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGa                                                 <  1:774834/148‑1 (MQ=255)
                                                                                              gTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACGCGGTATAGCGGCGTTACCTACACTGTACCGGATGACATCAATATAAGATGCATTATGACCCATCAGTAAGTTAATGTCATATTCAACTAGATAAGCaatta                                                 >  1:212481/1‑145 (MQ=255)
                                                                                               ttAATGATTGGTTGGTCGTTTACATAAACACTTACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTCTACCGGCTGAGATCAATATTAGATGCATTCTGCAACATCAGGAAGGTATGGTCATATTCAACGAGTTAAGCACTGaa                                                >  1:588529/1‑149 (MQ=255)
                                                                                               ttAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGaa                                                >  2:150129/1‑149 (MQ=255)
                                                                                               ttAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATACCACGAGTTCAGCACTGaa                                                >  1:923591/1‑149 (MQ=255)
                                                                                                 aaTGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAg                                              >  1:1167631/1‑149 (MQ=255)
                                                                                                    gATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAg                                           <  1:357756/149‑1 (MQ=255)
                                                                                                    gATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAg                                           <  1:823569/149‑1 (MQ=255)
                                                                                                    gATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGGATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATTTTAGATGCATTCTGCCCCATCAGGAAGGTATGGGCATATTCAACGAGTTCAGCACTGAAAGCAg                                           <  2:1070034/149‑1 (MQ=255)
                                                                                                     aTTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCACAATAAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAATCAGt                                          >  1:1076338/1‑149 (MQ=255)
                                                                                                     aTTGGTTGGTCGATTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATGCTGCCCCGTCAGGAAGGGATGGTCATATTCATCGAGTTCAGCACTGAAAGCAGt                                          >  1:234006/1‑149 (MQ=255)
                                                                                                        gtgtggtCGTTTACATAAACACTTCAGTCGTAAACAACCGGTATAGCGGGTTTACCTTCACTGTACCGTCTGAGATAAATATTAGATGCATTCTTCCCCAGCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAACGCAGTAtt                                       <  1:67358/146‑1 (MQ=255)
                                                                                                        ggtgggtCGGTTACATAAACACTGACGTGGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCGCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTAtt                                       <  1:958069/149‑1 (MQ=255)
                                                                                                         gttggtCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCACCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTg                                      >  1:247589/1‑149 (MQ=255)
                                                                                                          ttggtCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATAAATAATAGATGCATTCTGCCGCCTCAGGAAGGTATGGTCATAATCAACGAGTTCACCACTGCAAGCAGTATTga                                     >  2:322609/1‑148 (MQ=255)
                                                                                                            ggtCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCaa                                   <  1:485716/149‑1 (MQ=255)
                                                                                                                    aCATAAACACTGACGTCGTAAACAACCGGTATAGCGGGGTTACCTTCACTGTACCGGCAGAGATCAATATTAGATGAATTCTGCCCCATCATTAAGGTATGGTCATAATAAACGAGTTAAGCAATAAAAGCAATATAGCAATAcagcag                           >  2:532432/1‑149 (MQ=255)
                                                                                                                       tAAACACTGACGTCGTAAACACCCGGTGTAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAgcgc                        <  2:690716/149‑1 (MQ=255)
                                                                                                                       tAAACACTGACGTCGTAAACACCCGGTGTAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAgcgc                        <  2:720525/149‑1 (MQ=255)
                                                                                                                            acTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACAGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGAATCGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGAAATACATCAGCGCGCaaa                   >  2:235437/1‑149 (MQ=255)
                                                                                                                                  gtcgtATCCACCCGGTATAGCGCGGATACCTTCACTGTACCGGCTGCGATGAATATTAGTGGCTTGCTGGCCCTGCTGGAAGGGTTGTTCATCTTCTAAGAGTTCAGCACTGTATGCTGTATTGCATTAGAGCAGCGCGATAAACGTgg             <  2:791610/149‑1 (MQ=255)
                                                                                                                                  gtcgtAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGTTCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTgg             >  1:284435/1‑149 (MQ=255)
                                                                                                                                  gtcgtAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTgg             <  1:37086/149‑1 (MQ=255)
                                                                                                                                   tcgtAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGc            >  1:977657/1‑149 (MQ=255)
                                                                                                                                   tcgtAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGc            <  2:977657/149‑1 (MQ=255)
                                                                                                                                       aaaCACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCCGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATa        >  1:535095/1‑149 (MQ=255)
                                                                                                                                          cacCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCCTCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACgg     <  2:168799/149‑1 (MQ=255)
                                                                                                                                          cacCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACgg     <  1:349413/149‑1 (MQ=255)
                                                                                                                                            cccGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGt   <  1:119739/149‑1 (MQ=255)
                                                                                                                                            cccGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGt   >  1:5277/1‑149 (MQ=255)
                                                                                                                                             ccGGTATAGCGGTGGTACCTTCACTGTACCGGCAGAGATCAAAATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGACAGCAGTATTGCACTACAGCAGCGCGCAAAACGTGGCGATACGGGTc  <  2:989726/149‑1 (MQ=255)
                                                                                                                                             ccGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATAAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGAATTGCAATACAGCAGCGCGCAAACCATGTCGATACGTGTc  >  1:342864/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CGGGAGAATTAATATGAAACTGCAATAAATTCTTTAATGTGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTC  >  NC_000913/155102‑155391

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: ‑