Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,854,312 T→A 5.3% intergenic (+48/‑26) hypE → / → fhlA carbamoyl dehydratase, hydrogenases 1,2,3 maturation protein/formate hydrogenlyase transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,854,3120TA 5.3% 451.0 / 11.4 151intergenic (+48/‑26)hypE/fhlAcarbamoyl dehydratase, hydrogenases 1,2,3 maturation protein/formate hydrogenlyase transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (76/67);  new base A (4/4);  total (80/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01

Only 100 of 152 total aligned reads displayed.

GTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCG  >  NC_000913/2854164‑2854456
                                                                                                                                                    |                                                                                                                                                
        gtggAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCt                                                                                                                                         >  1:1800410/1‑150 (MQ=255)
          ggAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGTTTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCttt                                                                                                                                       >  2:1163063/1‑150 (MQ=255)
          ggAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCCCACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCttt                                                                                                                                       >  2:1680403/1‑150 (MQ=255)
          ggAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCttt                                                                                                                                       >  1:800548/1‑150 (MQ=255)
           gAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCAAAGCGGTGcaaaa                                                                                                                                      <  1:1607236/150‑5 (MQ=255)
           gAACGTAAAGGTGTTCGTCTTGACGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCAAAGCGGTGcaaaa                                                                                                                                      <  1:164057/150‑5 (MQ=255)
            aaCGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCAAAGCGGTGcaaaag                                                                                                                                     <  1:1282586/150‑6 (MQ=255)
            aaCGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGAGTGAAACGAACCCTCGAGTTACCACACTCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATCGTTAGTCAATGACCTTTTGCACCGCAAAGCGGTGcaaaag                                                                                                                                     <  1:1502679/150‑6 (MQ=255)
               gTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTgg                                                                                                                                  >  2:1053238/1‑150 (MQ=255)
                 aaaGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGgaa                                                                                                                                <  1:1783002/150‑1 (MQ=255)
                  aaGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGgaag                                                                                                                               <  2:708396/150‑1 (MQ=255)
                   aGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGgaaga                                                                                                                              <  1:167637/150‑1 (MQ=255)
                    ggTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGgaagaa                                                                                                                             >  1:371707/1‑150 (MQ=255)
                      tgtTCGTCTTGCCGATCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACa                                                                                                                           >  1:1144057/1‑150 (MQ=255)
                           gTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATg                                                                                                                      <  1:198181/150‑1 (MQ=255)
                           gTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATg                                                                                                                      <  2:551696/150‑1 (MQ=255)
                             cTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTc                                                                                                                    >  1:804922/1‑150 (MQ=255)
                             cTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTc                                                                                                                    <  2:606897/150‑1 (MQ=255)
                              ttGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCa                                                                                                                   <  2:1326341/150‑1 (MQ=255)
                               tGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCat                                                                                                                  <  1:1419967/150‑1 (MQ=255)
                                gCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCata                                                                                                                 >  1:1253286/1‑150 (MQ=255)
                                gCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGGTTATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAAGGACCTTTTGCACCGCTTTGCGGTGGGTGCCTGGAAGAACAAAATGTCatg                                                                                                                 >  2:1221702/1‑149 (MQ=255)
                                 ccGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCatat                                                                                                                >  2:1871695/1‑150 (MQ=255)
                                   ggTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATac                                                                                                              >  1:1312169/1‑150 (MQ=255)
                                      cTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACAcc                                                                                                           <  2:865863/150‑1 (MQ=255)
                                      cTGTATGGCGTGAAACGAACACTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAATTTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTACTGGAAGAACAAAATGTCATATACAcc                                                                                                           <  2:432245/150‑1 (MQ=255)
                                          aTGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATg                                                                                                       >  2:1016045/1‑150 (MQ=255)
                                              cGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTg                                                                                                   >  2:1794774/1‑150 (MQ=255)
                                              cGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTg                                                                                                   <  2:821239/150‑1 (MQ=255)
                                              cGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTg                                                                                                   >  1:731665/1‑150 (MQ=255)
                                               gTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGa                                                                                                  <  2:1740936/150‑1 (MQ=255)
                                                  aaaCGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGAtct                                                                                               <  1:315601/150‑1 (MQ=255)
                                                   aaCGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGAtctc                                                                                              >  1:1011271/1‑150 (MQ=255)
                                                    aCGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCg                                                                                             >  1:1482445/1‑150 (MQ=255)
                                                     cGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCgg                                                                                            <  2:648665/150‑1 (MQ=255)
                                                     cGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCgg                                                                                            <  1:1616247/150‑1 (MQ=255)
                                                       aaCCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGac                                                                                          <  1:683987/150‑1 (MQ=255)
                                                       aaCCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGac                                                                                          <  2:554078/150‑1 (MQ=255)
                                                       aaCCCTCGATTTACAACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACGGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGac                                                                                          <  2:1166339/150‑1 (MQ=255)
                                                          ccTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGacaac                                                                                       >  1:1354788/1‑150 (MQ=255)
                                                           cTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGacaaca                                                                                      <  1:863902/150‑1 (MQ=255)
                                                              gATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAgg                                                                                   <  2:671029/150‑1 (MQ=255)
                                                                   aCCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGgttgt                                                                              <  2:1273332/150‑1 (MQ=255)
                                                                     cacaCGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTc                                                                            >  2:1139451/1‑150 (MQ=255)
                                                                      acacGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTCTAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCg                                                                           <  1:317382/150‑1 (MQ=255)
                                                                      acacGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCg                                                                           <  2:895094/150‑1 (MQ=255)
                                                                      acacGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCg                                                                           >  2:385520/1‑150 (MQ=255)
                                                                        acGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGAc                                                                         <  2:143856/150‑1 (MQ=255)
                                                                         cGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACa                                                                        >  1:516595/1‑150 (MQ=255)
                                                                             gacccGCTTGCTCGTATATGCTAATATAATTCTAACTCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTCCGGTGCTTTCCTGGAATAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCAc                                                                    <  2:76234/147‑1 (MQ=255)
                                                                               aCCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTc                                                                  >  1:141945/1‑150 (MQ=255)
                                                                                  gCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGa                                                               >  2:1239093/1‑150 (MQ=255)
                                                                                    ttCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATGTCGAACAACCCTTGTTGTCCGAGATCACTcatcgg                                                             <  1:461888/150‑6 (MQ=37)
                                                                                    ttCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGaca                                                             <  1:549400/150‑1 (MQ=255)
                                                                                     tCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATGTCGAACAACCCTTGTTGTCCGAGATCACTcatcggt                                                            <  1:281091/150‑7 (MQ=37)
                                                                                     tCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGacac                                                            >  1:1755136/1‑150 (MQ=255)
                                                                                      ccGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACt                                                           >  1:1386335/1‑150 (MQ=255)
                                                                                        gcgTATATTCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTAt                                                         <  2:1242272/150‑1 (MQ=255)
                                                                                        gcgTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTAt                                                         >  1:928426/1‑150 (MQ=255)
                                                                                           tataTGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTgc                                                      >  2:1761509/1‑150 (MQ=255)
                                                                                            atatGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTgca                                                     <  1:836808/150‑1 (MQ=255)
                                                                                             tatGCTAATAAAATTCTAAATCTCCAATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTgcag                                                    <  2:262613/150‑1 (MQ=255)
                                                                                               tGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGTACACTATTgcagca                                                  >  1:945084/1‑150 (MQ=255)
                                                                                                gCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCCCCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATCTACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGcagcag                                                 <  1:372495/150‑1 (MQ=255)
                                                                                                    ataaAATTCTAAATCTACTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCAGATCCACCGATGACTGATCTCGGACAACTACGGTGGTTCGACAGCCGTCGGACACTATTGAagaagcagg                                             >  2:310821/1‑147 (MQ=255)
                                                                                                     taaAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGa                                            >  2:880651/1‑150 (MQ=255)
                                                                                                     taaAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGa                                            >  2:1815879/1‑150 (MQ=255)
                                                                                                       aaaTTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATc                                          <  2:175507/150‑1 (MQ=255)
                                                                                                       aaaTTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATc                                          >  2:31387/1‑150 (MQ=255)
                                                                                                        aaTTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCt                                         >  2:193864/1‑150 (MQ=255)
                                                                                                        aaTTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAACGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCt                                         <  2:101734/150‑1 (MQ=255)
                                                                                                           tCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGc                                      <  1:1163063/150‑1 (MQ=255)
                                                                                                             tAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCt                                    >  2:104539/1‑150 (MQ=255)
                                                                                                               aaTCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTACTGGAAGAACAAAAAGTCATATACACCGATCAGTGATCTCGTACATCAAGGGTTGTTCGACAACACTCGGAAACAATTACATCAGCCCGATCTGTCATCg                                  >  1:1364434/1‑150 (MQ=255)
                                                                                                               aaTCTCCTAAAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCTGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCg                                  <  2:350075/150‑1 (MQ=255)
                                                                                                                 tctcCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCt                                <  1:1397158/150‑1 (MQ=255)
                                                                                                                 tctcCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCt                                <  2:764636/150‑1 (MQ=255)
                                                                                                                  ctcCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCtg                               >  2:805769/1‑150 (MQ=255)
                                                                                                                    ccTATAGTTAGTCATTGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGCGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCAGATCTGGCCTCGCtgtg                             >  2:715043/1‑150 (MQ=255)
                                                                                                                    ccTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAATAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGTGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCtggg                             >  1:1008842/1‑148 (MQ=255)
                                                                                                                    ccTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGCCCTCGCtgtg                             >  1:830590/1‑150 (MQ=255)
                                                                                                                      tttagttagtCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCtgtgtg                           <  2:454444/148‑1 (MQ=255)
                                                                                                                      tttagtgagtCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGCACAAAATGTAATATACACCGATGCGTGATCTCGGACAACAAGGGTTGTCCGACATCACTCGGACACTAGTGCAGCAGCCCGATCTGGCCTCGCggtgtg                           <  2:263512/148‑1 (MQ=255)
                                                                                                                        tagttagtCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACACGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAg                         <  2:1431086/150‑1 (MQ=255)
                                                                                                                        tagttagtCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAg                         >  1:180321/1‑150 (MQ=255)
                                                                                                                          gttagtCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGc                       >  2:802900/1‑150 (MQ=255)
                                                                                                                          gttagtCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGc                       <  2:804922/150‑1 (MQ=255)
                                                                                                                          gttagtCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGc                       <  1:398229/150‑1 (MQ=255)
                                                                                                                                cAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGTCTCGTTc                 >  2:278337/1‑150 (MQ=255)
                                                                                                                                cAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTc                 <  2:1386335/150‑1 (MQ=255)
                                                                                                                                 aaTGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCg                <  1:883928/150‑1 (MQ=255)
                                                                                                                                   tGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCa              <  2:692541/150‑1 (MQ=255)
                                                                                                                                     caggaaagCACCGCAAAGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAAc            >  2:1617164/8‑150 (MQ=255)
                                                                                                                                     caggaaagCACCGCAAAGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGATCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAAc            >  2:819788/8‑150 (MQ=255)
                                                                                                                                      ccTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACt           >  2:1562763/1‑150 (MQ=255)
                                                                                                                                        ttttGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTgg         <  2:800548/150‑1 (MQ=255)
                                                                                                                                         aaagCACCGCAAAGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGt        >  1:1264523/4‑150 (MQ=255)
                                                                                                                                            gCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTaaa     >  2:1294369/1‑150 (MQ=255)
                                                                                                                                             cACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAg    >  1:1836941/1‑150 (MQ=255)
                                                                                                                                             cACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAg    <  2:1383148/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCACACGCCGAACCGCTTCCGCGTATATGCTAATAAAATTCTAAATCTCCTATAGTTAGTCAATGACCTTTTGCACCGCTTTGCGGTGCTTTCCTGGAAGAACAAAATGTCATATACACCGATGAGTGATCTCGGACAACAAGGGTTGTTCGACATCACTCGGACACTATTGCAGCAGCCCGATCTGGCCTCGCTGTGTGAGGCTCTTTCGCAACTGGTAAAGCG  >  NC_000913/2854164‑2854456

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: