Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 912,587 C→G 7.4% G410R (GGC→CGC)  hcp ← hybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913912,5870CG 7.4% 239.1 / 10.3 81G410R (GGC→CGC) hcphybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base C (40/35);  new base G (2/4);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.02e-01

GTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTG  >  NC_000913/912438‑912736
                                                                                                                                                     |                                                                                                                                                     
gTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCTCAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTcgcc                                                                                                                                                       <  1:237540/150‑1 (MQ=255)
     aTTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTg                                                                                                                                                  >  1:1424611/1‑150 (MQ=255)
     aTTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTg                                                                                                                                                  <  1:950471/150‑1 (MQ=255)
          ccTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACCGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCg                                                                                                                                             >  1:881453/1‑150 (MQ=255)
          ccTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCg                                                                                                                                             >  1:373338/1‑150 (MQ=255)
          ccTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGAGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCg                                                                                                                                             <  2:1645248/150‑1 (MQ=255)
           cTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCACGCGGCGAGCGccac                                                                                                                                            <  2:537656/150‑3 (MQ=255)
            tGCAGCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGGCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCACGCGGCGAGCGccact                                                                                                                                           <  2:716435/150‑4 (MQ=255)
             gCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGATGTAGTGGCGCTCGCCGCGTGCGCCGTcc                                                                                                                                          >  1:72577/1‑149 (MQ=255)
                tCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCACGCGGCGAGCGccactactt                                                                                                                                       <  2:706574/150‑8 (MQ=255)
                  tACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGcc                                                                                                                                     >  1:1640836/1‑150 (MQ=255)
                  tACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTAACACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACCGcc                                                                                                                                     >  2:85289/1‑150 (MQ=255)
                     cAGCCGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCCCAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGccacc                                                                                                                                  >  2:1330227/1‑150 (MQ=255)
                     cAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGccacc                                                                                                                                  >  2:772410/1‑150 (MQ=255)
                       caCGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCaa                                                                                                                                <  1:1480445/149‑1 (MQ=255)
                         cgcgGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGc                                                                                                                              <  2:891595/150‑1 (MQ=255)
                                gACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGa                                                                                                                       >  1:1884/1‑150 (MQ=255)
                                 aCCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGat                                                                                                                      <  1:1441369/150‑1 (MQ=255)
                                 aCCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGCCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAATGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGat                                                                                                                      >  1:440323/1‑150 (MQ=255)
                                   cTTCGATATCGCCAAACTCAAGTTTGATAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGACATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGatat                                                                                                                    >  1:411236/1‑150 (MQ=255)
                                    ttCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGGGGCGAAATCGGTGAAGTCGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATg                                                                                                                   >  2:253395/1‑150 (MQ=255)
                                     tCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGa                                                                                                                  >  1:909627/1‑150 (MQ=255)
                                       gATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACg                                                                                                                <  2:976399/150‑1 (MQ=255)
                                       gATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACg                                                                                                                >  2:1486572/1‑150 (MQ=255)
                                       gATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCACGCGGCGAGCGccactact                                                                                                                                        >  2:1442365/1‑120 (MQ=255)
                                       gATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCACGCGGCGAGCGccactact                                                                                                                                        <  1:1442365/126‑7 (MQ=255)
                                          atCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAgg                                                                                                                           <  2:1054318/136‑1 (MQ=255)
                                          atCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAgg                                                                                                                           >  1:1054318/1‑136 (MQ=255)
                                                aaaCTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCCGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTc                                                                                                       <  1:1324871/150‑1 (MQ=255)
                                                 aaCTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCa                                                                                                      <  1:222552/150‑1 (MQ=255)
                                                 aaCTCAAGTTTGTTAAAGCCATATTTACCACAGGCCAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGAGGCGAAATCGGTGAAGTAGAGGCGCTCCCCGCGTGCGCCGTCCCAGCCACCAAGCAGGAAGATATGACGCAGTGTTTCa                                                                                                      >  2:1611731/1‑150 (MQ=255)
                                                  aCTCAAGTTTGTTAAAGCGATATTTACCCCAGGCGAGGGGCCGGATCAGGCTGTCAGCCGGCACGCTGGGGGCGAAATCGGTGAAGTAGTGGCGCGCGCCGCGTTCGCCGTCACCTCCCCCAAGCAGGAAGATATGACGCAGTTTATCAc                                                                                                     <  1:281475/150‑1 (MQ=255)
                                                  aCTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCAc                                                                                                     <  1:36315/150‑1 (MQ=255)
                                                  aCTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCAc                                                                                                     >  1:493027/1‑150 (MQ=255)
                                                   cTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAAGCGGGGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACg                                                                                                    <  2:1208050/150‑1 (MQ=255)
                                                     cAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGc                                                                                                  >  1:501278/1‑150 (MQ=255)
                                                     cAAGTTTGTTAAAGCGATATTTACCACAGCCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGc                                                                                                  <  2:298724/150‑1 (MQ=255)
                                                           tGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCa                                                                                            >  2:7975/1‑150 (MQ=255)
                                                               aaaGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGatc                                                                                        <  2:208877/150‑1 (MQ=255)
                                                                 aGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCTGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGatcaa                                                                                      >  2:794569/1‑150 (MQ=255)
                                                                 aGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGatcaa                                                                                      >  2:778967/1‑150 (MQ=255)
                                                                    gATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGatcaatca                                                                                   <  1:1486572/150‑1 (MQ=255)
                                                                      tatTTACCACAGGCGAGGGTCATGATCAGGCAGTCATCCGGCACGCTGGTGGCGATATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGc                                                                                 <  1:1436472/150‑1 (MQ=255)
                                                                       atTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCg                                                                                >  1:1583536/1‑150 (MQ=255)
                                                                       atTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCg                                                                                <  2:80045/150‑1 (MQ=255)
                                                                         ttACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTa                                                                              >  2:945284/1‑150 (MQ=255)
                                                                            ccACAGGCGAGGGGCAGGATCAGGCTGTCATCCGGCACGCGGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATca                                                                           <  2:208934/150‑1 (MQ=255)
                                                                              acaGGCGAGGGTCAGGATCAGGTAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGCCTCACCAGAACAATCAGCGTATcagc                                                                         >  1:520456/1‑150 (MQ=255)
                                                                                aGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATcagcag                                                                       <  2:1276602/150‑1 (MQ=255)
                                                                                 ggCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCagcagc                                                                      >  2:1551569/1‑150 (MQ=255)
                                                                                 ggCGAGGGTCAGGATCAGGCAGTCATCCGCCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCg                                                                                <  1:441884/140‑1 (MQ=255)
                                                                                 ggCGAGGGTCAGGATCAGGCAGTCATCCGCCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCg                                                                                >  2:441884/1‑140 (MQ=255)
                                                                                  gCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAgcg                                                                     <  1:521311/150‑1 (MQ=255)
                                                                                      gggTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCaa                                                                 >  1:793106/1‑150 (MQ=255)
                                                                                        gTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCCCCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAagc                                                               >  2:1351801/1‑150 (MQ=255)
                                                                                        gTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGCTATGCCGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAagc                                                               >  1:582655/1‑150 (MQ=255)
                                                                                              aTCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTc                                                         <  1:1228205/150‑1 (MQ=255)
                                                                                                  ggCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGc                                                     >  2:1027/1‑150 (MQ=255)
                                                                                                    cAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGa                                                   <  1:1689127/150‑1 (MQ=255)
                                                                                                     aGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGAc                                                  <  1:586438/150‑1 (MQ=255)
                                                                                                     aGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGAc                                                  >  1:941611/1‑150 (MQ=255)
                                                                                                      gTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGAcc                                                 <  1:40356/150‑1 (MQ=255)
                                                                                                      gTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGAcc                                                 >  1:374247/1‑150 (MQ=255)
                                                                                                      gTCATCCGGCACGCTGGGGGCGAAATCGGTGAAGTAGTGGCTCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGAcc                                                 <  2:72577/150‑1 (MQ=255)
                                                                                                       tCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCTTATCAGCAGCGCCAAGCAGCGTCTGGCGACCa                                                >  1:1044280/1‑150 (MQ=255)
                                                                                                          tCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCaaaa                                             <  1:1167786/150‑1 (MQ=255)
                                                                                                            cGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCCCCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCCGCGTCTGGCGCCCAAAAcc                                           >  1:1533600/1‑150 (MQ=255)
                                                                                                               cACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCCAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGCTCAATCAGCGTATCAGCAGCGCCACGCAGCGTCTGGCGACCAAAACCCAc                                        >  1:1425072/1‑150 (MQ=255)
                                                                                                                aCGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACg                                       <  2:1583536/150‑1 (MQ=255)
                                                                                                                     ggtggCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAACCCCACGGTGAt                                  >  1:1657239/1‑150 (MQ=255)
                                                                                                                     ggtggCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGATCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGAt                                  <  1:89236/150‑1 (MQ=255)
                                                                                                                      gtggCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGCCGCAGTGTTTCACGGCTGACAAGATCAATAAGCGTATCAGCAGCGCCAAGCAGCGGCTGGCGACCAAACCCCACGGTtata                                 >  2:258659/1‑150 (MQ=255)
                                                                                                                      gtggCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATa                                 <  1:778967/150‑1 (MQ=255)
                                                                                                                           gAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTg                            <  2:493027/150‑1 (MQ=255)
                                                                                                                              aTCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGGGCCCAAAACCCACGGGGATAAGGTGCgg                         >  2:1039994/1‑150 (MQ=255)
                                                                                                                                      agtagtGGCGCTCGCCGCGTGCGCTGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGc                 <  1:794569/150‑1 (MQ=255)
                                                                                                                                      agtagtGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGc                 <  1:1303896/150‑1 (MQ=255)
                                                                                                                                       gtagtGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCCGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCCCGGTGCTAAGGTGCGGCATTTCGcg                >  1:313078/1‑149 (MQ=255)
                                                                                                                                        tagtGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTg               <  2:1226438/150‑1 (MQ=255)
                                                                                                                                        tagtGGCGCTCGCCGCGTGCGCCGGCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTg               <  2:346727/150‑1 (MQ=255)
                                                                                                                                         gacggCGCACGCGGCGAGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCTAGCAGCGTCTGGCGGCCAAAACCCACGGTCAGAAGGTGCGGAATTTCGCTGt              >  2:1478087/4‑150 (MQ=255)
                                                                                                                                            ggCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGAGGCAGATTTTCACGGCTCCCCAGAACAATCAGCGTCTCAGCAGCGCCAAGCCGCGTCTGGCGACCAAAACCCACGGTGATAAGGAGGGGAATTTCGCTGTACg           >  2:1268379/1‑150 (MQ=255)
                                                                                                                                               gcgcgccgcGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGCGCGGAATTTCGCTGTACGGGa        <  1:142484/147‑1 (MQ=255)
                                                                                                                                                  cgccgcGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACGGGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAAATTCGCTGTACGGGAAGc     >  1:563208/1‑150 (MQ=255)
                                                                                                                                                    ccgcGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGCCGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCCGCGCCACGCAGCGTCTGGCGACCACAACCGACGGGGATAAGGTGCGGAATTTCGCTGTACGGGAAGcc  >  1:1525300/1‑149 (MQ=255)
                                                                                                                                                     cgctTGCGCCGTCACAGGCTCCAAGGAGGAAGATATGACGGAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAAGCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTg  <  2:495114/150‑1 (MQ=255)
                                                                                                                                                     cgcgTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTg  <  2:411236/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTACATTGACCTGCATCTACCAGACGCGGCAGACCTTCGATATCGCCAAACTCAAGTTTGTTAAAGCGATATTTACCACAGGCGAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCCACCAAGCAGGAAGATATGACGCAGTTTTTCACGGCTCACCAGATCAATCAGCGTATCAGCAGCGCCAAGCAGCGTCTGGCGACCAAAACCCACGGTGATAAGGTGCGGAATTTCGCTGTACGGGAAGCCTG  >  NC_000913/912438‑912736

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: