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

GGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTT  >  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)
                        ccaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGtt                                                                                                                        >  8:123273/1‑148 (MQ=255)
                        ccaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCAGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGtt                                                                                                                        >  8:116573/1‑148 (MQ=255)
                         caGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGttt                                                                                                                       >  1:457088/1‑148 (MQ=255)
                         caGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGttt                                                                                                                       <  8:203108/148‑1 (MQ=255)
                            cgcTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCGGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCg                                                                                                                    >  8:127521/1‑148 (MQ=255)
                             gcTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGc                                                                                                                   <  5:288104/148‑1 (MQ=255)
                              cTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCGGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCt                                                                                                                  <  7:127521/148‑1 (MQ=255)
                                ccAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTAAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTTTTTCTCGGTAATATCTTTGGCGTTGTACTGGATATACACAGTTAAAATGTCCATTT‑CGTTTTCTCTCa                                                                                                                >  8:230443/1‑148 (MQ=255)
                                 cAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGCAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTc                                                                                                                 <  1:61751/146‑1 (MQ=255)
                                 cAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGCAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTc                                                                                                                 <  1:55003/146‑1 (MQ=255)
                                 cAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTc                                                                                                                 <  2:457088/146‑1 (MQ=255)
                                 cAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTc                                                                                                                 <  8:322420/146‑1 (MQ=255)
                                   gCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGc                                                                                                             >  1:195338/1‑148 (MQ=255)
                                         aTACCGCTGTGTTTAATATATTCACGTCTAGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTGCTCGGTAATATCTTTGGCGTTGTACTGGAAATAGCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCTTGCTGTTTc                                                                                                       <  6:20783/148‑1 (MQ=255)
                                         aTACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGGAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGGT‑CGTTTTCGCTCATGCTGGTTc                                                                                                       >  4:54507/1‑148 (MQ=255)
                                           aCCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTg                                                                                                           <  5:399479/142‑1 (MQ=255)
                                           aCCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTg                                                                                                           >  6:399479/1‑142 (MQ=255)
                                           aCCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGt                                                                                                     >  8:56406/1‑148 (MQ=255)
                                           aCCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATAACTGTTTATAGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTG‑AGTTTTCTCTCa                                                                                                                >  1:382443/1‑137 (MQ=255)
                                             cGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGc                                                                                                   <  2:195338/148‑1 (MQ=255)
                                                 gtgtTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGcccc                                                                                                >  8:306239/1‑147 (MQ=255)
                                                      taataaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGGTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTGTTCGCTCATGCTGGTTCGTGCCCCCATTcc                                                                                          >  2:113185/1‑148 (MQ=255)
                                                         taaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCTTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCttt                                                                                       >  4:360400/1‑148 (MQ=255)
                                                             ttCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCa                                                                                   <  5:280938/148‑1 (MQ=255)
                                                             ttCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCa                                                                                   <  3:15527/148‑1 (MQ=255)
                                                                 cGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATgg                                                                               <  2:428022/148‑1 (MQ=255)
                                                                          aaTGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCAc                                                                       <  8:389757/147‑1 (MQ=255)
                                                                           aTGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGt                                                                                                                                           <  4:332895/79‑1 (MQ=255)
                                                                           aTGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGt                                                                                                                                           >  3:332895/1‑79 (MQ=255)
                                                                              tCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCAtt                                                                  >  6:267478/1‑148 (MQ=255)
                                                                               cATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTg                                                                 >  2:47451/1‑148 (MQ=255)
                                                                                  aaCTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACg                                                              <  1:308087/148‑1 (MQ=255)
                                                                                   aCTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACgg                                                             >  6:416189/1‑148 (MQ=255)
                                                                                   aCTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGGACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACg                                                              <  3:365197/147‑1 (MQ=255)
                                                                                      tcctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTgg                                                                                                          >  6:350853/1‑100 (MQ=255)
                                                                                      tcctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTgg                                                                                                          <  5:350853/100‑1 (MQ=255)
                                                                                        ctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTaa                                                        <  3:61516/148‑1 (MQ=255)
                                                                                        ctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                         <  4:338540/147‑1 (MQ=255)
                                                                                        ctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                         <  7:116573/147‑1 (MQ=255)
                                                                                        ctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                         <  7:123273/147‑1 (MQ=255)
                                                                                         tccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAAt                                                       >  6:334832/1‑148 (MQ=255)
                                                                                          cccTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATg                                                      >  6:151495/1‑148 (MQ=255)
                                                                                           ccTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAAt                                                       <  2:21511/146‑1 (MQ=255)
                                                                                             tGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTc                                                   >  6:174217/1‑148 (MQ=255)
                                                                                              gATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcg                                                  >  8:159018/1‑148 (MQ=255)
                                                                                               aTCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcgg                                                 >  2:109777/1‑148 (MQ=255)
                                                                                                  aaCGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcggcg                                               <  5:174217/147‑1 (MQ=255)
                                                                                                   aCGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGTAAATACCGCGTTAAAATGTCCAGTT‑CGTATTCGCTCATCCTGGTTCGTGCACCCAATCCTATGGCAATGGAAGGCCACGCATTGACGGTGTAATGGTcggcgg                                              >  2:299002/1‑147 (MQ=255)
                                                                                                    cGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGGT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCggcggc                                             >  6:151802/1‑147 (MQ=255)
                                                                                                     gAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                         <  4:240878/134‑1 (MQ=255)
                                                                                                     gAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                         >  3:240878/1‑134 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCTCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCAAGAATTGACGGTGTAATCTTCGGCGTCAAt                                          <  6:440044/148‑1 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                          <  1:109777/148‑1 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                          <  7:306239/148‑1 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                          <  7:45741/148‑1 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                          <  5:151495/148‑1 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                          <  4:202286/148‑1 (MQ=255)
                                                                                                      aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGGCGGCGGCAAt                                          <  6:47097/148‑1 (MQ=255)
                                                                                                            aCTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGc                                    <  7:230443/148‑1 (MQ=255)
                                                                                                             cTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCa                                   >  3:411289/1‑148 (MQ=255)
                                                                                                               gTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGc                                                                     <  7:155656/112‑1 (MQ=255)
                                                                                                               gTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGc                                                                     >  8:155656/1‑112 (MQ=255)
                                                                                                               gTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCACCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCAt                                  >  7:247573/1‑147 (MQ=255)
                                                                                                                    cccgttataaCTTTGGCCGCGGACCGGCAATTCCGCGCTAAAATGTCCCCTTCCGTTTGCCCCCA‑GCCTGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGCCGGCAATAGGGGCATGACAAc                            <  4:149930/142‑1 (MQ=255)
                                                                                                                      ggTAATATCTTTGGCGTTGTACTGGAACTACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACaa                             <  7:429030/145‑1 (MQ=255)
                                                                                                                                         tccTGGTAATACCGCTTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGATCCTGCCCCCATTCCTTTGGCAATTGTCGGCCACGCTGTGCCGGTGTAATTGTCTGCGGCAATAGGGCCATGACAAGCAGCGCAATACGTATCGTCat       <  6:454796/146‑2 (MQ=255)
                                                                                                                                           cTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAg      <  4:334159/147‑1 (MQ=255)
                                                                                                                                            tGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGtt    >  1:253845/1‑148 (MQ=255)
                                                                                                                                             ggAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCTTCAAGttt   >  4:215975/1‑148 (MQ=255)
                                                                                                                                             ggAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGttt   >  4:209028/1‑148 (MQ=255)
                                                                                                                                              gAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGtttt  >  3:440441/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTT‑CGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTT  >  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: ‑