Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,852,048 C→T 5.5% A62V (GCC→GTC)  hypC → hydrogenase maturation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,852,0480CT 5.5% 233.6 / 3.5 73A62V (GCC→GTC) hypChydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base C (37/31);  new base T (2/2);  total (39/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.54e-01

GACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTG  >  NC_000913/2851900‑2852195
                                                                                                                                                    |                                                                                                                                                   
gACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGTCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGc                                                                                                                                                     <  2:799735/149‑1 (MQ=255)
gACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGc                                                                                                                                                     <  4:131245/149‑1 (MQ=255)
      aaCCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGTCGAAGc                                                                                                                                               >  3:362130/1‑149 (MQ=255)
      aaCCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGTCGAAGc                                                                                                                                               >  1:422846/1‑149 (MQ=255)
         cAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCAcg                                                                                                                                            >  3:701080/1‑149 (MQ=255)
            gCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGa                                                                                                                                         <  1:641932/149‑1 (MQ=255)
            gCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGa                                                                                                                                         <  3:61861/149‑1 (MQ=255)
                      aCGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGAc                                                                                                                               >  3:366192/1‑149 (MQ=255)
                        gTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGc                                                                                                                             <  4:716017/149‑1 (MQ=255)
                                cATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACaaa                                                                                                                     <  1:457076/149‑1 (MQ=255)
                                 aTTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACaaaa                                                                                                                    >  1:365071/1‑149 (MQ=255)
                                        gcgATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGttt                                                                                                             <  4:366192/149‑1 (MQ=255)
                                          gATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGa                                                                                                           <  1:263666/149‑1 (MQ=255)
                                             gTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACa                                                                                                                  >  1:529556/1‑139 (MQ=255)
                                             gTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACa                                                                                                                  <  2:529556/139‑1 (MQ=255)
                                                      aCGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGa                                                                                               <  1:872791/149‑1 (MQ=255)
                                                       cGTTAGTCGGCTGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGAt                                                                                              <  3:599564/149‑1 (MQ=255)
                                                        gTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCgg                                                                                                <  2:807731/146‑1 (MQ=255)
                                                        gTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCgg                                                                                                >  1:807731/1‑146 (MQ=255)
                                                        gTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATg                                                                                             >  3:700505/1‑149 (MQ=255)
                                                          tggTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTc                                                                                           >  4:627244/3‑149 (MQ=255)
                                                                    cTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCtgt                                                                                 >  3:484679/1‑149 (MQ=255)
                                                                      gCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCtgttg                                                                               >  4:324073/1‑149 (MQ=255)
                                                                       cGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGTTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCtgttgt                                                                              <  2:487049/149‑1 (MQ=255)
                                                                        gATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGTCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTa                                                                             <  4:362130/149‑1 (MQ=255)
                                                                        gATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGTCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTa                                                                             <  2:422846/149‑1 (MQ=255)
                                                                         aTGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGttt                                                                                                             <  2:22417/116‑1 (MQ=255)
                                                                         aTGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGttt                                                                                                             >  1:22417/1‑116 (MQ=255)
                                                                            aaaaCGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGc                                                                         <  2:176160/149‑1 (MQ=255)
                                                                             aaaCGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCg                                                                        >  1:136789/1‑149 (MQ=255)
                                                                              aaCGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGa                                                                       >  1:220394/1‑149 (MQ=255)
                                                                              aaCGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGa                                                                       <  4:596362/149‑1 (MQ=255)
                                                                                   tCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGaaa                                                                  >  3:420828/1‑149 (MQ=255)
                                                                                       ccGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAaata                                                              >  2:580117/1‑149 (MQ=255)
                                                                                        cgcgcgTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAaataa                                                             <  1:388846/149‑1 (MQ=255)
                                                                                         gcgcgTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAaataat                                                            <  4:420828/149‑1 (MQ=255)
                                                                                         gcgcgTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAaataat                                                            >  1:255067/1‑149 (MQ=255)
                                                                                             gTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGt                                                        <  4:459150/149‑1 (MQ=255)
                                                                                              tGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGtt                                                       >  3:7561/1‑149 (MQ=255)
                                                                                              tGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGtt                                                       <  1:295152/149‑1 (MQ=255)
                                                                                                ggCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGtttt                                                     >  4:657925/1‑149 (MQ=255)
                                                                                                     gTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTtga                                                <  2:229940/149‑1 (MQ=255)
                                                                                                      tGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTtgat                                               <  3:557630/149‑1 (MQ=255)
                                                                                                          tACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGtt                                                                                                       <  4:494710/89‑1 (MQ=255)
                                                                                                          tACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGtt                                                                                                       >  3:494710/1‑89 (MQ=255)
                                                                                                                 acacGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAg                                                                      >  4:345181/1‑115 (MQ=255)
                                                                                                                 acacGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAg                                                                      <  3:345181/115‑1 (MQ=255)
                                                                                                                 acacGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATcgcgc                                    >  4:899922/1‑149 (MQ=255)
                                                                                                                     gTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCg                                >  4:8100/1‑149 (MQ=255)
                                                                                                                     gTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCg                                >  4:237590/1‑149 (MQ=255)
                                                                                                                       tGGCTTTGCCATGAGCGTAATTAATGAAGCCGCAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGCCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCgca                              >  2:577436/1‑147 (MQ=255)
                                                                                                                          cTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACa                           >  1:446078/1‑149 (MQ=255)
                                                                                                                           tttGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTtgat                                               <  4:855986/128‑1 (MQ=255)
                                                                                                                           tttGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTtgat                                               >  3:855986/1‑128 (MQ=255)
                                                                                                                           tttGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAg                          >  3:646747/1‑149 (MQ=255)
                                                                                                                           tttGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAg                          >  2:258225/1‑149 (MQ=255)
                                                                                                                           tgtGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAg                          >  1:216948/3‑149 (MQ=255)
                                                                                                                               ccATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGa                      >  1:852988/1‑149 (MQ=255)
                                                                                                                                 aTGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATg                    >  2:91617/1‑149 (MQ=255)
                                                                                                                                 aTGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATg                    <  4:646747/149‑1 (MQ=255)
                                                                                                                                  tGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGc                   >  1:758655/1‑149 (MQ=255)
                                                                                                                                   gAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCa                  <  3:324073/149‑1 (MQ=255)
                                                                                                                                       gTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTTGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGACGAATATCGCGCGCCGGAACAGGTGATGCAGTTa              >  4:610738/1‑149 (MQ=255)
                                                                                                                                        taattaatGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTaa             <  4:624459/149‑1 (MQ=255)
                                                                                                                                        taattaatGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTaa             <  4:484679/149‑1 (MQ=255)
                                                                                                                                         aattaatGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAAt            >  4:314917/1‑149 (MQ=255)
                                                                                                                                          attaatGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAAtt           >  1:335977/1‑149 (MQ=255)
                                                                                                                                           ttaatGAAGCCGCAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTg          >  4:782834/1‑149 (MQ=255)
                                                                                                                                           ttaatGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTg          <  1:767003/149‑1 (MQ=255)
                                                                                                                                            taatGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGa         <  2:255067/149‑1 (MQ=255)
                                                                                                                                            taatGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGa         <  1:578649/149‑1 (MQ=255)
                                                                                                                                            taatGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGa         <  1:365541/149‑1 (MQ=255)
                                                                                                                                            taatGAAGACGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGa         <  2:766747/149‑1 (MQ=255)
                                                                                                                                              atGAAGCCGAAGCACGCGACACTCACGACGCCTTACAAAACATGTTTGACGGTGAGCCGGATGTCGGCGCGCAGCTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATGGCTCGCCGGAACAGGGGATGAAGTTAATTGAGc       >  4:676734/1‑149 (MQ=255)
                                                                                                                                                   gCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATACAGTTAATTGAGCATCTg  >  2:304373/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTCGACGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTG  >  NC_000913/2851900‑2852195

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