Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,607,935 C→A 6.5% A340E (GCG→GAG)  hyfF → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,607,9350CA 6.5% 195.7 / 3.4 61A340E (GCG→GAG) hyfFhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (25/31);  new base A (2/2);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

GGCATGTTGTCGGTTGCCGTGGCGGCATTTTTCATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGG  >  NC_000913/2607793‑2608079
                                                                                                                                              |                                                                                                                                                
ggCATGTTGTCGGTTGCCGTGGCGGCATTTTTCATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCgctgc                                                                                                                                             <  8:291307/148‑1 (MQ=255)
        gtCGGTTGCCGTGGCGGCATTTTTCATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATc                                                                                                                                     <  6:19595/148‑1 (MQ=255)
        gtCGGTTGCCGTGGCGGCATTTTTCATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCt                                                                                                                                    >  1:228497/1‑149 (MQ=255)
         tCGGTTGCCGTGGCGGCATTTTTCATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCt                                                                                                                                    <  4:123986/148‑1 (MQ=255)
           ggTTGCCGTGGCGGCATTTTTCATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTaa                                                                                                                                 <  6:314268/149‑1 (MQ=255)
             ttGCCGTGGCGGCATTTTTCATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAAc                                                                                                                               <  1:136406/149‑1 (MQ=255)
                 cGTGGCGGCATTTTTCATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACca                                                                                                                             >  3:353277/1‑147 (MQ=255)
                  gTGGCGGCATTTTTCATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAg                                                                                                                          >  6:209160/1‑149 (MQ=255)
                       ggcATTTTTCATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTgg                                                                                                                     >  8:217699/1‑149 (MQ=255)
                        gcATTTTTCATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTgg                                                                                                                     >  8:29952/1‑148 (MQ=255)
                                cATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACgc                                                                                                            >  5:338584/1‑149 (MQ=255)
                                cATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACg                                                                                                             <  5:209160/148‑1 (MQ=255)
                                 aTTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACgct                                                                                                           >  1:336636/1‑149 (MQ=255)
                                  ttCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACgctg                                                                                                          >  3:35848/1‑149 (MQ=255)
                                     tGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGctgctg                                                                                                       >  5:291261/1‑149 (MQ=255)
                                       gTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGtt                                                                                                     <  6:338584/149‑1 (MQ=255)
                                       gTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGtt                                                                                                     >  7:277436/1‑149 (MQ=255)
                                         aCAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCt                                                                                                   >  1:399908/1‑149 (MQ=255)
                                          cAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTg                                                                                                  <  6:291261/149‑1 (MQ=255)
                                          cAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTg                                                                                                  >  2:713/1‑149 (MQ=255)
                                           aGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGtt                                                                                                     <  4:267238/145‑1 (MQ=255)
                                           aGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGtt                                                                                                     >  3:267238/1‑145 (MQ=255)
                                             cGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCgg                                                                                               <  4:35848/149‑1 (MQ=255)
                                                gACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTc                                                                                            <  2:399908/149‑1 (MQ=255)
                                                 aCATTTAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTTGGCGCGGTGTTTCTAGGGATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACATTCTGGCAAAAACGCTGCTGTTCTGCGGTTcc                                                                                           <  8:96006/149‑1 (MQ=255)
                                                  cATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCg                                                                                          <  8:318398/149‑1 (MQ=255)
                                                      aaGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCa                                                                                       <  2:336636/148‑1 (MQ=255)
                                                         cGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATg                                                                                    >  7:43322/1‑148 (MQ=255)
                                                          gTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTa                                                                                  <  1:112994/149‑1 (MQ=255)
                                                           tCTGCTGTCGTACTTCAGCGTGGAGAAGATGGTGATGGTCGCGGTGGATCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGGGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTAc                                                                                 <  4:353277/149‑1 (MQ=255)
                                                           tCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGATAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTAc                                                                                 >  6:243665/1‑149 (MQ=255)
                                                                 ggCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCgcgc                                                                                                                                                <  1:192838/80‑1 (MQ=255)
                                                                 ggCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCgcgc                                                                                                                                                >  2:192838/1‑80 (MQ=255)
                                                                 ggCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTc                                                                            >  7:124427/1‑148 (MQ=255)
                                                                  gCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTc                                                                            <  3:299101/147‑1 (MQ=255)
                                                                    gTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAg                                                                         >  2:36544/1‑148 (MQ=255)
                                                                        tCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACgg                                                                    >  3:232656/1‑149 (MQ=255)
                                                                        tCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACgg                                                                    >  2:325110/1‑149 (MQ=255)
                                                                           aGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCAc                                                                 <  1:281835/149‑1 (MQ=255)
                                                                            gCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGAGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTCCTGCTCAAGTACGGCAcg                                                                >  2:43537/1‑149 (MQ=255)
                                                                            gCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGAGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCAcg                                                                >  2:43532/1‑149 (MQ=255)
                                                                                ggAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCAcgcgc                                                             >  8:82798/1‑148 (MQ=255)
                                                                                   gaACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGAtc                                                         <  7:29952/149‑1 (MQ=255)
                                                                                       aTGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCaa                                                     >  2:413023/1‑149 (MQ=255)
                                                                                             cTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAAcgtcgt                                               <  8:43322/149‑1 (MQ=255)
                                                                                             cTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGAGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAAcgtcgt                                               <  1:43537/149‑1 (MQ=255)
                                                                                             cTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGAGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAAcgtcgt                                               <  1:43532/149‑1 (MQ=255)
                                                                                               ggTCGCGGGGGAGCTAGGCATTGGCGGGCCGCTGGGAATTGTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCt                                             <  4:282819/149‑1 (MQ=255)
                                                                                                 tCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCtgt                                           <  1:713/149‑1 (MQ=255)
                                                                                                 tCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCtgt                                           <  1:325110/149‑1 (MQ=255)
                                                                                                       tGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCtgtg                                          <  7:121210/144‑1 (MQ=255)
                                                                                                       tGGAGCTAGGCATAGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCtgtg                                          >  8:121210/1‑144 (MQ=255)
                                                                                                                gCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTg                                                                                                  <  5:283425/79‑1 (MQ=255)
                                                                                                                gCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTAAAAAt                             >  6:128828/1‑148 (MQ=255)
                                                                                                                gCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCAGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTg                                                                                                  >  6:283425/1‑79 (MQ=255)
                                                                                                                 cATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCa                           <  4:232656/149‑1 (MQ=255)
                                                                                                                     ggCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGAGCTCAACGTCGTCTGTGGGATGCTCAAAATCATGc                        <  6:373394/148‑1 (MQ=255)
                                                                                                                      gCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGc                                                                                                 <  1:159619/74‑1 (MQ=255)
                                                                                                                      gCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGc                                                                                                 >  2:159619/1‑74 (MQ=255)
                                                                                                                      gCGGGCCGCGGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCa                      <  6:243845/149‑1 (MQ=255)
                                                                                                                                 ggAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTAccgccg            <  7:82798/148‑1 (MQ=255)
                                                                                                                                         tGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGttt    <  5:243665/148‑1 (MQ=255)
                                                                                                                                           ccGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTgg  <  2:228497/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GGCATGTTGTCGGTTGCCGTGGCGGCATTTTTCATTCTGGTACAGCGGGACATTAAGCGTCTGCTGGCGTACTCCAGCGTGGAGAACATGGGGCTGGTCGCGGTGGAGCTAGGCATTGGCGGGCCGCTGGGAATTTTTGCCGCGCTGCTGCACATCTTAAACCACAGTCTGGCAAAAACGCTGCTGTTCTGCGGTTCCGGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGG  >  NC_000913/2607793‑2608079

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: ‑