Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 119,295 A→T 5.7% T5T (ACA→ACT ampE → ampicillin resistance inner membrane protein; putative signaling protein in beta‑lactamase regulation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913119,2950AT 5.7% 243.7 / 3.8 80T5T (ACA→ACTampEampicillin resistance inner membrane protein; putative signaling protein in beta‑lactamase regulation
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base T (2/3);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

TGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGT  >  NC_000913/119151‑119438
                                                                                                                                                 |                                                                                                                                               
tgattgatTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAAcc                                                                                                                                              >  8:226792/1‑148 (MQ=255)
 gattgattGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCt                                                                                                                                             >  7:2774/1‑148 (MQ=255)
 gattgattGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATACCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCt                                                                                                                                             >  3:230375/1‑148 (MQ=255)
    tgattgCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTAc                                                                                                                                          >  8:119917/1‑148 (MQ=255)
    tgattgCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTAc                                                                                                                                          >  2:332294/1‑148 (MQ=255)
            aTCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGtt                                                                                                                                   >  2:375290/1‑147 (MQ=255)
               cGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAAt                                                                                                                                <  8:10395/147‑1 (MQ=255)
                  atatCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTc                                                                                                                            >  7:112358/1‑148 (MQ=255)
                     acGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGTATCGGAAAACCGATCCCGGTCCTGCATTTGACTGGGCACGGTCTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAg                                                                                                                         <  6:204846/147‑1 (MQ=255)
                      cGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAgc                                                                                                                        >  3:107929/1‑148 (MQ=255)
                       gCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAgcg                                                                                                                       >  1:125462/1‑148 (MQ=255)
                        cTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCCATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAgcgc                                                                                                                      >  7:81004/1‑148 (MQ=255)
                          aaaaaCCTGAGGGTGCATTGTGATATTGCGGCGGATCGGATAAGCGATCCCGGTCCTGCTTTTGATTGGGCACGTTTTCGTGTGCTGGTCAGCAAGGAGACATCATGACGCTATTTACAACCtt                                                                                                                                            <  7:318731/124‑1 (MQ=255)
                              aCATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGtt                                                                                                                 >  3:357475/1‑147 (MQ=255)
                               cATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTa                                                                                                               <  4:357475/148‑1 (MQ=255)
                                 tGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAg                                                                                                             >  3:395514/1‑148 (MQ=255)
                                   aCGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGt                                                                                                            <  1:421321/147‑1 (MQ=255)
                                   aCGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGt                                                                                                            <  6:164992/147‑1 (MQ=255)
                                         cATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGa                                                                                                     <  7:226792/148‑1 (MQ=255)
                                            tgtgATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTaa                                                                                                              >  8:55176/1‑136 (MQ=255)
                                            tgtgATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTaa                                                                                                              <  7:55176/136‑1 (MQ=255)
                                            tgtgATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCa                                                                                                  >  1:89544/1‑148 (MQ=255)
                                                  atTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACTACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCg                                                                                                      >  8:229565/1‑138 (MQ=255)
                                                  atTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACTACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCg                                                                                                      <  7:229565/138‑1 (MQ=255)
                                                  atTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACTACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCg                                                                                                      <  7:235621/138‑1 (MQ=255)
                                                  atTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACTACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCg                                                                                                      >  8:235621/1‑138 (MQ=255)
                                                  atTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTgg                                                                                              >  1:231728/1‑146 (MQ=255)
                                                   ttGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAg                                                                                           >  4:70843/1‑148 (MQ=255)
                                                      cgcCGGATCGGAAAACCGATCCCTGTCCTGCATTTGATTGTGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCtt                                                                                        <  8:18032/148‑1 (MQ=255)
                                                        ccGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGa                                                                                      <  2:89544/148‑1 (MQ=255)
                                                          ggATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACATCATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGat                                                                                     >  3:393028/1‑147 (MQ=255)
                                                          ggATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGat                                                                                     >  6:107069/1‑147 (MQ=255)
                                                             tCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGatcatc                                                                                 >  6:259862/1‑148 (MQ=255)
                                                              cGGAAAACCGATCCCGGTCCTGCATTTGATTAGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTTATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGatcat                                                                                  <  4:26604/147‑1 (MQ=255)
                                                              cGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGc                                                                                                   <  7:162239/129‑1 (MQ=255)
                                                              cGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGc                                                                                                   >  8:162239/1‑129 (MQ=255)
                                                               ggAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGt                                                                               <  2:95759/148‑1 (MQ=255)
                                                                 aaaaCCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCt                                                                             <  3:152437/148‑1 (MQ=255)
                                                                       gATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGc                                                                       >  7:408255/1‑148 (MQ=255)
                                                                          cccGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGtt                                                                    >  4:380095/1‑148 (MQ=255)
                                                                           ccGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGtt                                                                    <  5:259862/147‑1 (MQ=255)
                                                                             ggTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCtt                                                                 <  8:112358/148‑1 (MQ=255)
                                                                                  tGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGt                                                                                                            >  7:346872/1‑100 (MQ=255)
                                                                                  tGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGt                                                                                                            <  8:346872/100‑1 (MQ=255)
                                                                                   gCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTcggcg                                                           >  3:293956/1‑148 (MQ=255)
                                                                                    cATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTcggcg                                                           <  1:375290/147‑1 (MQ=255)
                                                                                      tttGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTAGAACCTTACTGGTGTTAATTTTCGTGGGGCTGTTTTAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGt                                                        <  1:135266/148‑1 (MQ=255)
                                                                                       ttGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGgtg                                                                                                                                     >  5:352498/1‑70 (MQ=255)
                                                                                       ttGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGgtg                                                                                                                                     <  6:352498/70‑1 (MQ=255)
                                                                                          aTTGGGCAAGTTTTCGTGTGCTGGTCAGCAAGGAGACAAGATGACGCTATTTACATCCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGGTCTTTCGGCGTGTGaa                                                     <  3:389157/147‑1 (MQ=255)
                                                                                                 cACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACAtttt                                              >  7:131334/1‑148 (MQ=255)
                                                                                                     gatgtgtgtGGTGGTGAGCAAGGAGACTACATGTGGATATTTACTAGATTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCAATGGCAGCTTGATCATCGTCTTGAAGAGTTATTTCGGCGGGTGAAACATTTTtct                                           <  4:371751/142‑1 (MQ=255)
                                                                                                      tttCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTtctct                                         <  1:332294/148‑1 (MQ=255)
                                                                                                      tttCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTtctct                                         <  7:119917/148‑1 (MQ=255)
                                                                                                        tCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTtctctc                                        <  8:2774/147‑1 (MQ=255)
                                                                                                         cGTGTGCTGGTCAGCAAGGAGTCTACATGACGCTTTTTAGAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTtctctgg                                       <  8:81004/147‑3 (MQ=255)
                                                                                                         cGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCg                                       <  8:408255/147‑1 (MQ=255)
                                                                                                            gtgCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGgcg                                   >  6:284360/1‑148 (MQ=255)
                                                                                                            gtgCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCATGTTTAAGTTGGGCGAGCACTGGCATCTTGATCATCGTCATGAAGAGTTCTTTCGGCGGGTGAAACATTTTtctctagatct                                   >  3:371751/1‑142 (MQ=255)
                                                                                                                 ggTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCa                                 >  6:45942/1‑145 (MQ=255)
                                                                                                                 ggTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCa                                 <  5:45942/145‑1 (MQ=255)
                                                                                                                 ggTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGt                              >  3:17455/1‑148 (MQ=255)
                                                                                                                      gCAAGGAGACAACATGACGCTATTTACAACCTTACTGGAGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGTCGGGTGAAACATTTTTCTCTCGGGCGCACGTAAgg                          >  6:51766/1‑147 (MQ=255)
                                                                                                                        aaGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCAt                       >  8:110839/1‑148 (MQ=255)
                                                                                                                           gagaCAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTATAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTTATCATCGTCTTGAAGCGTTATTTCGTCTGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGAATGa                     >  5:53008/1‑147 (MQ=255)
                                                                                                                            agaAAACATCACGCTATTAACAACCTAACTGGTGATAATTTACGAGCGCCTGTTTAAGTACGGCGAGCACTGGCAGCTTGATCATCCTCTTCAAGCGATCTTTCGGCGGGTGAAACAATTTACTCTCGGGCGCACGATACGCATGAc                    >  4:80505/1‑147 (MQ=255)
                                                                                                                               caacaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGt                                                                     >  6:113185/1‑95 (MQ=255)
                                                                                                                               caacaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGt                                                                     <  5:113185/95‑1 (MQ=255)
                                                                                                                                 acaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCattat              <  2:66237/148‑1 (MQ=255)
                                                                                                                                 acaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCattat              >  8:19648/1‑148 (MQ=255)
                                                                                                                                 acaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGAACATCTTCTTGAAGCTTTATTTCGGCGGGTGAAACAATTTTATATCTGGCGCACGTTAGGCATGACAattat              >  3:84134/1‑148 (MQ=255)
                                                                                                                                      aCGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCg          <  4:393028/147‑1 (MQ=255)
                                                                                                                                      aCGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCg          <  2:231728/147‑1 (MQ=255)
                                                                                                                                      aCGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCg          <  3:380095/147‑1 (MQ=255)
                                                                                                                                      aCGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCg          <  4:17455/147‑1 (MQ=255)
                                                                                                                                       cGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGAt        <  3:70843/148‑1 (MQ=255)
                                                                                                                                           aTTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGc    >  1:118066/1‑148 (MQ=255)
                                                                                                                                             ttACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGt  >  4:46823/1‑148 (MQ=255)
                                                                                                                                              ttcAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGc                         >  3:10823/3‑124 (MQ=255)
                                                                                                                                              ttcAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGc                         <  4:10823/122‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGT  >  NC_000913/119151‑119438

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