Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,140,844 G→A 6.2% A148V (GCG→GTG)  hybD ← maturation protease for hydrogenase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,140,8440GA 6.2% 202.0 / 2.8 64A148V (GCG→GTG) hybDmaturation protease for hydrogenase 2
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base A (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

CAGAAAAGAGAGATCGTGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGT  >  NC_000913/3140698‑3140991
                                                                                                                                                  |                                                                                                                                                   
cagtcaaGAGAGATCGTGCATCGAACGCCGGGCATTTTCTTCAAACGCTGACTGTACTTGCGCCTGCGGGGAGGTCGGGAAACCTGCTCTCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCAcgcagcgc                                                                                                                                                    <  4:259520/143‑1 (MQ=255)
   aaaaGAGAGATCGTGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTAGGGGAGGTCTGGAATCCTGGTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACgcagcgcagc                                                                                                                                                 >  2:267707/1‑148 (MQ=255)
          agaTCGTGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAAcc                                                                                                                                          <  8:368074/148‑1 (MQ=255)
                tGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTc                                                                                                                                     >  6:348294/1‑147 (MQ=255)
                   aTCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAgc                                                                                                                                 <  1:400011/148‑1 (MQ=255)
                      gAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAgcgc                                                                                                                               <  4:348215/147‑1 (MQ=255)
                      gAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCa                                                                                                                              <  5:32679/148‑1 (MQ=255)
                      gAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCa                                                                                                                              <  6:402866/148‑1 (MQ=255)
                       aaCGCCGGGCAATTTCTTCAAACGCTGCCTGTAGTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAAGCTGCTCAAGCGCAg                                                                                                                             <  8:190071/148‑1 (MQ=255)
                            cGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCa                                                                                                                        <  7:75893/148‑1 (MQ=255)
                              ggCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGAATCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAAt                                                                                                                      <  1:991/148‑1 (MQ=255)
                               gCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGCTTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAAt                                                                                                                      >  2:31550/1‑147 (MQ=255)
                                         cAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCaa                                                                                                           >  1:455475/1‑148 (MQ=255)
                                         cAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCa                                                                                                            <  5:327512/147‑1 (MQ=255)
                                            aCGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAAcc                                                                                                         >  5:339569/1‑147 (MQ=255)
                                             cGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGt                                                                                                       <  7:330046/148‑1 (MQ=255)
                                               cTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCATGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCg                                                                                                     >  7:11835/1‑148 (MQ=255)
                                                  ccTGTACTTGCGCCTTCGGGGAGGTCTGGAAACGTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCg                                                                                                  <  7:92419/148‑1 (MQ=255)
                                                       aCTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCACAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAg                                                                                             >  8:100426/1‑148 (MQ=255)
                                                       aCTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCACAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAg                                                                                             >  8:107038/1‑148 (MQ=255)
                                                        cTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGc                                                                                            <  3:145341/148‑1 (MQ=255)
                                                          tGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGcc                                                                                           >  6:296336/1‑147 (MQ=255)
                                                          tGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGcc                                                                                           >  7:248078/1‑147 (MQ=255)
                                                           gcgcCTTC‑GGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGAt                                                                                        >  6:285571/1‑148 (MQ=255)
                                                            cgcCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGa                                                                                         <  7:306173/147‑1 (MQ=255)
                                                              ccTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGAtgt                                                                                      >  8:402399/1‑148 (MQ=255)
                                                               cTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGAtgtg                                                                                     >  5:269259/1‑148 (MQ=255)
                                                                 tCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGACAGATTCACGCAGCGCAGCCAGa                                                                                                                                             >  7:70023/1‑90 (MQ=255)
                                                                 tCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGACAGATTCACGCAGCGCAGCCAGa                                                                                                                                             <  8:70023/90‑1 (MQ=255)
                                                                              gAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGa                                                                         <  8:301926/145‑1 (MQ=255)
                                                                              gAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTc                                                                      <  5:20615/148‑1 (MQ=255)
                                                                                  ccTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCgg                                                                   >  2:202375/1‑147 (MQ=255)
                                                                                   cTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGa                                                                 >  4:35795/1‑148 (MQ=255)
                                                                                     gcgaTCTCTGCCGGCATGAATCGCCTCCCGGGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCGCAAGCGCAGGTTCAATCATTGCTTCAACCGGCGGCGTTAAGCCGATGGGTGGCTCCAGCGATTCCGGGATc                                                               <  6:309119/145‑1 (MQ=255)
                                                                                      cTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATc                                                               <  1:187805/147‑1 (MQ=255)
                                                                                       tATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCAc                                                             >  2:147716/1‑148 (MQ=255)
                                                                                          ctctTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGcc                                                          >  1:292055/1‑148 (MQ=255)
                                                                                            ctTCAGTCATGAATCACCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGa                                                        >  7:21608/1‑148 (MQ=255)
                                                                                             ttCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCATGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGc                                                           <  8:11835/144‑1 (MQ=255)
                                                                                             ttCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCACAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCg                                                         <  7:100426/146‑1 (MQ=255)
                                                                                             ttCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCACAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCg                                                         <  7:107038/146‑1 (MQ=255)
                                                                                              tCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGAc                                                       >  4:299144/1‑147 (MQ=255)
                                                                                                aGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGt                                                                                                       >  1:175824/1‑97 (MQ=255)
                                                                                                aGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGt                                                                                                       <  2:175824/97‑1 (MQ=255)
                                                                                                   cATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAg                                                    >  1:449973/1‑145 (MQ=255)
                                                                                                   cATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAg                                                    <  2:449973/145‑1 (MQ=255)
                                                                                                   cATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGt                                                 >  1:224011/1‑148 (MQ=255)
                                                                                                   cATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGt                                                 >  7:80333/1‑148 (MQ=255)
                                                                                                    aTGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTc                                                <  1:234931/148‑1 (MQ=255)
                                                                                                         tCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCtt                                           >  3:371555/1‑148 (MQ=255)
                                                                                                             cTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTc                                       >  1:146058/1‑148 (MQ=255)
                                                                                                             cTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTc                                       >  8:216894/1‑148 (MQ=255)
                                                                                                               cccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCACCCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGCTTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCg                                      >  2:273001/1‑147 (MQ=255)
                                                                                                                ccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGACCCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCg                                      <  7:407192/146‑1 (MQ=255)
                                                                                                                ccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCCAGCGCAGGTTCAATCATTGCTTGAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCg                                      <  1:267707/146‑1 (MQ=255)
                                                                                                                         ggCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCg                           >  1:431323/1‑148 (MQ=255)
                                                                                                                             tCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTg                        >  3:427500/1‑147 (MQ=255)
                                                                                                                              ccACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTg                        >  4:259787/1‑146 (MQ=255)
                                                                                                                                aCGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAgc                    >  5:38742/1‑148 (MQ=255)
                                                                                                                                  gCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGTCCAGGGTCAGCTTTTTCTGAAACTCGCCGGTGAAgcgc                  <  6:339569/148‑1 (MQ=255)
                                                                                                                                  gCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAgcgc                  <  2:292055/148‑1 (MQ=255)
                                                                                                                                   ccAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCa                 <  4:427500/148‑1 (MQ=255)
                                                                                                                                     aGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAgg               <  7:216894/148‑1 (MQ=255)
                                                                                                                                        ttCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGcc            >  6:29244/1‑148 (MQ=255)
                                                                                                                                          cACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTGACGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTGTTTCGGAACCTCGCCGGTGAAGCGCAGGGCCGa          >  6:360499/1‑148 (MQ=255)
                                                                                                                                           aCGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGAGGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGAc         <  5:296336/148‑1 (MQ=255)
                                                                                                                                            cgcagcgcagCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGAc         >  8:132099/1‑147 (MQ=255)
                                                                                                                                            cgcagcgcagCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGAc         >  1:209528/1‑147 (MQ=255)
                                                                                                                                             gcagcgcagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACaa       <  4:247059/148‑1 (MQ=255)
                                                                                                                                              cagcgcagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAg      <  1:202375/148‑1 (MQ=255)
                                                                                                                                                gcgcagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGAcgtctt    <  2:455475/148‑6 (MQ=255)
                                                                                                                                                  gcagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGt  <  2:209528/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAGAAAAGAGAGATCGTGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGT  >  NC_000913/3140698‑3140991

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