Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,843,719 C→T E164E (GAG→GAA hycG ← hydrogenase 3 and formate hydrogenase complex, HycG subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,843,7190CT98.7% 269.6 / ‑3.8 75E164E (GAG→GAAhycGhydrogenase 3 and formate hydrogenase complex, HycG subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (42/31);  total (44/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

GGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTAGAGATCGTGGA  >  NC_000913/2843575‑2843846
                                                                                                                                                |                                                                                                                               
ggCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCcagc                                                                                                                           <  1:288142/151‑1 (MQ=255)
 gCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCcagca                                                                                                                          <  2:945835/151‑1 (MQ=255)
  cAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCcagcag                                                                                                                         <  2:284312/151‑1 (MQ=255)
  cAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCcagcag                                                                                                                         <  2:1052327/151‑1 (MQ=255)
  cAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCcagc                                                                                                                           >  1:777370/1‑149 (MQ=255)
   aaTCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAgg                                                                                                                        >  1:301815/1‑151 (MQ=255)
   aaTCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAgg                                                                                                                        <  1:1213254/151‑1 (MQ=255)
     tCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGcc                                                                                                                      >  1:77136/1‑151 (MQ=255)
        gacgacCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCg                                                                                                                     >  1:936420/1‑149 (MQ=255)
          cgacCGTAACGATAACCCGCCAGGCGACGTGCTTCGGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAgcgc                                                                                                                <  2:527755/151‑1 (MQ=255)
             ccGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGcc                                                                                                               >  1:302376/1‑150 (MQ=255)
             ccGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCa                                                                                                              <  2:156179/151‑1 (MQ=255)
             ccGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTCCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGGTCCAGCAGGCCGAGCGcc                                                                                                               >  1:221066/1‑150 (MQ=255)
                       aaCCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGc                                                                                                     <  2:1184487/150‑1 (MQ=255)
                         cccGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGcc                                                                                                    <  1:1054698/149‑1 (MQ=255)
                         cccGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCg                                                                                                   <  1:132528/150‑1 (MQ=255)
                           cGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTAc                                                                                                <  1:367939/151‑1 (MQ=255)
                            gCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTAc                                                                                                >  2:914227/1‑150 (MQ=255)
                              cAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGc                                                                                             >  1:153821/1‑151 (MQ=255)
                                 gCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTg                                                                                          >  2:857067/1‑151 (MQ=255)
                                  cGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTgg                                                                                         >  1:232113/1‑151 (MQ=255)
                                  cGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTgg                                                                                         <  2:1135680/151‑1 (MQ=255)
                                   gACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTggc                                                                                        >  2:627521/1‑151 (MQ=255)
                                       tGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTggcggc                                                                                     >  1:1104469/1‑150 (MQ=255)
                                       tGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTggcggc                                                                                     >  2:509025/1‑150 (MQ=255)
                                          ttCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGc                                                                                 >  2:267447/1‑151 (MQ=255)
                                          ttCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGc                                                                                 >  2:1148704/1‑151 (MQ=255)
                                            cgcgATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCg                                                                                >  2:347609/1‑150 (MQ=255)
                                            cgcgATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGt                                                                               >  2:855394/1‑151 (MQ=255)
                                             gcgATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGtt                                                                              <  2:730555/151‑1 (MQ=255)
                                                   cgcTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCgg                                                                         <  1:460423/148‑1 (MQ=255)
                                                   cACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCggg                                                                        >  2:561985/1‑151 (MQ=255)
                                                   cACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCggg                                                                        >  1:834234/1‑151 (MQ=255)
                                                         cACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGACATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCa                                                                  <  1:337362/151‑1 (MQ=255)
                                                          aCGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGcc                                                                   >  1:630495/1‑149 (MQ=255)
                                                           cgcgCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAgg                                                                <  2:545478/151‑1 (MQ=255)
                                                                aGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGatat                                                           <  1:565084/151‑1 (MQ=255)
                                                                aGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGatat                                                           <  2:478182/151‑1 (MQ=255)
                                                                aGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGata                                                            >  1:1101288/1‑150 (MQ=255)
                                                                  cGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGAtata                                                          >  1:693766/1‑150 (MQ=255)
                                                                    gCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATaaa                                                        <  1:732058/150‑1 (MQ=255)
                                                                       gCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACa                                                      <  1:639764/149‑1 (MQ=255)
                                                                       gCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATc                                                    >  1:170766/1‑151 (MQ=255)
                                                                         aCCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATcc                                                   >  1:15377/1‑150 (MQ=255)
                                                                         aCCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCca                                                  <  2:1080336/151‑1 (MQ=255)
                                                                           cATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCcac                                                 >  2:885202/1‑150 (MQ=255)
                                                                                cACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGGTATAAACATCCACAGGGAc                                           >  1:349459/1‑151 (MQ=255)
                                                                                cACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAggg                                             >  2:786334/1‑149 (MQ=255)
                                                                                cACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGa                                            <  2:568258/150‑1 (MQ=255)
                                                                                cACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGAc                                           <  1:592948/151‑1 (MQ=255)
                                                                                 aCCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCTGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAACCATCCACAGGGAc                                           <  2:167546/150‑1 (MQ=255)
                                                                                 aCCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGAc                                           >  1:672632/1‑150 (MQ=255)
                                                                                  ccATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACaa                                         <  1:988184/151‑1 (MQ=255)
                                                                                   cATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAAt                                        >  2:16151/1‑151 (MQ=255)
                                                                                   cATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAAt                                        >  2:1176718/1‑151 (MQ=255)
                                                                                           aTCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCg                                >  1:266816/1‑151 (MQ=255)
                                                                                             ctcCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGt                               >  1:560469/1‑150 (MQ=255)
                                                                                              tcCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTAc                             >  1:661567/1‑151 (MQ=255)
                                                                                               ccgccgGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCg                                <  1:1094704/147‑1 (MQ=255)
                                                                                               ccgccgGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTa                              <  2:1042894/149‑1 (MQ=255)
                                                                                                cgccgGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCg                           <  1:519574/151‑1 (MQ=255)
                                                                                                cgccgGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCg                           >  1:22635/1‑151 (MQ=255)
                                                                                                cgccgGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCg                           >  1:159035/1‑151 (MQ=255)
                                                                                                   cgGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGccc                        >  1:521704/1‑151 (MQ=255)
                                                                                                   cgGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGACATAAACATCCACAGGGACAATTTTATCCGTACCGccc                        >  1:843525/1‑151 (MQ=255)
                                                                                                          tCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACg                  <  1:1157306/150‑1 (MQ=255)
                                                                                                             tCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAg               <  1:451105/150‑1 (MQ=255)
                                                                                                             tCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGt              >  2:298840/1‑151 (MQ=255)
                                                                                                              ccAGTTCACCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAg               >  2:809010/1‑149 (MQ=255)
                                                                                                                     aCCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTag            >  2:1200140/1‑145 (MQ=255)
                                                                                                                     aCCCGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTag            >  2:919082/1‑145 (MQ=255)
                                                                                                                      cccGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTAGAGATCg      <  1:1033612/150‑1 (MQ=255)
                                                                                                                      cccGGCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTAGAGATCGt     >  2:813146/1‑151 (MQ=255)
                                                                                                                        cGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAACCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTAGAGATCGTgg   >  1:461873/1‑151 (MQ=255)
                                                                                                                          gCCCACGGGCGTGAATTTTCTGTTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTAGAGATCGTGGa  >  1:289313/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                               
GGCAATCTGACGACCGTAACGATAACCCGCCAGGCGACGTGCTTCGCGATCCACTTTCACGCGCAGCGGCTGCACCATATCACCATGCAGGATCTCCGCCGGTTGTTCATCCAGTTCACCCGGCCCACGGGCGTGAATTTTCTGCTCCAGCAGGCCGAGCGCCATTGCAAAGCCGTACAGCGTGGCGGCAGGCGTTGGCGGGCAGCCAGGGATATAAACATCCACAGGGACAATTTTATCCGTACCGCCCCACACGCAGTAGAGATCGTGGA  >  NC_000913/2843575‑2843846

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