Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,604,991 T→A 9.3% L61M (TTG→ATG) ‡ hyfD → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,604,9910TA 9.3% 211.3 / 10.8 75L61M (TTG→ATG) ‡hyfDhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (31/36);  new base A (4/3);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAAC  >  NC_000913/2604844‑2605129
                                                                                                                                                   |                                                                                                                                          
ccTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCacgttgacgtt                                                                                                                                           <  4:100236/149‑1 (MQ=255)
  ttttATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTgg                                                                                                                                         >  6:44023/1‑149 (MQ=255)
    ttATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGAGCTGACCACGCTGTGCATGTTGTCACTGATATCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGt                                                                                                                                        <  7:328235/148‑1 (MQ=255)
       tCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAAc                                                                                                                                    <  5:44023/149‑1 (MQ=255)
           cgcACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTgg                                                                                                                                <  2:58004/149‑1 (MQ=255)
             cACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTgggg                                                                                                                              <  7:120887/149‑1 (MQ=255)
             cACTGGTCGTTTCGTGTTCGCCACAACGGCGGGCCGCCGAATGGGGGGTGGTGTGCGCCGCGCTGCCCACGCGGTGCATGTTGGCACTGATCTCCGCGTTTTATCAGGCCGATAAAGGTGCCGTCACGTTGACGGTGGTCAACGTgggg                                                                                                                              <  8:14497/149‑1 (MQ=255)
                  ggcgtttcgtttTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGAtgt                                                                                                                         >  5:287920/3‑149 (MQ=255)
                  gTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGAtgt                                                                                                                         >  8:27025/1‑149 (MQ=255)
                  gTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGAtgt                                                                                                                         >  2:160494/1‑149 (MQ=255)
                   tcgtttcgttTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGAtgtg                                                                                                                        <  8:263549/149‑1 (MQ=255)
                             ttCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGAGCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGttt                                                                                                               <  4:178250/148‑1 (MQ=255)
                              tCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGGCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTg                                                                                                              <  7:28420/148‑1 (MQ=255)
                               cGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTgg                                                                                                             <  4:86454/148‑1 (MQ=255)
                               cGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGc                                                                                                            >  3:327992/1‑149 (MQ=255)
                                   acaACGTCGGGCCGCCGAATGTGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTgg                                                                                                        >  8:106500/1‑149 (MQ=255)
                                   acaACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTgg                                                                                                        >  7:212265/1‑149 (MQ=255)
                                     aaCGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGt                                                                                                       >  8:43975/1‑148 (MQ=255)
                                        gTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCAtt                                                                                                   <  7:229642/149‑1 (MQ=255)
                                             gccgccGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGAtgtg                                                                                                                        >  4:74804/1‑123 (MQ=255)
                                             gccgccGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGAtgtg                                                                                                                        <  3:74804/123‑1 (MQ=255)
                                                gccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGt                                                                                           <  3:98329/149‑1 (MQ=255)
                                                gccgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGt                                                                                           <  1:15707/149‑1 (MQ=255)
                                                  cgAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGa                                                                                         <  1:160494/149‑1 (MQ=255)
                                                      tGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTAc                                                                                     >  7:151637/1‑149 (MQ=255)
                                                      tGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTAc                                                                                     <  6:287920/149‑1 (MQ=255)
                                                           ggTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCAGTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGa                                                                                >  7:61998/1‑149 (MQ=255)
                                                                  ttCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGt                                                                          <  4:309664/148‑1 (MQ=255)
                                                                    cgccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATATAGTTTCCGTAACGTTTACGTTGTTCAAGTTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGGTGATTCTGttt                                                                        <  2:216963/148‑1 (MQ=255)
                                                                     gccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGTATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTtgt                                                                      >  6:36866/1‑149 (MQ=255)
                                                                     gccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTtgt                                                                      <  7:27025/149‑1 (MQ=255)
                                                                      ccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTtgtg                                                                     <  4:105849/149‑1 (MQ=255)
                                                                        gcgcTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGAtt                                                                              <  7:78746/138‑1 (MQ=255)
                                                                        gcgcTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGAtt                                                                              >  8:78746/1‑138 (MQ=255)
                                                                          gcTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGGGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTgtggtgg                                                                 <  1:17587/149‑1 (MQ=255)
                                                                              aCCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGtt                                                             >  2:165362/1‑149 (MQ=255)
                                                                              aCCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGtt                                                             >  6:106370/1‑149 (MQ=255)
                                                                                cACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTtc                                                           >  8:261922/1‑149 (MQ=255)
                                                                                 aCGCTGGGCTTGTTGTCACTGATCTGCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCTTTGATCGCGTGAGTACGCTGATTCTGGTTGTGGTGGTGGTtct                                                          <  8:128036/149‑1 (MQ=255)
                                                                                   gCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTTATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCg                                                        >  6:90720/1‑149 (MQ=255)
                                                                                   gCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCg                                                        <  5:106370/149‑1 (MQ=255)
                                                                                   gCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCg                                                        >  1:290780/1‑149 (MQ=255)
                                                                                       tgCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGttt                                                    >  2:149377/1‑149 (MQ=255)
                                                                                       tgCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGt                                                      <  7:106500/147‑1 (MQ=255)
                                                                                         cATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGc                                                  >  2:327963/1‑149 (MQ=255)
                                                                                           tgttgtCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGAt                                                                                                <  5:86505/101‑1 (MQ=255)
                                                                                           tgttgtCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGAt                                                                                                >  6:86505/1‑101 (MQ=255)
                                                                                             ttgtCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGt                                              <  5:209595/149‑1 (MQ=255)
                                                                                                   cTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATAGATCGCGTGAGTACGCTGATTCTGTATGTGGTGGTGTTTCTCGGTTAGCTGGTCACGAt                                        >  1:257365/1‑149 (MQ=255)
                                                                                                    tGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCt                                                                            >  5:285839/1‑112 (MQ=255)
                                                                                                    tGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCt                                                                            <  6:285839/112‑1 (MQ=255)
                                                                                                      aTCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATc                                       <  1:149377/147‑1 (MQ=255)
                                                                                                          ccGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGGGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTCTGTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTc                                  >  1:160213/1‑148 (MQ=255)
                                                                                                           cgcgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCCCGATCTACTCCa                                >  5:107358/1‑149 (MQ=255)
                                                                                                           cgcgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTAc                                    <  7:43975/145‑1 (MQ=255)
                                                                                                           cgcgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACt                                   <  4:327992/146‑1 (MQ=255)
                                                                                                            gcgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGTTCACGATCTACTCCAc                               >  6:232000/1‑149 (MQ=255)
                                                                                                            gcgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTTACGATCTACTCCAc                               >  4:242919/1‑149 (MQ=255)
                                                                                                              gTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACgg                             <  8:61998/149‑1 (MQ=255)
                                                                                                                  tATCAGGCCGATAAAGTTGCCGTCAGGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGtt                          <  2:290780/148‑1 (MQ=255)
                                                                                                                    tCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATc                       >  5:109134/1‑149 (MQ=255)
                                                                                                                      aGGCCGATAATGTTGCCGTCACGTTGACGTTTGTCAAGGTGGGGGATGTGGCGTTGTTTGGCCTGGTCTTTGATCGCGTGATGACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCt                      <  2:160213/148‑1 (MQ=255)
                                                                                                                      aGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGGTGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACgg                             >  1:86325/1‑141 (MQ=255)
                                                                                                                      aGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGGTGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACgg                             <  2:86325/141‑1 (MQ=255)
                                                                                                                        gCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                   <  7:58763/149‑1 (MQ=255)
                                                                                                                           gATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGa                >  6:153614/1‑149 (MQ=255)
                                                                                                                                      ccccACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATcgc                                                                                             <  3:212608/58‑1 (MQ=21)
                                                                                                                                      ccccACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATcgc                                                                                             >  4:212608/4‑61 (MQ=21)
                                                                                                                                      ccccACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATcgc                                                                                             >  8:19100/4‑61 (MQ=21)
                                                                                                                                      ccccACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATcgc                                                                                             <  7:19100/58‑1 (MQ=21)
                                                                                                                                      cGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATcgcg     <  1:165362/149‑1 (MQ=255)
                                                                                                                                      cGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATcgc      >  3:34198/1‑148 (MQ=255)
                                                                                                                                        tCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGaa   <  1:327963/149‑1 (MQ=255)
                                                                                                                                        ccACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGGTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                   <  4:179994/132‑1 (MQ=255)
                                                                                                                                        ccACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGAGTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                   >  3:179994/2‑133 (MQ=255)
                                                                                                                                         cACGTTGACCAACGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAac  >  6:202524/1‑149 (MQ=255)
                                                                                                                                            gttgacgttgGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGAtt                                                                              >  8:11666/1‑70 (MQ=255)
                                                                                                                                            gttgacgttgGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGAtt                                                                              <  7:11666/70‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGGGCCGCCGAATGGGGGGTTTTGTTCGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAAC  >  NC_000913/2604844‑2605129

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