Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,058,696 A→G 5.7% T205A (ACG→GCG)  torC → trimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,058,6960AG 5.7% 196.7 / 2.7 71T205A (ACG→GCG) torCtrimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base A (36/30);  new base G (2/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

GCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAG  >  NC_000913/1058552‑1058841
                                                                                                                                                |                                                                                                                                                 
gcaCGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGcg                                                                                                                                               >  2:1015812/1‑148 (MQ=255)
    gTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTAc                                                                                                                                           <  1:106080/149‑1 (MQ=255)
        gATGAATGTGGCAGCGAAAGATAATCACTCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGATCGATGAGCTGCGCGCCAGTCCTAATGACAGTGGTGACACGCTGTACTCTa                                                                                                                                       >  2:1169212/1‑149 (MQ=255)
        gATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTa                                                                                                                                       <  2:879022/149‑1 (MQ=255)
          tGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTAtt                                                                                                                                     <  1:31756/149‑1 (MQ=255)
             tggaggCAGCGAATTATACTCAATCCTGCATCGTCTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGat                                                                                                                                  <  2:248674/145‑1 (MQ=255)
              ggTGTCAGCGAAAGATAATCAATCCTGCATCGACTGGCATAAAGGTATTGCCCACCAGGTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTGCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGGGACACGCTGTACTCTATTGata                                                                                                                                 <  2:893464/149‑1 (MQ=255)
               gTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGatat                                                                                                                                >  2:17726/1‑149 (MQ=255)
                tGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATAtt                                                                                                                               <  1:478619/149‑1 (MQ=255)
                  gCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTaa                                                                                                                             >  1:93996/1‑149 (MQ=255)
                    aGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGc                                                                                                                           >  1:479262/1‑149 (MQ=255)
                    aGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGc                                                                                                                           <  1:525792/149‑1 (MQ=255)
                    aGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGAGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACATTGGTGACACGCTGTACTCTATTGATATTAAGc                                                                                                                           >  2:341009/1‑149 (MQ=255)
                      cGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCa                                                                                                                         >  1:224314/1‑149 (MQ=255)
                          aGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACGCGCTGTACTCTATTGATATTAAGCCAAttt                                                                                                                     >  1:624898/1‑149 (MQ=255)
                          aGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACGCGCTGTACTCTATTGATATTAAGCCAAttt                                                                                                                     >  1:232657/1‑149 (MQ=255)
                               atcaatcCTGCATCGACTGTCATAAAGGTATTGCCAACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTTCTAATGACAGTCGTGACACGCTGTAATCTATTGATATTAAGCCAATTTATgcg                                                                                                                >  2:729392/1‑149 (MQ=255)
                                tcaatcCTGCATCGACTGTCATAAAGGTATTGCCCACTAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATgcgg                                                                                                               <  1:679943/149‑1 (MQ=255)
                                   atcCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACCCGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGa                                                                                                            >  1:1021265/1‑149 (MQ=255)
                                       tGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAgg                                                                                                        <  1:58042/149‑1 (MQ=255)
                                           tCGACTGTCATAAAGGTATTGCCCACCAGTTCCCGGATATGAGTAGCGGCTTCCGTAAGCAGGTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGAt                                                                                                    <  2:107840/149‑1 (MQ=255)
                                            cGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATa                                                                                                   <  1:17726/149‑1 (MQ=255)
                                                    aTAAAGGTAGTGCACACAAGTTACAGGATATGATTAGAGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATGTTAAGCCAATTTATGCGGCGACAGGCAATAAAGAAGca                                                                                           >  1:985672/1‑148 (MQ=255)
                                                     tAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCct                                                                                          >  1:571988/1‑149 (MQ=255)
                                                        aGGTATTGCCCACCAGTTACCGGCTATGAGTCGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGGGGTGACACGCTGGACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTg                                                                                       <  2:290438/149‑1 (MQ=255)
                                                        aGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTg                                                                                       <  1:1075154/149‑1 (MQ=255)
                                                             ttGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCaa                                                                                                                        >  2:725966/1‑111 (MQ=255)
                                                             ttGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCaa                                                                                                                        <  1:965139/111‑1 (MQ=255)
                                                             ttGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCaa                                                                                                                        >  2:965139/1‑111 (MQ=255)
                                                             ttGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCaa                                                                                                                        <  1:725966/111‑1 (MQ=255)
                                                               gCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACGCGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTctct                                                                                <  2:232657/149‑1 (MQ=255)
                                                               gCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACGCGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTctct                                                                                <  2:624898/149‑1 (MQ=255)
                                                                  caccaGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTctgct                                                                             >  1:450030/1‑149 (MQ=255)
                                                                  caccaGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTATctgct                                                                             >  2:293042/1‑149 (MQ=255)
                                                                        ttACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGctgcctgc                                                                       >  2:851655/1‑149 (MQ=255)
                                                                        ttACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATAGATATTAAGCCAATTTATGCCGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGctgcctgc                                                                       >  2:404844/1‑149 (MQ=255)
                                                                           ccGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTc                                                                    >  2:393314/1‑149 (MQ=255)
                                                                              gATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGa                                                                 >  2:502369/1‑149 (MQ=255)
                                                                                  tGAGTAGCGGCGTCCGGAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTg                                                             <  1:759304/149‑1 (MQ=255)
                                                                                   gAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGa                                                            >  2:23459/1‑149 (MQ=255)
                                                                                   gAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGa                                                            >  2:323792/1‑149 (MQ=255)
                                                                                        gCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATTAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTc                                                       >  2:725498/1‑149 (MQ=255)
                                                                                         cGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGATAGTcc                                                      <  1:201199/149‑1 (MQ=255)
                                                                                           gCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCtt                                                    >  1:1102453/1‑149 (MQ=255)
                                                                                               ccGTAAGCAGTTCGATGAGCTGCGCGCCAGTTCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAAc                                                <  1:237940/149‑1 (MQ=255)
                                                                                                 gTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGt                                              >  2:2203/1‑149 (MQ=255)
                                                                                                  tAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTTGAGACACTCTGTACTCTATTTATATTAAGCCAATTTATGCTGCGAAAGGCGATAAATAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTg                                             <  1:1091536/149‑1 (MQ=255)
                                                                                                        gTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTTTGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCg                                       >  2:100483/1‑149 (MQ=255)
                                                                                                         ttCGATGCGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGATGCCTGCTTCGGAAGTGAAAGTCCTTAAACTTGACGGCGa                                      >  1:572866/1‑149 (MQ=255)
                                                                                                         ttCGATGAGCTGCGCGCCAGTTCTAATGACAGTGGTGACACGATGTACTCTATTGATATTAAGTCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGa                                      <  2:1122338/149‑1 (MQ=255)
                                                                                                         ttCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGa                                      <  2:904750/149‑1 (MQ=255)
                                                                                                             aTGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTgg                                  <  1:868780/149‑1 (MQ=255)
                                                                                                             aTGAGCTGCGCGCCAGAGCTACTGACAGTGGTGACACGCTGGACTCTATTGCTATTAAGCCAATTTATGCCGCGAAAGGCGATAAAGAAGCCTCTGGGTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTgg                                  <  1:253612/149‑1 (MQ=255)
                                                                                                               gAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCt                                <  2:158850/149‑1 (MQ=255)
                                                                                                                  cTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGATGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCa                             >  1:358935/1‑149 (MQ=255)
                                                                                                                  cTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCa                             >  2:119114/1‑149 (MQ=255)
                                                                                                                   tGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAATGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTTCTTCCGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCaa                            >  2:707632/1‑149 (MQ=255)
                                                                                                                    gcgcgcCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTTCTGACTTCTTCTGAATGGAAAGTCCGGAAACAGGACTGAGACTGGCTgaata                           >  2:1187191/1‑145 (MQ=255)
                                                                                                                    gcgcgcCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACgg                                         >  2:1165315/1‑135 (MQ=255)
                                                                                                                    gcgcgcCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACgg                                         <  1:1165315/135‑1 (MQ=255)
                                                                                                                         ccAGTTCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCGTCTGGTTCTCTGCTGCCTGCTTCGGATGTGAAAGTCCTTAATCGTGACGGCGACTGGCTGCAAATTGaa                      <  2:278266/149‑1 (MQ=255)
                                                                                                                         ccAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGaa                      >  2:455575/1‑149 (MQ=255)
                                                                                                                              gCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTAc                 <  2:1066654/149‑1 (MQ=255)
                                                                                                                              gCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAAATTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTAc                 <  2:620409/149‑1 (MQ=255)
                                                                                                                                 aaTGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCgg              >  1:750078/1‑149 (MQ=255)
                                                                                                                                     aCAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGAGACGGCGACTGGCTGCTAATTGAAATTACCGGCTgg          >  1:845301/1‑149 (MQ=255)
                                                                                                                                       cgtggtgACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCCGCAAATTGAAATTACCGGCTGGAc        <  2:1127388/148‑1 (MQ=255)
                                                                                                                                         tggtgACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACgg      >  1:103330/1‑149 (MQ=255)
                                                                                                                                           gtgACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCATCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGaa    >  1:955899/1‑149 (MQ=255)
                                                                                                                                            tgACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGTTTCTCTGCTGCCTGCTTCGGAAGTGAAAGGCCTTAAACGTTAAGGCGAATGCATGCAAATTCAAATTACCGCCTGTACGGaaa   >  2:1162999/1‑149 (MQ=255)
                                                                                                                                             gccacGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAATCGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGaaac  >  1:388411/3‑148 (MQ=255)
                                                                                                                                             gACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAg  <  1:5118/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAG  >  NC_000913/1058552‑1058841

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: