Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,143,980 C→A 5.9% G62V (GGC→GTC)  hybB ← hydrogenase 2 membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,143,9800CA 5.9% 309.4 / 4.9 91G62V (GGC→GTC) hybBhydrogenase 2 membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (47/38);  new base A (2/4);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCCACGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTAC  >  NC_000913/3143831‑3144124
                                                                                                                                                     |                                                                                                                                                
gACCCACGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGc                                                                                                                                                  <  1:294567/150‑1 (MQ=255)
  cccACGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGccc                                                                                                                                                >  2:41970/1‑149 (MQ=255)
  cccACGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCa                                                                                                                                                >  1:551851/1‑150 (MQ=255)
    cACGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAAt                                                                                                                                              >  1:1238528/1‑150 (MQ=255)
      cGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATca                                                                                                                                            <  1:1214193/150‑1 (MQ=255)
      cGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATca                                                                                                                                            >  1:450511/1‑150 (MQ=255)
               tGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCaa                                                                                                                                   >  1:1196831/1‑150 (MQ=255)
                gATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCaaa                                                                                                                                  >  2:989851/1‑150 (MQ=255)
                gATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCaaa                                                                                                                                  >  2:448981/1‑150 (MQ=255)
                   cGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAAcgc                                                                                                                               >  1:845311/1‑150 (MQ=255)
                      cAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGAt                                                                                                                            <  1:498601/150‑1 (MQ=255)
                        aGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATcc                                                                                                                          <  1:688401/150‑1 (MQ=255)
                         gCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCca                                                                                                                         <  2:113133/150‑1 (MQ=255)
                           cACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCcaca                                                                                                                       <  1:397774/150‑1 (MQ=255)
                               cAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGcc                                                                                                                   >  1:75680/1‑150 (MQ=255)
                               cAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGcc                                                                                                                   <  1:1602582/150‑1 (MQ=255)
                                  tGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCa                                                                                                                >  2:1528094/1‑150 (MQ=255)
                                  tGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCa                                                                                                                >  2:1209192/1‑150 (MQ=255)
                                   gAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCAc                                                                                                               >  2:87176/1‑150 (MQ=255)
                                     gTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACgg                                                                                                             <  1:1144876/150‑1 (MQ=255)
                                         ccAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCTCGCAGCCGCCACAGGCAAAGCCGGTGCCATTCAGCATGTCAAACGCGATCCACACGCCCCACGGGAAg                                                                                                         <  2:458521/150‑1 (MQ=255)
                                          cAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATc                                                                                                                                             <  2:29698/113‑1 (MQ=255)
                                          cAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATc                                                                                                                                             >  1:29698/1‑113 (MQ=255)
                                          cAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAgc                                                                                                        >  1:1083516/1‑150 (MQ=255)
                                           aaaCAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCTGTGCAAATCAGCTGGTCAAACGAGATCCACACGCCCCACGGGAAgcc                                                                                                       <  2:263902/150‑1 (MQ=255)
                                            aaCAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAgccg                                                                                                      <  1:558586/150‑1 (MQ=255)
                                               aGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGccgccg                                                                                                   >  2:276238/1‑150 (MQ=255)
                                                gACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGt                                                                                                  <  2:845311/150‑1 (MQ=255)
                                                  cTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTc                                                                                                >  2:1261810/1‑150 (MQ=255)
                                                  cTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTc                                                                                                >  1:628530/1‑150 (MQ=255)
                                                  cTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTc                                                                                                >  1:1602574/1‑150 (MQ=255)
                                                  cTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTc                                                                                                <  2:313956/150‑1 (MQ=255)
                                                    cGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAg                                                                                              >  2:474453/1‑150 (MQ=255)
                                                    cGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAg                                                                                              >  2:303992/1‑150 (MQ=255)
                                                     gCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGACAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAgg                                                                                             >  2:272550/1‑150 (MQ=255)
                                                        aaCAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCa                                                                                          >  1:138435/1‑150 (MQ=255)
                                                         aCAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACCCGCCCCACGGGACGCCGCCGTTCAGGTCag                                                                                         >  1:541851/1‑150 (MQ=255)
                                                         aCAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCag                                                                                         <  2:918857/150‑1 (MQ=255)
                                                           aGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCagag                                                                                       <  1:31707/150‑1 (MQ=255)
                                                               ccGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACag                                                                                   <  1:1260111/150‑1 (MQ=255)
                                                                     gCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCTGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCa                                                                             >  2:15137/1‑150 (MQ=255)
                                                                      cACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAg                                                                            >  2:731768/1‑150 (MQ=255)
                                                                         cAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGAcc                                                                         <  1:1122130/150‑1 (MQ=255)
                                                                         cAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGAcc                                                                         <  1:1072719/150‑1 (MQ=255)
                                                                         cAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGACAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGAcc                                                                         <  1:272550/150‑1 (MQ=255)
                                                                          aGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGAGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCg                                                                        >  2:1573582/1‑150 (MQ=255)
                                                                          aGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGCCAAACCCGGTGCCAATCAGCAGGTCAAACGCAATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCg                                                                        >  2:96696/1‑150 (MQ=255)
                                                                           gCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGa                                                                       <  1:626378/150‑1 (MQ=255)
                                                                              gATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAaca                                                                    <  1:989851/150‑1 (MQ=255)
                                                                              gATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAaca                                                                    <  1:1261810/150‑1 (MQ=255)
                                                                                tGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACAcc                                                                  >  2:383066/1‑150 (MQ=255)
                                                                                 ggTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGTAATTTTGTTTAAATTTAGGTGGTCAAACGAGATCTACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCa                                                                 <  2:592882/150‑1 (MQ=255)
                                                                                 ggTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGc                                                                                                                    <  1:702851/99‑1 (MQ=255)
                                                                                 ggTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGc                                                                                                                    >  2:702851/1‑99 (MQ=255)
                                                                                 ggTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCa                                                                 >  2:778522/1‑150 (MQ=255)
                                                                                     ttGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACg                                                             <  1:383066/150‑1 (MQ=255)
                                                                                      tGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGc                                                            >  2:376398/1‑150 (MQ=255)
                                                                                       gCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCt                                                           <  1:1345417/150‑1 (MQ=255)
                                                                                        cccACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGCTGCAAATCAGCAGGTCAAACGCGATCCACCCGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCta                                                          >  2:1300128/1‑149 (MQ=255)
                                                                                          cACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGCTCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTca                                                        >  1:709688/1‑150 (MQ=255)
                                                                                          cACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTca                                                        <  2:1083516/150‑1 (MQ=255)
                                                                                          cACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTca                                                        <  1:1528094/150‑1 (MQ=255)
                                                                                            cGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGACAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTcaca                                                      >  2:35256/1‑150 (MQ=255)
                                                                                                tAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCa                                                  >  2:1192162/1‑150 (MQ=255)
                                                                                                  aaGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAATTTTTTTTAATTTTAGGTGGTCAAACGAGATCTACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAgg                                                <  2:628530/150‑1 (MQ=255)
                                                                                                  aaGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAgg                                                >  2:17843/1‑150 (MQ=255)
                                                                                                  aaGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAgg                                                >  1:1442448/1‑150 (MQ=255)
                                                                                                  aaGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAgg                                                <  2:75680/150‑1 (MQ=255)
                                                                                                   aGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGa                                               >  1:1520381/1‑150 (MQ=255)
                                                                                                   aGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGa                                               <  1:87176/150‑1 (MQ=255)
                                                                                                    gACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAg                                              >  1:928892/1‑150 (MQ=255)
                                                                                                    gACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAg                                              <  1:303992/150‑1 (MQ=255)
                                                                                                      cGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCa                                            >  2:1563604/1‑150 (MQ=255)
                                                                                                        tataCCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATa                                          >  1:1459389/1‑150 (MQ=255)
                                                                                                             cGCCCACGCCAGCGCCCAGCCGCCCCAGGCAAAGCCGGTGCCAATCAGCAGGTCAACCGCGATCCCCACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGCCCGAACACCAGCCGCTTCCCAATCAGGAGCAGACAGAt                                     >  1:1356287/1‑150 (MQ=255)
                                                                                                                  aCGCCAGCGCCCAGCCGCCACAGGCAATTTTTTTTAATTTTAGGTGGTCAAACGAGATCTACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGa                                <  2:13215/150‑1 (MQ=255)
                                                                                                                  aCGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGa                                <  1:276238/150‑1 (MQ=255)
                                                                                                                     ccAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGCCAGAGCCCAGACCGAACACCAGACGCTTCCCAATCAGGAGCATACAGATGACGATTa                             >  1:1260349/1‑150 (MQ=255)
                                                                                                                         cgcCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACgg                         >  2:1392770/1‑150 (MQ=255)
                                                                                                                         cgcCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACgg                         <  2:551851/150‑1 (MQ=255)
                                                                                                                           cccAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTc                       >  2:1045696/1‑150 (MQ=255)
                                                                                                                            ccAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTcc                      <  1:1625861/150‑1 (MQ=255)
                                                                                                                             cAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTGACAATCAGGAGCATACAGATGACGATTAACGGTCCa                     <  1:474453/150‑1 (MQ=255)
                                                                                                                                   ccACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATc               <  1:17843/150‑1 (MQ=255)
                                                                                                                                   ccACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATc               <  1:1209192/150‑1 (MQ=255)
                                                                                                                                    cacaGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCa              <  2:450511/150‑1 (MQ=255)
                                                                                                                                    cacaGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCa              >  1:566810/1‑150 (MQ=255)
                                                                                                                                    cacaGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCa              <  1:1192162/150‑1 (MQ=255)
                                                                                                                                      caGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATg            <  1:448981/150‑1 (MQ=255)
                                                                                                                                        ggCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACCGAGCCCAGACCGAACACCAGACGCTTCACCATCAGGCGCATACAGATGCCCATTAACGGACCAAAAAAAATGAc          >  1:712211/1‑150 (MQ=255)
                                                                                                                                        ggCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGAc          >  2:556318/1‑150 (MQ=255)
                                                                                                                                              gCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGttt    >  2:1347187/1‑150 (MQ=255)
                                                                                                                                                cGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTAc  >  1:236552/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GACCCACGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTAC  >  NC_000913/3143831‑3144124

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