Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,835,766 A→C 5.7% C554G (TGC→GGC)  hypF ← carbamoyl phosphate phosphatase and [NiFe] hydrogenase maturation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,835,7660AC 5.7% 282.2 / 4.3 93C554G (TGC→GGC) hypFcarbamoyl phosphate phosphatase and [NiFe] hydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base A (34/53);  new base C (2/4);  total (36/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGA  >  NC_000913/2835621‑2835915
                                                                                                                                                 |                                                                                                                                                     
gccACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTgc                                                                                                                                                   <  1:1063772/150‑1 (MQ=255)
gccACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTgc                                                                                                                                                   <  2:1374871/150‑1 (MQ=255)
  cACTGCAGCGAAACAACGCCCACACGATGACTCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTgcgc                                                                                                                                                 <  1:356269/150‑1 (MQ=255)
     tGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCcag                                                                                                                                              >  2:1265557/1‑150 (MQ=255)
      gCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGTGCCTgcgctttg                                                                                                                                             <  1:17987/150‑5 (MQ=255)
       cATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGTGCCTgcgctttgt                                                                                                                                            <  2:347916/150‑6 (MQ=255)
       cATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCcagca                                                                                                                                            >  2:1593493/1‑150 (MQ=255)
       cATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCcagca                                                                                                                                            <  1:1229580/150‑1 (MQ=255)
       cATCGAAAAAACGCCCACACGATGACACCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGTGCCTgcgctttgt                                                                                                                                            <  1:1273322/150‑6 (MQ=255)
        aTCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGTGCCTgcgctttgtg                                                                                                                                           <  2:1056112/150‑7 (MQ=255)
           gAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGtt                                                                                                                                        >  2:1362823/1‑150 (MQ=255)
              aaaaCGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCg                                                                                                                                     <  2:1388379/150‑1 (MQ=255)
               aaaCGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGc                                                                                                                                    <  1:448777/150‑1 (MQ=255)
                   gCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACg                                                                                                                                <  2:1135582/150‑1 (MQ=255)
                   gCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACg                                                                                                                                <  1:646057/150‑1 (MQ=255)
                      cacaCGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGgct                                                                                                                             >  1:593384/1‑150 (MQ=255)
                      cacaCGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTGTCCGCGTAATTCTGCCACTCCGGCTCAACGCGCAGGCACTGCGCCAGCAGGTTTCGCCCCGgct                                                                                                                             >  1:1101601/1‑150 (MQ=255)
                        cacGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGgctgc                                                                                                                           <  2:1444629/150‑1 (MQ=255)
                         acGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCTGGCACTGCGCCAGCAGGTTTCGCCACGgctgct                                                                                                                          <  1:1642483/150‑1 (MQ=255)
                         acGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGgctgct                                                                                                                          <  2:1181895/150‑1 (MQ=255)
                           gATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTc                                                                                                                        <  1:376658/150‑1 (MQ=255)
                           gATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTc                                                                                                                        >  2:817842/1‑150 (MQ=255)
                            aTGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCg                                                                                                                       >  1:1018076/1‑150 (MQ=255)
                            aTGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCg                                                                                                                       <  1:1709735/150‑1 (MQ=255)
                            aTGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCg                                                                                                                       >  1:1024434/1‑150 (MQ=255)
                            aTGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGACACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCg                                                                                                                       <  1:505884/150‑1 (MQ=255)
                               aCGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCAGTTTCAGAGTAATTCTGCCACTCCGGCACACAGCGCAGACACTGCGCCAGAAGTTTTCGCCACCGCTGCTTCGCTg                                                                                                                    >  1:1226250/1‑150 (MQ=255)
                                 gCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGcc                                                                                                                  <  2:1012932/150‑1 (MQ=255)
                                  ccAGCGGCGCGTTAATTCCAGGCTCAATGGCCCGCGCCAGCACTCTCCAGTTTTTCTGTAGCTCACTTGCTGTTTCAGAGTAATTCAGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCa                                                                                                                 <  1:449009/150‑1 (MQ=255)
                                  ccAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCa                                                                                                                 <  1:1716775/150‑1 (MQ=255)
                                  ccAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCa                                                                                                                 >  2:1568114/1‑150 (MQ=255)
                                  ccAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTATCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCa                                                                                                                 <  1:765788/150‑1 (MQ=255)
                                   cAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCaa                                                                                                                <  2:1362307/150‑1 (MQ=255)
                                     gcggcgCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAAt                                                                                                              >  2:175929/1‑150 (MQ=255)
                                       ggcgCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCg                                                                                                            <  1:1370173/150‑1 (MQ=255)
                                       ggcgCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCg                                                                                                            <  1:143587/150‑1 (MQ=255)
                                       ggcgCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCg                                                                                                            >  1:248173/1‑150 (MQ=255)
                                         cgcgTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGcc                                                                                                          <  1:1061044/150‑1 (MQ=255)
                                                  tCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAg                                                                                                 >  2:972240/1‑150 (MQ=255)
                                                    cACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACCCTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCCCTGCGCCCGCAGGTTTCGCCACGGCTGCTTCGCTGCCAACTCGCCACCCGGAAgcg                                                                                               >  2:299577/1‑150 (MQ=255)
                                                          cAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTg                                                                                         >  1:127218/1‑150 (MQ=255)
                                                              ggCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAgg                                                                                     <  2:912064/150‑1 (MQ=255)
                                                               gCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGc                                                                                    <  2:1133141/150‑1 (MQ=255)
                                                                cccGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCa                                                                                   >  2:1592558/1‑150 (MQ=255)
                                                                  cgcgCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTTCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAg                                                                                 <  1:842779/150‑1 (MQ=255)
                                                                  cgcgCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAg                                                                                 <  2:1312564/150‑1 (MQ=255)
                                                                  cgcgCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAg                                                                                 <  1:1593493/150‑1 (MQ=255)
                                                                   gcgcCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAgc                                                                                <  2:1712586/150‑1 (MQ=255)
                                                                       cAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAgccgcc                                                                            <  2:1022066/150‑1 (MQ=255)
                                                                          cACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAg                                                                         <  1:158586/150‑1 (MQ=255)
                                                                          cACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAg                                                                         <  1:171181/150‑1 (MQ=255)
                                                                           aCGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAgg                                                                        <  1:418167/150‑1 (MQ=255)
                                                                               tCCAGTTTTGCTGTTGCACACTTGCTGTATCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCt                                                                    <  1:1547187/150‑1 (MQ=255)
                                                                                  aGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTcac                                                                 >  1:1550983/1‑150 (MQ=255)
                                                                                  aGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTcac                                                                 <  1:1362823/150‑1 (MQ=255)
                                                                                     tttGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACAtt                                                              <  1:1568114/150‑1 (MQ=255)
                                                                                        gCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTcgc                                                           <  2:584115/150‑1 (MQ=255)
                                                                                         cTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCCGGCAAGCCGCCCAGGTGCTCACATTcgcg                                                          >  2:1246037/1‑150 (MQ=255)
                                                                                         cTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTcgcg                                                          <  1:350570/150‑1 (MQ=255)
                                                                                           gTTGCACACTTGCTGTTTCAGAGTAATTCTGCCCCTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCACATCGCCACCCGGAAGCGCCACTGCCGGCAAGCCCCCCAGGTGCTCACCTTCGCGag                                                        >  2:1141787/1‑149 (MQ=255)
                                                                                           gTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGAt                                                        <  1:1592558/150‑1 (MQ=255)
                                                                                             tGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCCCTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCCCTGCAGGCCAGCCGCCCAGGTGCTCACATTCGCGATAg                                                      >  2:850203/1‑150 (MQ=255)
                                                                                               cacaCTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGtt                                                    >  2:928489/1‑150 (MQ=255)
                                                                                                  acTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACACAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCAc                                                 <  2:973322/150‑1 (MQ=255)
                                                                                                  acTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCAc                                                 <  2:979919/150‑1 (MQ=255)
                                                                                                   cTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCTCCACTGCAGGCATGCCGCCCAGGTGCTCACATTCGCGATAGTTCAcc                                                <  2:826475/150‑1 (MQ=255)
                                                                                                   cTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCAcc                                                <  1:1265557/150‑1 (MQ=255)
                                                                                                       cagTTTCAGAGTAATTCTGCCAATCCGGCACAAAGCCCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCAAGGTGCTCACATTCGCGATAGTTCAACCGCa                                            >  1:1337642/3‑150 (MQ=255)
                                                                                                       cTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGTGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCa                                            <  2:834122/150‑1 (MQ=255)
                                                                                                       cTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCAAGCAGGTTTCGCCACGGCTGATTCGCTGCCAAAACGCCACCCGGCAGCGCCACTGCCGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGcc                                            >  1:665970/1‑149 (MQ=255)
                                                                                                         gTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCCCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATACGCCATAGTAACCCCGCAgg                                          >  1:825094/1‑150 (MQ=255)
                                                                                                            tCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCAc                                       <  1:832654/150‑1 (MQ=255)
                                                                                                               gagTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTcg                                    <  1:768329/150‑1 (MQ=255)
                                                                                                               gagTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTcc                                    >  2:944824/1‑149 (MQ=255)
                                                                                                                     ttCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGACGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCCCCCGCAGGCACTCGCCGccc                              >  1:580497/1‑150 (MQ=255)
                                                                                                                       cTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCAAATTCGCGAGAGTGCACCCGCAGGCACTCGCCGccccc                            >  1:1531380/1‑149 (MQ=255)
                                                                                                                         gCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCcaca                          >  1:1255833/1‑150 (MQ=255)
                                                                                                                         gCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCcaca                          >  1:228656/1‑150 (MQ=255)
                                                                                                                           cACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACaaa                        <  1:1245014/150‑1 (MQ=255)
                                                                                                                           cACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACaaa                        >  1:866313/1‑150 (MQ=255)
                                                                                                                             cTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCCCCCGCAGGCACTCGCCGCCCCACAAAgc                      >  1:664862/1‑150 (MQ=255)
                                                                                                                              tCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAgcg                     <  2:248173/150‑1 (MQ=255)
                                                                                                                               ccGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAgcgc                    >  2:1045795/1‑150 (MQ=255)
                                                                                                                               ccGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAgcgc                    >  1:867558/1‑150 (MQ=255)
                                                                                                                                     caAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGtgcg              >  1:969373/1‑147 (MQ=255)
                                                                                                                                      aaaGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTtctc             >  1:1281444/1‑150 (MQ=255)
                                                                                                                                        aGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTccc           <  2:983975/150‑1 (MQ=255)
                                                                                                                                        aGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTccc           >  1:1019242/1‑150 (MQ=255)
                                                                                                                                          cgcAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTccccc         <  2:78268/150‑1 (MQ=255)
                                                                                                                                            cAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCat       <  1:1106256/150‑1 (MQ=255)
                                                                                                                                             aGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCata      <  1:398741/150‑1 (MQ=255)
                                                                                                                                               gCACTGCGCCAGCAGGGTTCGCCACGGCGGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCatacc    <  2:1101601/150‑1 (MQ=255)
                                                                                                                                                cACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCataccg   <  2:1550983/150‑1 (MQ=255)
                                                                                                                                                 aCTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCCCAAAGCGCCGTTCTCCCCCataccgc  >  2:24290/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCCACTGCATCGAAAAAACGCCCACACGATGACGCCAGCGGCGCGTTAATTCCACGCTCAATGGCCCGCGCCAGCACGCTCCAGTTTTGCTGTTGCACACTTGCTGTTTCAGAGTAATTCTGCCACTCCGGCACAAAGCGCAGGCACTGCGCCAGCAGGTTTCGCCACGGCTGCTTCGCTGCCAAATCGCCACCCGGAAGCGCCACTGCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGA  >  NC_000913/2835621‑2835915

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