Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,140,861 A→T 6.3% L142L (CTT→CTA) ‡ 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,8610AT 6.3% 340.4 / 8.3 120L142L (CTT→CTA) ‡hybDmaturation protease for hydrogenase 2
Reads supporting (aligned to +/- strand):  ref base A (57/52);  new base T (5/3);  total (63/57)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 123 total aligned reads displayed.

CGTGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGATG  >  NC_000913/3140712‑3141010
                                                                                                                                                     |                                                                                                                                                     
cGTGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCa                                                                                                                                                       >  1:3623/1‑150 (MQ=255)
  tGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAg                                                                                                                                                     <  2:914867/150‑1 (MQ=255)
  tGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAg                                                                                                                                                     <  2:1079209/150‑1 (MQ=255)
  tGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAg                                                                                                                                                     <  1:1597763/150‑1 (MQ=255)
    cATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAgcg                                                                                                                                                   >  1:390371/1‑150 (MQ=255)
    cATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAgcg                                                                                                                                                   >  1:510560/1‑150 (MQ=255)
     aTCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAgcgc                                                                                                                                                  >  1:479833/1‑150 (MQ=255)
     aTCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAgcgc                                                                                                                                                  >  2:353378/1‑150 (MQ=255)
     aTCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCGGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAgcgc                                                                                                                                                  <  2:1122953/150‑1 (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)
         aaCTCCGGGCATTTTCTTCAAACGCTGCGTTTAGATCCGGCTGCGGGGAGGGCTGGAACCCTGCTATCTCTTCATTCATGATTCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCTCAGCCAGAACCTGCTCCATCGCAGGt                                                                                                                                              <  2:667911/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)
                  cAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTg                                                                                                                                    <  2:795596/150‑1 (MQ=255)
                   aaTTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTg                                                                                                                                    <  2:618545/150‑1 (MQ=255)
                      ttcttcAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCtt                                                                                                                                 <  2:1259070/150‑1 (MQ=255)
                        cttcAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTGAATCATTGCTTCa                                                                                                                               >  2:1414311/1‑150 (MQ=255)
                           cAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAAcc                                                                                                                            >  2:582555/1‑150 (MQ=255)
                             aaCGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCGCGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCATTCATGGCTTCACCCGt                                                                                                                          >  1:137493/1‑150 (MQ=255)
                             aaCGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGt                                                                                                                          >  1:1081925/1‑150 (MQ=255)
                             aaCGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGt                                                                                                                          >  2:786137/1‑150 (MQ=255)
                             aaCGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGt                                                                                                                          >  1:935019/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)
                                          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)
                                                      gggAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTcc                                                                                                 >  2:403687/1‑150 (MQ=255)
                                                         aGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGc                                                                                              <  1:654577/150‑1 (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                                                                                                                                    <  2:543582/99‑1 (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)
                                                                          aTCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGcc                                                                             <  2:1081925/150‑1 (MQ=255)
                                                                          aTCTCTTCAGTAATGAATCGCCTCCCGTGGGATGGCTACCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGcc                                                                             <  2:137493/150‑1 (MQ=255)
                                                                            ctctTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAGGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGa                                                                           >  1:677672/1‑150 (MQ=255)
                                                                            ctctTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGa                                                                           >  2:491603/1‑150 (MQ=255)
                                                                               ttCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCATCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCa                                                                        <  2:1512976/150‑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)
                                                                                               cTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGAGACCGGGATCACGCCGACCAGGGTCAGCTTTTTCgg                                                        >  1:1540067/1‑150 (MQ=255)
                                                                                                 cccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCCAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAACCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGCCCAGGGTCAGCTTTTTCGGaa                                                      >  1:1286473/1‑150 (MQ=255)
                                                                                                 cccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGaa                                                      >  1:582383/1‑150 (MQ=255)
                                                                                                  ccGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGaaa                                                     <  2:71020/150‑1 (MQ=255)
                                                                                                   cGTGGGATGGCTTCCACACCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTTCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAAc                                                    >  1:1170823/1‑150 (MQ=255)
                                                                                                            gCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTCAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTg                                           >  2:358553/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)
                                                                                                                   cGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCa                                    <  2:1557129/150‑1 (MQ=255)
                                                                                                                   cGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCa                                    >  1:1389650/1‑150 (MQ=255)
                                                                                                                   cGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGGGAAGCGCa                                    >  1:1198489/1‑150 (MQ=255)
                                                                                                                    gCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAg                                   >  2:1119287/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)
                                                                                                                                  gcgcagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCACCCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGt                     <  2:1559834/150‑1 (MQ=255)
                                                                                                                                     cagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCgg                  >  1:1475366/1‑150 (MQ=255)
                                                                                                                                     cagcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTCAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCgg                  <  1:358553/150‑1 (MQ=255)
                                                                                                                                      agcCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggc                 <  1:403687/150‑1 (MQ=255)
                                                                                                                                        ccAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggcca               >  1:153786/1‑150 (MQ=255)
                                                                                                                                        ccAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggcca               >  1:797780/1‑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:1668785/1‑150 (MQ=255)
                                                                                                                                            aaCCTGCGCTTGAGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCggccaggcc           >  1:1406898/1‑150 (MQ=255)
                                                                                                                                              ccTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCaa         >  2:197688/1‑150 (MQ=255)
                                                                                                                                                tGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGc       <  1:746543/150‑1 (MQ=255)
                                                                                                                                                tGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGc       <  2:804981/150‑1 (MQ=255)
                                                                                                                                                    cAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGAGCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGAt   <  1:353378/150‑1 (MQ=255)
                                                                                                                                                     aaGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGATg  >  1:1038586/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTGCATCGAACGCCGGGCAATTTCTTCAAACGCTGCCTGTACTTGCGCCTTCGGGGAGGTCTGGAAACCTGCTATCTCTTCAGTCATGAATCGCCTCCCGTGGGATGGCTTCCACGCCAGATTCACGCAGCGCAGCCAGAACCTGCTCAAGCGCAGGTTCAATCATTGCTTCAACCGTCGGCGTTAAGCCGATGTGTGGCTCCAGCGATTCCGGGATCACGCCGACCAGGGTCAGCTTTTTCGGAAACTCGCCGGTGAAGCGCAGGGCCGACAAGACGTCGGCCAGGCCAAGCTGATG  >  NC_000913/3140712‑3141010

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