Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,140,862 A→T 6.6% L142H (CTT→CAT)  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,8620AT 6.6% 321.2 / 9.2 121L142H (CTT→CAT) hybDmaturation protease for hydrogenase 2
Reads supporting (aligned to +/- strand):  ref base A (55/54);  new base T (3/5);  total (59/62)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

Only 100 of 121 total aligned reads displayed.

TGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGAT  >  NC_000913/3140714‑3141009
                                                                                                                                                    |                                                                                                                                                   
tGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAg                                                                                                                                                    >  2:893815/1‑150 (MQ=255)
  cATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAgcg                                                                                                                                                  <  2:1674767/150‑1 (MQ=255)
     cGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAg                                                                                                                                               <  1:218009/150‑1 (MQ=255)
     cGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAg                                                                                                                                               >  1:815349/1‑150 (MQ=255)
         cGCCGGGCAATTTCCTAAAACGCTGCCTGTACTTGCGCGTGCGGGGAGGGCTGGAACCCTGCTATCTCTTCAGTACTGAATCGCCTCCAGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCGCTTGAGCAGGTTc                                                                                                                                           <  1:1071613/150‑1 (MQ=255)
           ccGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCaa                                                                                                                                         >  2:638833/1‑150 (MQ=255)
             gggCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCGCTTGAGCAGGTTctggc                                                                                                                                       <  1:1559161/150‑5 (MQ=255)
              ggCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCa                                                                                                                                      >  2:1242435/1‑150 (MQ=255)
              ggCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCGCTTGAGCAGGTTctggct                                                                                                                                      <  1:613350/150‑6 (MQ=255)
                  aTTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGc                                                                                                                                  >  1:200402/1‑150 (MQ=255)
                  aTTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGc                                                                                                                                  >  1:32845/1‑150 (MQ=255)
                   tttCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGTATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCt                                                                                                                                 <  1:1194654/150‑1 (MQ=255)
                      cttcAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCa                                                                                                                              >  2:1304189/1‑150 (MQ=255)
                      cttcAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCa                                                                                                                              >  2:380725/1‑150 (MQ=255)
                      cttcAAAAGATGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCAAGAATCGCCTCCCGTGGGATGGCTTCAACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCa                                                                                                                              >  2:1304196/1‑150 (MQ=255)
                        tcAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCACGCGCAGGTTCAATCATTGCTTCAAc                                                                                                                            >  2:352178/1‑150 (MQ=255)
                          aaaCGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCg                                                                                                                          >  2:1215145/1‑150 (MQ=255)
                           aaCGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGt                                                                                                                         <  1:1734084/150‑1 (MQ=255)
                              gCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCgg                                                                                                                      <  1:814318/150‑1 (MQ=255)
                                 gCCTGAACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGt                                                                                                                   >  1:1725032/1‑150 (MQ=255)
                                  ccggTACTTGCGCCTTCGGGGAGGTCTGGAAACATGCTATCTCTTCAGTCATGAATCTCCTCCCGGGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGtt                                                                                                                  <  2:308294/147‑1 (MQ=255)
                                     ttACTTTCGCCTTCGGGGAGGTCTTGAAACCTGCTATCTCTTCAGTCAAGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGTCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAg                                                                                                               <  2:1071197/149‑1 (MQ=255)
                                      tACTTGAGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCTCCTCCCGTGGGATGGCTTCCACGCCAGATTCAAGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCTATCGTTGCTTCAACCGTCGGCGTTAAGc                                                                                                              <  1:1547021/150‑1 (MQ=255)
                                       aCTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGcc                                                                                                             <  2:777953/150‑1 (MQ=255)
                                       aCTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGcc                                                                                                             <  2:1656867/150‑1 (MQ=255)
                                         ttGCGCCTGCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGa                                                                                                           <  1:1118040/150‑1 (MQ=255)
                                           gcgcCTTCGGGGAGGTCTGGAAACCTGCTAGCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGAtg                                                                                                         <  1:1047582/150‑1 (MQ=255)
                                             gcCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGAtgtg                                                                                                       >  2:1674132/1‑150 (MQ=255)
                                              ccTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGAtgtgt                                                                                                      >  1:1410753/1‑150 (MQ=255)
                                               cTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGAtgtgtg                                                                                                     >  1:835173/1‑150 (MQ=255)
                                                   ggggAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTc                                                                                                 >  1:1364343/1‑150 (MQ=255)
                                                       aGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGc                                                                                             <  1:423352/150‑1 (MQ=255)
                                                          gcTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGAt                                                                                          <  2:1392322/149‑1 (MQ=255)
                                                           cTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAAAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGAtt                                                                                         >  1:342960/1‑150 (MQ=255)
                                                               aaaCCTGCTATCTCTTCAGTAATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGACCCTGCTCAAGCGCAGGTTCAATCTTTGCTGCTACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCgg                                                                                     <  2:1041707/150‑1 (MQ=255)
                                                                  ccTGCTATCTCTTCAGTCATGAATCGCCTGCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGAt                                                                                  <  1:107118/150‑1 (MQ=255)
                                                                   cTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATc                                                                                 >  1:389431/1‑150 (MQ=255)
                                                                     gCTATCTCTTCAGTCATGAATCGCCTCCATAGTGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCCGAACCTGCAAAAGCGCACGTACAATCATTGCTTCAACCGTCGGCTTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATAAc                                                                               >  2:1131683/1‑150 (MQ=255)
                                                                          ctctTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGa                                                                          >  2:215716/1‑150 (MQ=255)
                                                                           tctTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGAc                                                                         >  1:994114/1‑150 (MQ=255)
                                                                            ctTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCGACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGAcc                                                                        >  1:1068631/1‑150 (MQ=255)
                                                                            ctTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGccc                                                                        >  2:457855/1‑150 (MQ=255)
                                                                            ctTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGAcc                                                                        >  1:445166/1‑150 (MQ=255)
                                                                              tCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAg                                                                      <  1:1358869/150‑1 (MQ=255)
                                                                              tCAGTCATGAATCGCCACCCGTCGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCACGCGCAGGTTCAATCATTGCATCACCCGGCGGCGCTAAGCCGATGTGTGGCTACAGCGATTCCGGGATCACGCCGACCAg                                                                      <  1:841392/150‑1 (MQ=255)
                                                                                aGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAggg                                                                    <  1:1674132/150‑1 (MQ=255)
                                                                                   cATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCCCGCAGCGCAGCCAGAACCTCCTCAAGCGCAGGTACAATCATTGCTTCCACCGTCGGCGTTAAGCCCATGAGTGGCTCCAGCGATTCCGGGATCACGCCGACCAggggcc                                                                 >  1:1336579/1‑147 (MQ=255)
                                                                                    aTGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAg                                                                >  1:999868/1‑150 (MQ=255)
                                                                                      taaTCGCCTCCCGTGGGATGGCTGCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGTCGACCAGGGTCAGCt                                                              <  1:1085288/149‑1 (MQ=255)
                                                                                      caaTCGCCTCCCGGGGGATGGCGTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCt                                                              <  1:76827/149‑1 (MQ=255)
                                                                                       aaTCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCtt                                                             >  2:1343268/1‑150 (MQ=255)
                                                                                       aaTCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCtt                                                             >  2:1583669/1‑150 (MQ=255)
                                                                                           gCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGcc                                                                            >  1:419376/1‑131 (MQ=255)
                                                                                           gCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGcc                                                                            <  2:419376/131‑1 (MQ=255)
                                                                                            ccTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCg                                                        <  1:1207558/150‑1 (MQ=255)
                                                                                              tCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGa                                                      <  2:1246604/150‑1 (MQ=255)
                                                                                               cccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAGGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGaa                                                     >  2:1116637/1‑150 (MQ=255)
                                                                                               cccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGaa                                                     <  2:162435/150‑1 (MQ=255)
                                                                                               cccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGaa                                                     >  2:1275094/1‑150 (MQ=255)
                                                                                                ccGTGGGATGGGTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGaaa                                                    <  2:1123145/150‑1 (MQ=255)
                                                                                                  gTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACt                                                  <  1:1467883/150‑1 (MQ=255)
                                                                                                  gTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACt                                                  <  2:843734/150‑1 (MQ=255)
                                                                                                   tGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGTGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACta                                                 >  2:53781/1‑149 (MQ=255)
                                                                                                   ggggATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTc                                                 <  2:165470/149‑1 (MQ=255)
                                                                                                    gggATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCg                                                >  1:1541469/1‑150 (MQ=255)
                                                                                                      gATGGCTTCCACGCCAGATTCACGCAGCTCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGcc                                              <  1:1343268/150‑1 (MQ=255)
                                                                                                      gATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGcc                                              <  1:1583669/150‑1 (MQ=255)
                                                                                                         ggCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGt                                           <  1:179694/150‑1 (MQ=255)
                                                                                                         ggCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGGTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGt                                           <  2:193826/150‑1 (MQ=255)
                                                                                                          gCTTCCACGCCTGATTCACGCAGCGCAGCCAGAACGTGCTCAAGCGCAGGTTCAATCATTGCTTCAAACGTCGGCGTTAAGCCGATGGGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTg                                          <  1:1161387/150‑1 (MQ=255)
                                                                                                          gCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGACGGTg                                          >  1:864072/1‑150 (MQ=255)
                                                                                                           cTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCACCCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGCTTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTgc                                         >  1:1321623/1‑149 (MQ=255)
                                                                                                           cTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGa                                         >  2:1444891/1‑150 (MQ=255)
                                                                                                           cTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGa                                         >  1:50500/1‑150 (MQ=255)
                                                                                                           cTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGa                                         >  2:1591823/1‑150 (MQ=255)
                                                                                                             tCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAg                                       <  1:1242435/150‑1 (MQ=255)
                                                                                                                     aGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGc                               <  1:452962/150‑1 (MQ=255)
                                                                                                                           aCGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACaa                         <  1:1242285/150‑1 (MQ=255)
                                                                                                                           aCGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACaa                         >  2:68256/1‑150 (MQ=255)
                                                                                                                            cgcagcgcagCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGATGCGCCGGGCCGACAAg                        >  1:467097/1‑150 (MQ=255)
                                                                                                                             gcagcgcagcCAGAACCTGCTGATGCGCAGGTTCAATCATTGCTTCAACCGTCAGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGa                       <  1:1601898/150‑1 (MQ=255)
                                                                                                                               agcgcagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATAACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACg                     <  2:1559793/150‑1 (MQ=255)
                                                                                                                                     gcCAGAACCTGCTCAAGCGCAGGTTCAATTATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggcc               <  2:100556/150‑1 (MQ=255)
                                                                                                                                     gcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGCCCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggcc               >  2:64574/1‑150 (MQ=255)
                                                                                                                                     gattgAACCTGCGCTTGAGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggcc               >  2:1462757/5‑150 (MQ=255)
                                                                                                                                      ccAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggcca              <  1:1275094/150‑1 (MQ=255)
                                                                                                                                       cAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccag             >  1:1240253/1‑150 (MQ=255)
                                                                                                                                       cAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccag             <  2:431885/150‑1 (MQ=255)
                                                                                                                                        aGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGATACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccagg            >  2:697340/1‑150 (MQ=255)
                                                                                                                                          acccTGCTCACGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCGCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAACCTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccaggcc          <  1:1695819/148‑1 (MQ=255)
                                                                                                                                            ccTGCTCAAGCGCAAGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCaa        >  1:733232/1‑150 (MQ=255)
                                                                                                                                              tGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGA‑‑TGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTg    <  2:1728257/150‑1 (MQ=255)
                                                                                                                                              tGCTCAAGCGCAGGTTCAATCATGGCTTCACCCGGCGGCGTGAAGCCGATGGGTGGCTCCAGCGATGCCGGGATCACGCCGACCAGGGGCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGc      <  1:1290006/150‑1 (MQ=255)
                                                                                                                                               gCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACCCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCt     <  1:380725/150‑1 (MQ=255)
                                                                                                                                               gCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGGCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCt     <  1:680227/150‑1 (MQ=255)
                                                                                                                                                cTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTg    <  2:32845/150‑1 (MQ=255)
                                                                                                                                                cTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTg    <  2:200402/150‑1 (MQ=255)
                                                                                                                                                cTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTg    >  2:779131/1‑150 (MQ=255)
                                                                                                                                                cTCAAGCGCAGGTTCAATCATTGCTTCAACCGGCGGCGTTAAGCCGATGTGGGGCTCCAGCGATTCCGGGAGCACGCCGACCAGGGTCAGCTTTATCGGACACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTg    <  1:1153314/150‑1 (MQ=255)
                                                                                                                                                  cAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGAt  <  2:1725032/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGAT  >  NC_000913/3140714‑3141009

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: