Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,957,069 A→C 5.5% F355C (TTC→TGC)  torY ← TMAO reductase III (TorYZ), cytochrome c‑type subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,957,0690AC 5.5% 283.9 / 6.3 91F355C (TTC→TGC) torYTMAO reductase III (TorYZ), cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base A (41/44);  new base C (2/3);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTT  >  NC_000913/1956926‑1957215
                                                                                                                                               |                                                                                                                                                   
ggAGCGGTGCTGCCTATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATa                                                                                                                                                 <  2:82421/148‑1 (MQ=255)
   gCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGCAATACcg                                                                                                                                              >  2:55003/1‑148 (MQ=255)
   gCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGCAATACcg                                                                                                                                              >  2:61751/1‑148 (MQ=255)
   gCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACcg                                                                                                                                              <  6:24114/148‑1 (MQ=255)
   gCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACcg                                                                                                                                              >  7:322420/1‑148 (MQ=255)
        gctgcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGtt                                                                                                                                          >  3:334159/1‑147 (MQ=255)
           gcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATa                                                                                                                                                 <  8:235930/137‑1 (MQ=255)
           gcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATa                                                                                                                                                 >  7:235930/1‑137 (MQ=255)
            ccGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATGACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTaaa                                                                                                                                       <  2:215116/146‑1 (MQ=255)
             cGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATg                                                                                                                                    <  6:306845/148‑1 (MQ=255)
              gATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGt                                                                                                                                   <  6:251367/148‑1 (MQ=255)
               aTGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAAGGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTc                                                                                                                                  <  5:206130/148‑1 (MQ=255)
                 gTCACCACCAGCGCTCCAGCGGCAATAACGCTGTGTTTAATAAATTAACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTc                                                                                                                                  >  5:440044/1‑146 (MQ=255)
                   caccacCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGt                                                                                                                              >  7:389757/1‑148 (MQ=255)
                   caccacCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGt                                                                                                                              >  7:112261/1‑148 (MQ=255)
                    accaccAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGGTGTACTGGAAATACCGCGTTAAAATGTCCAGt                                                                                                                              <  2:302335/147‑1 (MQ=255)
                     ccaccaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGt                                                                                                                              >  2:310390/1‑146 (MQ=255)
                     ccaccaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGt                                                                                                                              <  1:310390/146‑1 (MQ=255)
                      caccaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTaaa                                                                                                                                       >  3:212008/1‑136 (MQ=255)
                      caccaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTaaa                                                                                                                                       <  4:212008/136‑1 (MQ=255)
                        ccaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGtt                                                                                                                        >  8:123273/1‑148 (MQ=255)
                        ccaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCAGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGtt                                                                                                                        >  8:116573/1‑148 (MQ=255)
                         caGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGttt                                                                                                                       >  1:457088/1‑148 (MQ=255)
                         caGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGttt                                                                                                                       <  8:203108/148‑1 (MQ=255)
                            cgcTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCGGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCg                                                                                                                    >  8:127521/1‑148 (MQ=255)
                             gcTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGc                                                                                                                   <  5:288104/148‑1 (MQ=255)
                              cTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCGGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCt                                                                                                                  <  7:127521/148‑1 (MQ=255)
                                ccAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTAAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTTTTTCTCGGTAATATCTTTGGCGTTGTACTGGATATACACAGTTAAAATGTCCATTTCGTTTTCTCTCa                                                                                                                >  8:230443/1‑148 (MQ=255)
                                 cAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGCAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTc                                                                                                                 <  1:61751/146‑1 (MQ=255)
                                 cAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGCAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTc                                                                                                                 <  1:55003/146‑1 (MQ=255)
                                 cAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTc                                                                                                                 <  2:457088/146‑1 (MQ=255)
                                 cAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTc                                                                                                                 <  8:322420/146‑1 (MQ=255)
                                   gCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGc                                                                                                             >  1:195338/1‑148 (MQ=255)
                                         aTACCGCTGTGTTTAATATATTCACGTCTAGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTGCTCGGTAATATCTTTGGCGTTGTACTGGAAATAGCGCGTTAAAATGTCCAGTTCGTTTTCGCTCTTGCTGTTTc                                                                                                       <  6:20783/148‑1 (MQ=255)
                                         aTACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGGAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGGTCGTTTTCGCTCATGCTGGTTc                                                                                                       >  4:54507/1‑148 (MQ=255)
                                           aCCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTg                                                                                                           <  5:399479/142‑1 (MQ=255)
                                           aCCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTg                                                                                                           >  6:399479/1‑142 (MQ=255)
                                           aCCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGt                                                                                                     >  8:56406/1‑148 (MQ=255)
                                           aCCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATAACTGTTTATAGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTGAGTTTTCTCTCa                                                                                                                >  1:382443/1‑137 (MQ=255)
                                             cGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGc                                                                                                   <  2:195338/148‑1 (MQ=255)
                                                 gtgtTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGcccc                                                                                                >  8:306239/1‑147 (MQ=255)
                                                      taataaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGGTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTGTTCGCTCATGCTGGTTCGTGCCCCCATTcc                                                                                          >  2:113185/1‑148 (MQ=255)
                                                         taaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCTTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCttt                                                                                       >  4:360400/1‑148 (MQ=255)
                                                             ttCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCa                                                                                   <  5:280938/148‑1 (MQ=255)
                                                             ttCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCa                                                                                   <  3:15527/148‑1 (MQ=255)
                                                                 cGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATgg                                                                               <  2:428022/148‑1 (MQ=255)
                                                                          aaTGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCAc                                                                       <  8:389757/147‑1 (MQ=255)
                                                                           aTGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGt                                                                                                                                           <  4:332895/79‑1 (MQ=255)
                                                                           aTGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGt                                                                                                                                           >  3:332895/1‑79 (MQ=255)
                                                                              tCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCAtt                                                                  >  6:267478/1‑148 (MQ=255)
                                                                               cATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTg                                                                 >  2:47451/1‑148 (MQ=255)
                                                                                  aaCTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACg                                                              <  1:308087/148‑1 (MQ=255)
                                                                                   aCTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACgg                                                             >  6:416189/1‑148 (MQ=255)
                                                                                   aCTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGGACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACg                                                              <  3:365197/147‑1 (MQ=255)
                                                                                      tcctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTgg                                                                                                          >  6:350853/1‑100 (MQ=255)
                                                                                      tcctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTgg                                                                                                          <  5:350853/100‑1 (MQ=255)
                                                                                        ctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTaa                                                        <  3:61516/148‑1 (MQ=255)
                                                                                        ctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                         <  4:338540/147‑1 (MQ=255)
                                                                                        ctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                         <  7:116573/147‑1 (MQ=255)
                                                                                        ctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                         <  7:123273/147‑1 (MQ=255)
                                                                                         tccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAAt                                                       >  6:334832/1‑148 (MQ=255)
                                                                                          cccTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATg                                                      >  6:151495/1‑148 (MQ=255)
                                                                                           ccTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAAt                                                       <  2:21511/146‑1 (MQ=255)
                                                                                             tGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTc                                                   >  6:174217/1‑148 (MQ=255)
                                                                                              gATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcg                                                  >  8:159018/1‑148 (MQ=255)
                                                                                               aTCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcgg                                                 >  2:109777/1‑148 (MQ=255)
                                                                                                  aaCGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcggcg                                               <  5:174217/147‑1 (MQ=255)
                                                                                                   aCGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGTAAATACCGCGTTAAAATGTCCAGTTCGTATTCGCTCATCCTGGTTCGTGCACCCAATCCTATGGCAATGGAAGGCCACGCATTGACGGTGTAATGGTcggcgg                                              >  2:299002/1‑147 (MQ=255)
                                                                                                    cGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGGTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCggcggc                                             >  6:151802/1‑147 (MQ=255)
                                                                                                     gAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                         <  4:240878/134‑1 (MQ=255)
                                                                                                     gAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                         >  3:240878/1‑134 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCTCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCAAGAATTGACGGTGTAATCTTCGGCGTCAAt                                          <  6:440044/148‑1 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                          <  1:109777/148‑1 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                          <  7:306239/148‑1 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                          <  7:45741/148‑1 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                          <  5:151495/148‑1 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                          <  4:202286/148‑1 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGGCGGCGGCAAt                                          <  6:47097/148‑1 (MQ=255)
                                                                                                            aCTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGc                                    <  7:230443/148‑1 (MQ=255)
                                                                                                             cTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCa                                   >  3:411289/1‑148 (MQ=255)
                                                                                                               gTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGc                                                                     <  7:155656/112‑1 (MQ=255)
                                                                                                               gTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGc                                                                     >  8:155656/1‑112 (MQ=255)
                                                                                                               gTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCACCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCAt                                  >  7:247573/1‑147 (MQ=255)
                                                                                                                    cccgttataaCTTTGGCCGCGGACCGGCAATTCCGCGCTAAAATGTCCCCTTCCGTTTGCCCCCAGCCTGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGCCGGCAATAGGGGCATGACAAc                            <  4:149930/142‑1 (MQ=255)
                                                                                                                      ggTAATATCTTTGGCGTTGTACTGGAACTACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACaa                             <  7:429030/145‑1 (MQ=255)
                                                                                                                                         tccTGGTAATACCGCTTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGATCCTGCCCCCATTCCTTTGGCAATTGTCGGCCACGCTGTGCCGGTGTAATTGTCTGCGGCAATAGGGCCATGACAAGCAGCGCAATACGTATCGTCat       <  6:454796/146‑2 (MQ=255)
                                                                                                                                           cTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAg      <  4:334159/147‑1 (MQ=255)
                                                                                                                                            tGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGtt    >  1:253845/1‑148 (MQ=255)
                                                                                                                                             ggAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCTTCAAGttt   >  4:215975/1‑148 (MQ=255)
                                                                                                                                             ggAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGttt   >  4:209028/1‑148 (MQ=255)
                                                                                                                                              gAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGtttt  >  3:440441/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTT  >  NC_000913/1956926‑1957215

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: