Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 155,285 T→C 6.1% I48V (ATT→GTT)  htrE ← putative outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913155,2850TC 6.1% 204.6 / 3.8 65I48V (ATT→GTT) htrEputative outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base C (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

TGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATG  >  NC_000913/155141‑155432
                                                                                                                                                |                                                                                                                                                   
tgATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTa                                                                                                                                                 <  3:532029/149‑1 (MQ=255)
 gATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAg                                                                                                                                                >  4:477020/1‑149 (MQ=255)
      aaGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCa                                                                                                                                           >  2:510294/1‑149 (MQ=255)
                 ttCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCAt                                                                                                                                >  4:170130/1‑149 (MQ=255)
                  tCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGTGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATc                                                                                                                               <  1:510294/149‑1 (MQ=255)
                  tCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATc                                                                                                                               <  3:534784/149‑1 (MQ=255)
                   cTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCa                                                                                                                              <  2:280884/149‑1 (MQ=255)
                   cTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCa                                                                                                                              <  4:402384/149‑1 (MQ=255)
                       tttCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGaa                                                                                                                          <  4:214894/149‑1 (MQ=255)
                           cTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGc                                                                                                                                     <  1:108059/134‑1 (MQ=255)
                           cTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGc                                                                                                                                     >  2:108059/1‑134 (MQ=255)
                            ttCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTAt                                                                                                                     >  1:43132/1‑149 (MQ=255)
                             tCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATg                                                                                                                    <  3:5425/149‑1 (MQ=255)
                               aaTTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGt                                                                                                                  >  1:426894/1‑149 (MQ=255)
                                      aCAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTc                                                                                                           <  4:325402/149‑1 (MQ=255)
                                      aCAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTc                                                                                                           <  3:486478/149‑1 (MQ=255)
                                        aaaTGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCaa                                                                                                         <  3:265671/149‑1 (MQ=255)
                                         aaTGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAAc                                                                                                        <  1:384236/149‑1 (MQ=255)
                                           tGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGa                                                                                                      <  4:137776/149‑1 (MQ=255)
                                           tGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATGCTGCCCCATCAGGAAGGTATGGTCATATTCAACGa                                                                                                      <  2:316975/149‑1 (MQ=255)
                                             tAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGt                                                                                                    <  1:556250/149‑1 (MQ=255)
                                              aaTACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCAGTGTAGCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGtt                                                                                                   >  3:536388/1‑149 (MQ=255)
                                                   tttGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCa                                                                                              >  4:221376/1‑149 (MQ=255)
                                                            tGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCa                                                                                     <  1:493749/149‑1 (MQ=255)
                                                                 ggttggtCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTAtt                                                                                <  2:129039/149‑1 (MQ=255)
                                                                   ttggtCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGc                                                                              >  2:183119/1‑149 (MQ=255)
                                                                             aCATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATAcagcag                                                                    >  3:269906/1‑149 (MQ=255)
                                                                             aCATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATAcagcag                                                                    >  1:41726/1‑149 (MQ=255)
                                                                               aTAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAgcg                                                                  >  4:379590/1‑149 (MQ=255)
                                                                                 aaaCACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGc                                                                                                                                     >  3:270058/1‑80 (MQ=255)
                                                                                 aaaCACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGc                                                                                                                                     <  4:270058/80‑1 (MQ=255)
                                                                                 aaaCACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAgcgcg                                                                >  3:340714/1‑149 (MQ=255)
                                                                                     acTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCaaa                                                            >  1:407306/1‑149 (MQ=255)
                                                                                     acTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCaaa                                                            >  1:483105/1‑149 (MQ=255)
                                                                                         aCGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGt                                                        >  1:382594/1‑149 (MQ=255)
                                                                                         aCGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGt                                                        <  3:457873/149‑1 (MQ=255)
                                                                                           gtcgtAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTgg                                                      <  1:192173/149‑1 (MQ=255)
                                                                                            tcgtAAACACCCGGTATAACGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTAAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGc                                                     >  2:247470/1‑149 (MQ=255)
                                                                                              gtAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGa                                                   <  4:536388/149‑1 (MQ=255)
                                                                                              gtAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGa                                                   >  2:1308/1‑149 (MQ=255)
                                                                                                aaaCACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATa                                                 <  1:381404/149‑1 (MQ=255)
                                                                                                 aaCACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATAc                                                <  4:32653/149‑1 (MQ=255)
                                                                                                 aaCACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATAc                                                >  2:527131/1‑149 (MQ=255)
                                                                                                    acCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACggg                                             >  1:22728/1‑149 (MQ=255)
                                                                                                    acCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAACATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACggg                                             >  2:251200/1‑149 (MQ=255)
                                                                                                    acCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAACATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACggg                                             >  4:10822/1‑149 (MQ=255)
                                                                                                             aGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGAAGCGCGCAAAACGTGGCGATACGGGTCAgatgat                                    <  2:440320/149‑1 (MQ=255)
                                                                                                                 gggTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATa                                >  1:199532/1‑149 (MQ=255)
                                                                                                                  ggTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATaa                               <  3:255359/149‑1 (MQ=255)
                                                                                                                  ggTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATaa                               >  4:130119/1‑149 (MQ=255)
                                                                                                                    ttACCTTCACTGTACCGGCTGAGATCAACATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAAtt                             <  1:251200/149‑1 (MQ=255)
                                                                                                                    ttACCTTCACTGTACCGGCTGAGATCAACATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAAtt                             <  3:10822/149‑1 (MQ=255)
                                                                                                                      aCCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAAtttt                           <  1:198748/149‑1 (MQ=255)
                                                                                                                        cTTCACTGTAGCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTa                         <  4:408769/149‑1 (MQ=255)
                                                                                                                        cTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTa                         >  3:315898/1‑149 (MQ=255)
                                                                                                                        cTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTa                         <  2:41726/149‑1 (MQ=255)
                                                                                                                         ttCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAg                        >  2:506794/1‑149 (MQ=255)
                                                                                                                                tACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTc                 <  3:170130/149‑1 (MQ=255)
                                                                                                                                  ccGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCta               <  3:221376/149‑1 (MQ=255)
                                                                                                                                     gCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATATGGGTCAGATGATGATAATTTTTAGTATATTCTATAg            <  1:183119/149‑1 (MQ=255)
                                                                                                                                      cTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGt           >  2:19300/1‑149 (MQ=255)
                                                                                                                                           aTCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTTAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGc      >  2:129560/1‑149 (MQ=255)
                                                                                                                                           aTCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGc      >  2:506883/1‑149 (MQ=255)
                                                                                                                                           aTCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGc      >  1:458119/1‑149 (MQ=255)
                                                                                                                                            tCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCt     >  1:368944/1‑149 (MQ=255)
                                                                                                                                               atatTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACgg                                              >  2:500160/1‑105 (MQ=255)
                                                                                                                                               atatTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACgg                                              <  1:500160/105‑1 (MQ=255)
                                                                                                                                               atatTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATg  >  1:494313/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TGATACAAGCCTGGGCGTTCTTTTTTCCTTCAATTGCGACAAATGTAATACTTTGGTTAATGATTGGTTGGTCGTTTACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATG  >  NC_000913/155141‑155432

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑