Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,850,816 Δ1 bp 6.2% coding (170/351 nt) hypA → hydrogenase 3 nickel incorporation protein HypA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,850,8160C. 6.2% 210.7 / 7.5 64coding (170/351 nt)hypAhydrogenase 3 nickel incorporation protein HypA
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base . (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

CAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTA  >  NC_000913/2850668‑2850964
                                                                                                                                                    |                                                                                                                                                    
cAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCg                                                                                                                                                     <  2:491221/150‑1 (MQ=255)
 aaCGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCgg                                                                                                                                                    <  1:507617/150‑1 (MQ=255)
       cACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGtt                                                                                                                                              >  2:542576/1‑150 (MQ=255)
         cTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGt                                                                                                                                            >  1:520316/1‑150 (MQ=255)
           ggAATTGATCGATCAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTTGCTCAAAATTGGCGCATTTTCTTTTTTCGAAACCAGCTATCTTGCCTTTTGTTTTGATCTGGTTTGGCGAGGCAGCGTGGCGGAAGTTTGTaa                                                                                                                                          <  2:674854/150‑1 (MQ=255)
                tGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGG‑GGAAGGTTGTAAACTGCa                                                                                                                                    >  2:359440/1‑150 (MQ=255)
                tGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGG‑GGAAGGTTGTAAACTGCa                                                                                                                                    >  2:364562/1‑150 (MQ=255)
                 gATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCTCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCa                                                                                                                                    <  2:1263148/150‑1 (MQ=255)
                 gATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCa                                                                                                                                    <  1:583864/150‑1 (MQ=255)
                       aCAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCga                                                                                                                              <  2:797787/150‑1 (MQ=255)
                        cagcagGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCgaa                                                                                                                             <  1:927475/150‑1 (MQ=255)
                          gcagGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCgaaga                                                                                                                           >  1:1063742/1‑150 (MQ=255)
                          gcagGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCgaaga                                                                                                                           <  2:117992/150‑1 (MQ=255)
                          gcagGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCgaaga                                                                                                                           >  2:788664/1‑150 (MQ=255)
                          gcagGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTACCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCgaaga                                                                                                                           <  1:176636/150‑1 (MQ=255)
                             ggCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACa                                                                                                                        >  2:1075689/1‑150 (MQ=255)
                               ccGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAag                                                                                                                      >  2:24545/1‑150 (MQ=255)
                                 gCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAagag                                                                                                                    >  2:1204339/1‑150 (MQ=255)
                                  cAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAgg                                                                                                                   <  2:418539/150‑1 (MQ=255)
                                   aaaaCACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGc                                                                                                                  <  2:529640/150‑1 (MQ=255)
                                         cGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATg                                                                                                            <  1:1033549/150‑1 (MQ=255)
                                           gcgcAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGAGCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCt                                                                                                          >  2:489819/1‑150 (MQ=255)
                                               aaaaCGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGgtg                                                                                                      <  1:1177671/150‑1 (MQ=255)
                                                aaaCGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGgtgt                                                                                                     <  1:600030/150‑1 (MQ=255)
                                                 aaCGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGtgtg                                                                                                    <  1:524155/150‑1 (MQ=255)
                                                 aaCGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGtgtg                                                                                                    <  2:451361/150‑1 (MQ=255)
                                                     cgTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAAc                                                                                                >  1:273/1‑150 (MQ=255)
                                                             gggTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAAc                                                                                        <  2:670630/150‑1 (MQ=255)
                                                                   ggCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTAtg                                                                                  <  1:1075689/150‑1 (MQ=255)
                                                                     cTCAAAATTGTCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTTCAGCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTAtgtg                                                                                <  2:1164220/150‑1 (MQ=255)
                                                                     cTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTAtgtg                                                                                <  1:54894/150‑1 (MQ=255)
                                                                       cAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGAc                                                                              >  1:940074/1‑150 (MQ=255)
                                                                       cAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGAc                                                                              >  2:204239/1‑150 (MQ=255)
                                                                         aaaTTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGc                                                                            >  1:87694/1‑150 (MQ=255)
                                                                              ggCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTg                                                                       >  2:1160263/1‑150 (MQ=255)
                                                                              ggCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTg                                                                       >  2:1289686/1‑150 (MQ=255)
                                                                                 gcATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGAcc                                                                    <  1:1403446/150‑1 (MQ=255)
                                                                                      ttCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAgcg                                                               <  1:1135326/150‑1 (MQ=255)
                                                                                        cTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGcgcg                                                             >  1:1247147/1‑150 (MQ=255)
                                                                                           gtgtCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTcc                                                          >  1:525153/1‑150 (MQ=255)
                                                                                                        ctctctTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGG‑GGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAg                                            <  1:364562/150‑1 (MQ=255)
                                                                                                        ctctctTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGG‑GGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAg                                            <  1:359440/150‑1 (MQ=255)
                                                                                                          ctctTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAgt                                           <  1:1204339/150‑1 (MQ=255)
                                                                                                            ctTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGcgcg                                                             >  2:990058/1‑130 (MQ=255)
                                                                                                            ctTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGcgcg                                                             <  1:990058/130‑1 (MQ=255)
                                                                                                            ctTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAgtgt                                         >  2:865350/1‑150 (MQ=255)
                                                                                                                        tttgATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACAtg                             <  2:180606/150‑1 (MQ=255)
                                                                                                                         ttgATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACAtgc                            <  2:1063742/150‑1 (MQ=255)
                                                                                                                         ttgATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACAtgc                            <  2:525153/150‑1 (MQ=255)
                                                                                                                          tgATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACAtgct                           >  1:819872/1‑150 (MQ=255)
                                                                                                                            aTCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACAtgctgc                         <  2:520316/150‑1 (MQ=255)
                                                                                                                             tCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCa                        >  1:575141/1‑150 (MQ=255)
                                                                                                                              cTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAg                       >  1:505718/1‑150 (MQ=255)
                                                                                                                                ggTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGAt                     >  1:1365464/1‑150 (MQ=255)
                                                                                                                                ggTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGAt                     <  2:1247147/150‑1 (MQ=255)
                                                                                                                                    tGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGATACTGAACGAGCGCGTCCGACGGTGTCAAGAGTGTAATGGTGAGATGCTGAAGATtgtg                 <  2:940074/150‑1 (MQ=255)
                                                                                                                                     gCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTgg                >  1:795240/1‑150 (MQ=255)
                                                                                                                                     gCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTgg                >  2:1406243/1‑150 (MQ=255)
                                                                                                                                      ccGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGc               >  1:410414/1‑150 (MQ=255)
                                                                                                                                             aGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAgacgac        >  1:857354/1‑150 (MQ=255)
                                                                                                                                             aGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAgacgac        >  2:977049/1‑150 (MQ=255)
                                                                                                                                               cGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACgg      <  1:489819/150‑1 (MQ=255)
                                                                                                                                               cGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACgg      >  2:971463/1‑150 (MQ=255)
                                                                                                                                                   gCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTa  <  1:971463/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAACGGGCACTGGAATTGATCGAACAGCAGGCCGCAAAACACGGCGCAAAACGCGTAACTGGGGTCTGGCTCAAAATTGGCGCATTTTCTTGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTA  >  NC_000913/2850668‑2850964

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