Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,060,648 A→G 6.0% I465V (ATC→GTC)  torA → trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,060,6480AG 6.0% 199.6 / 3.2 66I465V (ATC→GTC) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (28/34);  new base G (2/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

ACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGCCATAGATAACCAGTG  >  NC_000913/1060501‑1060796
                                                                                                                                                   |                                                                                                                                                    
aCGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGAt                                                                                                                                                     <  4:317777/149‑1 (MQ=255)
     ggggCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACt                                                                                                                                                <  7:113087/149‑1 (MQ=255)
       ggCGTAAAGGCGTTATTCTGAGTGGCGTCTCCGGCGCTACGTCGAATCCGCCGGGTCACGACAACAGTGACTCCAAAGGCTCCAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGACCCGGGGAAAGGGATCAACTgg                                                                                                                                              <  4:29447/149‑1 (MQ=255)
         cGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGaa                                                                                                                                            >  7:191595/1‑149 (MQ=255)
           tAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACg                                                                                                                                          <  3:298594/149‑1 (MQ=255)
                   ttATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCg                                                                                                                                  >  1:297784/1‑149 (MQ=255)
                   ttATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCACGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCg                                                                                                                                  <  1:124178/149‑1 (MQ=255)
                     aTTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCgg                                                                                                                                 >  3:281925/1‑148 (MQ=255)
                     aTTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCgg                                                                                                                                 >  7:340599/1‑148 (MQ=255)
                      ttCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGt                                                                                                                                <  4:281925/148‑1 (MQ=255)
                      ttCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGt                                                                                                                                <  8:340599/148‑1 (MQ=255)
                      ttCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGt                                                                                                                                <  3:61160/148‑1 (MQ=255)
                      ttCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTa                                                                                                                               <  8:191595/149‑1 (MQ=255)
                                tttCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTgccgc                                                                                                                     >  5:290080/1‑149 (MQ=255)
                                  tctcCGGCTCTACGTCGATTCCGCCTGTTCACGCCAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGccgccg                                                                                                                   >  2:205510/1‑149 (MQ=255)
                                          tctACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTTCAGCAGCAGTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAAt                                                                                                           >  3:303154/1‑149 (MQ=255)
                                             aCGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAAtgtg                                                                                                        >  5:77226/1‑149 (MQ=255)
                                              cGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAAtgt                                                                                                         <  8:203393/147‑1 (MQ=255)
                                                   aTTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTAtttt                                                                                                  <  2:269268/149‑1 (MQ=255)
                                                        gCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCg                                                                                             <  4:169828/149‑1 (MQ=255)
                                                          cttttCACGACCACATGGACTACAAAGGCTAAAGCAGCACTATTCCGATTTCCCGTTTTATCGAGTCGATCCTCGAACCGGGGAAATTGATCAACTGGTACGGTAAATCGGTAAAACTGCCGCCGCGGAAAATGTGTATTTTTGCCTGa                                                                                           <  6:111229/146‑1 (MQ=255)
                                                            gTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAAc                                                                                         >  1:146034/1‑149 (MQ=255)
                                                             ttCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGaa                                                                                          <  6:77226/147‑1 (MQ=255)
                                                                aCGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGACGCCGCTGAAAATGTGTATTTTTGCCGGAACTAAc                                                                                     <  2:305003/149‑1 (MQ=255)
                                                                aCGACAACAGTGACTACAAAGCATACAGCAGCACTATTCCGATTGCCCGTTTTATCTATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACTGTAAATCGGTAACACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAAc                                                                                     <  6:335607/149‑1 (MQ=255)
                                                                  gACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGGTCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAAc                                                                                     >  2:153143/1‑147 (MQ=255)
                                                                  gACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGGTCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAAc                                                                                     <  1:153143/147‑1 (MQ=255)
                                                                     aacaGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCAtt                                                                                <  5:233241/149‑1 (MQ=255)
                                                                           gACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCg                                                                          >  5:205394/1‑149 (MQ=255)
                                                                             cTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGcc                                                                        <  7:97504/149‑1 (MQ=255)
                                                                              tACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTaa                                                                                      <  1:344726/134‑1 (MQ=255)
                                                                              tACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTaa                                                                                      >  2:344726/1‑134 (MQ=255)
                                                                               aCAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCAt                                                                      >  7:186270/1‑149 (MQ=255)
                                                                                cAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTAttttt                                                                                                 >  4:173608/1‑121 (MQ=255)
                                                                                cAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTAttttt                                                                                                 <  3:173608/121‑1 (MQ=255)
                                                                                cAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCa                                                                       <  8:69852/147‑1 (MQ=255)
                                                                                cAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCAt                                                                      >  2:213296/1‑148 (MQ=255)
                                                                                  aaGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATcag                                                                   <  8:25896/149‑1 (MQ=255)
                                                                                  aaGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATcag                                                                   >  7:168981/1‑149 (MQ=255)
                                                                                        aCAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAATGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGAt                                                              >  8:337465/1‑148 (MQ=255)
                                                                                            cagcaCTATTCCGATTGTCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGGTCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCg                                                                          >  7:8383/1‑132 (MQ=255)
                                                                                            cagcaCTATTCCGATTGTCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGGTCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCg                                                                          <  8:8383/132‑1 (MQ=255)
                                                                                            cagcaCTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGACCCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAAcc                                                         >  2:205177/1‑149 (MQ=255)
                                                                                                    ttCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCattatt                                                 <  7:39094/149‑1 (MQ=255)
                                                                                                      ccGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTg                                                >  7:285517/1‑148 (MQ=255)
                                                                                                      ccGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTg                                                >  7:280562/1‑148 (MQ=255)
                                                                                                       cGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAATGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGaa                                              <  7:337465/149‑1 (MQ=255)
                                                                                                         aTTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCACCTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCAGCAGCAGATCAACCGCATTATTGAAgg                                            >  6:116440/1‑149 (MQ=255)
                                                                                                             cccGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTg                                        >  2:96339/1‑149 (MQ=255)
                                                                                                              ccGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCTCCGCTGAAAATGTGTATTTTTCCCGGAACTACCCCATTGCATCGCCATCAGCAGATCAACCGCATAATAGATGGCTTgc                                       >  2:219612/1‑149 (MQ=255)
                                                                                                              ccGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTgc                                       <  3:80503/149‑1 (MQ=255)
                                                                                                              ccGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTgc                                       <  6:290080/149‑1 (MQ=255)
                                                                                                                    tATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAg                                  <  2:146034/148‑1 (MQ=255)
                                                                                                                            cGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAAc                          <  4:303154/148‑1 (MQ=255)
                                                                                                                             gATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACgg                        >  1:166888/1‑149 (MQ=255)
                                                                                                                             gATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACgg                        <  8:288313/149‑1 (MQ=255)
                                                                                                                               tCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGtt                      <  8:285517/149‑1 (MQ=255)
                                                                                                                               tCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGtt                      <  8:280562/149‑1 (MQ=255)
                                                                                                                                 cTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTa                     <  1:205510/148‑1 (MQ=255)
                                                                                                                                   cGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGCAACGCTTATCg                  >  2:168504/1‑149 (MQ=255)
                                                                                                                                    gAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGc                 >  8:210926/1‑149 (MQ=255)
                                                                                                                                     aaCCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGc                 >  6:161852/1‑148 (MQ=255)
                                                                                                                                          gggAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGCCa               >  7:62607/1‑145 (MQ=255)
                                                                                                                                           ggAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGCCATAGAt          <  8:168981/149‑1 (MQ=255)
                                                                                                                                           ggAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTAATCGCCATAGAt          <  1:205177/149‑1 (MQ=255)
                                                                                                                                                  gATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGATACGGTTATCGCCATAGATAACCAGt   >  5:249023/1‑149 (MQ=255)
                                                                                                                                                   aTCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGACGGCTTGCGCAAGCTGGAAACGGTTATCGCCATAGATAACCAGTg  >  1:217250/1‑149 (MQ=255)
                                                                                                                                                   aTCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGCCATAGATAACCAGt   <  1:213296/148‑1 (MQ=255)
                                                                                                                                                   aTCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGCCATAGATAACCAGTg  <  7:195020/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGCCATAGATAACCAGTG  >  NC_000913/1060501‑1060796

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