Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,140,864 C→A 7.0% A141A (GCG→GCT)  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,8640CA 7.0% 341.9 / 12.6 114A141A (GCG→GCT) hybDmaturation protease for hydrogenase 2
Reads supporting (aligned to +/- strand):  ref base C (55/51);  new base A (5/3);  total (60/54)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

Only 100 of 119 total aligned reads displayed.

CATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGATG  >  NC_000913/3140716‑3141010
                                                                                                                                                    |                                                                                                                                                  
 aTCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAgcgc                                                                                                                                                  >  1:479833/1‑150 (MQ=255)
 aTCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAgcgc                                                                                                                                                  >  2:353378/1‑150 (MQ=255)
 aTCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCGGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAgcgc                                                                                                                                                  <  2:1122953/150‑1 (MQ=255)
  tCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCCCGCCAGATTCACGCAGCGCAGCCAGAACCTGCGCAAGCGCa                                                                                                                                                 >  2:1089860/1‑150 (MQ=255)
  tCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCTGATTCACGCAGCGCAGCCAGAACCAGCTCAAGCGCa                                                                                                                                                 >  1:297341/1‑150 (MQ=255)
  tCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCa                                                                                                                                                 <  2:1287708/150‑1 (MQ=255)
   cGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCAAGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAg                                                                                                                                                <  2:976306/150‑1 (MQ=255)
    gAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAgg                                                                                                                                               <  1:1251974/150‑1 (MQ=255)
    gAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAgg                                                                                                                                               <  1:543070/150‑1 (MQ=255)
       cGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTc                                                                                                                                            <  2:1153828/150‑1 (MQ=255)
         ccGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCGCTTGAGCAGGTTctg                                                                                                                                          <  2:1484178/150‑3 (MQ=255)
         ccGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCAATCTCTTCAGTCATGACTCGCCTCCCGTTGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCGCTTGAGCAGGTTctg                                                                                                                                          <  1:1565538/150‑3 (MQ=255)
          cGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAAt                                                                                                                                         >  2:654577/1‑150 (MQ=255)
          cGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCGCTTGAGCAGGTTctgg                                                                                                                                         <  1:1174051/150‑4 (MQ=255)
           tagaaaattcttcAAACGCTGAAAGTACTTGCGCCTTCGGGGAGGTAAGGAAAACTGATATCTCTTCAGTCATCAATCGCCTCCCTTGGTATGTCTACCAAGCCAGATTCACGCATCGCAGCCAGAACCTGCTCAAGCGCAGGTACAATc                                                                                                                                        >  2:198423/8‑150 (MQ=255)
           gggCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATc                                                                                                                                        <  2:535453/150‑1 (MQ=255)
               aaTTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTg                                                                                                                                    <  2:618545/150‑1 (MQ=255)
                  ttcttcAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCtt                                                                                                                                 <  2:1259070/150‑1 (MQ=255)
                       cAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAAcc                                                                                                                            >  2:582555/1‑150 (MQ=255)
                         aaCGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCGCGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCATTCATGGCTTCACCCGt                                                                                                                          >  1:137493/1‑150 (MQ=255)
                         aaCGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGt                                                                                                                          >  1:1081925/1‑150 (MQ=255)
                         aaCGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGt                                                                                                                          >  1:935019/1‑150 (MQ=255)
                         aaCGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGt                                                                                                                          >  2:786137/1‑150 (MQ=255)
                         aaCGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTAAACCgg                                                                                                                          >  1:1051766/1‑149 (MQ=255)
                           cGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCg                                                                                                                        <  2:540796/150‑1 (MQ=255)
                            gCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCgg                                                                                                                       <  1:1691331/150‑1 (MQ=255)
                             cTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGc                                                                                                                      <  2:1477067/150‑1 (MQ=255)
                                  tGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTaa                                                                                                                 >  1:406506/1‑150 (MQ=255)
                                    tACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGc                                                                                                               >  2:74298/1‑150 (MQ=255)
                                      cTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCg                                                                                                             >  1:799605/1‑150 (MQ=255)
                                       ttGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGa                                                                                                            <  2:3623/150‑1 (MQ=255)
                                          cgcCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGAtgt                                                                                                         >  2:746543/1‑150 (MQ=255)
                                           gcCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCATCGGCTTAACGCCGacggt                                                                                                        <  2:430595/150‑5 (MQ=255)
                                             cTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCATCGCAGCCAGAACCTGCTCAATCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGAtgtgtg                                                                                                      >  2:90062/1‑150 (MQ=255)
                                                cGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCt                                                                                                   <  1:423679/150‑1 (MQ=255)
                                                  gggAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTcc                                                                                                 >  2:403687/1‑150 (MQ=255)
                                                      ggTCTGGAAACCTGCTATCTCTTCAGTCATGAAACGCCTCCCGTGGGATGGCTTCCACGCCTGAATCACGCAGCGCAGCCAGAACCTGCTCAACCGCAGGTTCAATCATTCATTCAACCGTCCGCGTTAAGCCCAGCTGTGGATCCAGcc                                                                                             >  1:300230/1‑149 (MQ=255)
                                                            gAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGAGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCg                                                                                       >  1:991208/1‑150 (MQ=255)
                                                              aaCCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCggg                                                                                     <  1:7260/150‑1 (MQ=255)
                                                               aCCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGa                                                                                    <  1:800966/150‑1 (MQ=255)
                                                                ccTGCTATCTCTCCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGAt                                                                                   >  2:1116210/1‑150 (MQ=255)
                                                                  tGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTg                                                                                                                                    >  1:543582/1‑99 (MQ=255)
                                                                  tGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCa                                                                                 <  2:420303/150‑1 (MQ=255)
                                                                  tGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCa                                                                                 >  1:420303/1‑150 (MQ=255)
                                                                    cTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACg                                                                               >  2:1653345/1‑150 (MQ=255)
                                                                     tATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGc                                                                              <  1:1632268/150‑1 (MQ=255)
                                                                     tATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGc                                                                              >  1:657540/1‑150 (MQ=255)
                                                                      aTCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCATGCCAGATTCACGCAGCGCAGCCAGAACCTGCACAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGcc                                                                             <  1:696352/150‑1 (MQ=255)
                                                                      aTCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGcc                                                                             <  2:1081925/150‑1 (MQ=255)
                                                                      aTCTCTTCAGTAATGAATCGCCTCCCGTGGGATGGCTACCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGcc                                                                             <  2:137493/150‑1 (MQ=255)
                                                                        ctctTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAGGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGa                                                                           >  1:677672/1‑150 (MQ=255)
                                                                        ctctTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGa                                                                           <  1:570092/150‑1 (MQ=255)
                                                                        ctctTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGa                                                                           >  2:491603/1‑150 (MQ=255)
                                                                           ttCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCATCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCa                                                                        <  2:1512976/150‑1 (MQ=255)
                                                                              cgTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAggg                                                                     <  1:369206/149‑1 (MQ=255)
                                                                              aGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAggg                                                                     >  2:566392/1‑150 (MQ=255)
                                                                               gTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGt                                                                    >  2:682704/1‑150 (MQ=255)
                                                                                tCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTc                                                                   <  1:1644962/150‑1 (MQ=255)
                                                                                 cATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTcc                                                                  >  2:987314/1‑149 (MQ=255)
                                                                                 cATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGAGTCACGCAGCGCAGCCAGAACCTGCGCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGGTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCa                                                                  >  1:836175/1‑150 (MQ=255)
                                                                                    gAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTGCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCt                                                               <  1:58758/150‑1 (MQ=255)
                                                                                         gCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTc                                                          >  1:1258938/1‑150 (MQ=255)
                                                                                          ccTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCg                                                         <  2:390371/150‑1 (MQ=255)
                                                                                          ccTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCg                                                         <  1:74298/150‑1 (MQ=255)
                                                                                             cccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCCAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAACCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGCCCAGGGTCAGCTTTTTCGGaa                                                      >  1:1286473/1‑150 (MQ=255)
                                                                                             cccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGaa                                                      >  1:582383/1‑150 (MQ=255)
                                                                                              ccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGaaa                                                     <  2:71020/150‑1 (MQ=255)
                                                                                               cGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGAGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAAc                                                    <  2:158210/150‑1 (MQ=255)
                                                                                               cGTGGGATGGCTTCCACACCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTTCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAAc                                                    >  1:1170823/1‑150 (MQ=255)
                                                                                                        gCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTCAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTg                                           >  2:358553/1‑150 (MQ=255)
                                                                                                        gCTTCCACGCCAGACTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGCCCAGGGTCAGCTTTTTCGGAAACTCGCCGGTg                                           >  1:1604613/1‑150 (MQ=255)
                                                                                                         cTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGa                                          >  2:905940/1‑150 (MQ=255)
                                                                                                           gccACGCCAGATGCACGCAGCGCAGCCAGTACCTGCTCAAGCGCAGGTTCAATCTTTGCTTCTACCGTTGGCGTTATGCCGATGTGTGGCTCCAGCGATTCGGGGATCACGCCGACCAGGGTCAGCTTTTTCTGAAACTCGCCGGTGAAg                                        <  1:79372/149‑1 (MQ=255)
                                                                                                            ccACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCGTTTTCGGAAACTCGCCGGTGAAgc                                       <  2:510560/150‑1 (MQ=255)
                                                                                                              aCGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAgcgc                                     <  2:1258938/150‑1 (MQ=255)
                                                                                                               cGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCa                                    <  2:1557129/150‑1 (MQ=255)
                                                                                                                gCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCACCCGTCGGCGTTACGCCGATGTGTGGCTCCAGCGATTCCGGGCTCACGCCGACCAGGGTCCGCTTTTTCGGAAACTCGCCGGTGAAGCGCAg                                   >  1:74679/1‑150 (MQ=255)
                                                                                                                         aCGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACaa                          >  1:951394/1‑150 (MQ=255)
                                                                                                                         aCGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACaa                          >  1:266213/1‑150 (MQ=255)
                                                                                                                          cgcagcgcagCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAACCTCGCCGGTGAAGCGCAGGGCCGACAAg                         >  1:932087/1‑150 (MQ=255)
                                                                                                                          cgcagcgcagCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAg                         >  1:144734/1‑150 (MQ=255)
                                                                                                                          cgcagcgcagCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAg                         >  1:1500176/1‑150 (MQ=255)
                                                                                                                                 cagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCgg                  >  1:1475366/1‑150 (MQ=255)
                                                                                                                                 cagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTCAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCgg                  <  1:358553/150‑1 (MQ=255)
                                                                                                                                  agcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggc                 <  1:403687/150‑1 (MQ=255)
                                                                                                                                   gcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGATACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggcc                >  1:1186083/1‑150 (MQ=255)
                                                                                                                                    ccAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggcca               >  1:153786/1‑150 (MQ=255)
                                                                                                                                    ccAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggcca               >  1:797780/1‑150 (MQ=255)
                                                                                                                                    attgAACCTGCGCTTGAGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggcca               >  2:1479615/4‑150 (MQ=255)
                                                                                                                                      aGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTGCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccagg             >  1:87555/1‑150 (MQ=255)
                                                                                                                                      aGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGGTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGGCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccagg             <  1:92736/150‑1 (MQ=255)
                                                                                                                                       gAACCTGCTCAAGCGCAGGGTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccaggc            <  1:705272/150‑1 (MQ=255)
                                                                                                                                       gAACCTGCGCTTGAGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccaggc            >  2:590919/1‑150 (MQ=255)
                                                                                                                                       gAACCTGCGCTTGAGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccaggc            >  2:1668785/1‑150 (MQ=255)
                                                                                                                                       gAACCTGCGCTTGAGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccaggc            >  1:1369049/1‑150 (MQ=255)
                                                                                                                                        aaCCTGCGCTTGAGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccaggcc           >  1:1406898/1‑150 (MQ=255)
                                                                                                                                          ccTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCaa         >  2:197688/1‑150 (MQ=255)
                                                                                                                                            tGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGc       <  2:804981/150‑1 (MQ=255)
                                                                                                                                                cAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGAGCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGAt   <  1:353378/150‑1 (MQ=255)
                                                                                                                                                 aaGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGATg  >  1:1038586/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGATG  >  NC_000913/3140716‑3141010

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