Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,845,792 C→T 5.4% R227H (CGT→CAT)  hycE ← hydrogenase 3, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,845,7920CT 5.4% 305.7 / 5.4 93R227H (CGT→CAT) hycEhydrogenase 3, large subunit
Reads supporting (aligned to +/- strand):  ref base C (44/44);  new base T (3/2);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

AGAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACGTTGTTTT  >  NC_000913/2845648‑2845939
                                                                                                                                                |                                                                                                                                                   
agAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGt                                                                                                                                                 <  4:281683/149‑1 (MQ=255)
     ggCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCg                                                                                                                                            >  3:229405/1‑149 (MQ=255)
      gcgcGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGc                                                                                                                                           >  1:2696/1‑149 (MQ=255)
        gcgGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCa                                                                                                                                         >  2:650172/1‑149 (MQ=255)
         cgGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAg                                                                                                                                        <  1:418362/149‑1 (MQ=255)
           gATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGtt                                                                                                                                      >  2:250592/1‑149 (MQ=255)
             tcatcTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGtttt                                                                                                                                    <  2:218000/149‑1 (MQ=255)
             tcatcTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGtgtt                                                                                                                                    >  3:481603/1‑149 (MQ=255)
              catcTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGttttt                                                                                                                                   <  3:595938/149‑1 (MQ=255)
               atcTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTc                                                                                                                                  <  1:315616/149‑1 (MQ=255)
               atcTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTc                                                                                                                                  <  2:79679/149‑1 (MQ=255)
               atcTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTc                                                                                                                                  >  4:164901/1‑149 (MQ=255)
                 cTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCa                                                                                                                                <  2:310589/149‑1 (MQ=255)
                   gcgcACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATg                                                                                                                              <  2:81476/149‑1 (MQ=255)
                    cgcACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGc                                                                                                                             <  1:683290/149‑1 (MQ=255)
                    cgcACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGc                                                                                                                             <  4:118048/149‑1 (MQ=255)
                         gTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCTCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGa                                                                                                                        <  4:845088/149‑1 (MQ=255)
                             tGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGa                                                                                                                    <  1:250592/149‑1 (MQ=255)
                                 aCCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCGTATGGGttata                                                                                                                                             >  3:681978/1‑117 (MQ=255)
                                 aCCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCGTATGGGttata                                                                                                                                             <  4:681978/120‑4 (MQ=255)
                                 aCCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTa                                                                                                                >  1:656987/1‑149 (MQ=255)
                                 aCCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTa                                                                                                                <  2:624623/149‑1 (MQ=255)
                                     gAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAAc                                                                                                            <  4:589782/149‑1 (MQ=255)
                                        tACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGa                                                                                                         <  4:848577/149‑1 (MQ=255)
                                         aCCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGAc                                                                                                        <  4:832565/149‑1 (MQ=255)
                                           ccATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACgg                                                                                                      >  1:61181/1‑149 (MQ=255)
                                            cATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGt                                                                                                     <  2:567958/149‑1 (MQ=255)
                                              tCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAg                                                                                                   >  4:307222/1‑149 (MQ=255)
                                              tCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAg                                                                                                   >  2:736965/1‑149 (MQ=255)
                                                gcgTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTc                                                                                                 <  1:812622/149‑1 (MQ=255)
                                                   ttttCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGc                                                                                              >  3:430531/1‑149 (MQ=255)
                                                   ttttCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGc                                                                                              >  1:656204/1‑149 (MQ=255)
                                                     ttCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGt                                                                                            <  4:259830/149‑1 (MQ=255)
                                                     ttCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGt                                                                                            <  2:853568/149‑1 (MQ=255)
                                                     ttCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGt                                                                                            <  4:583486/149‑1 (MQ=255)
                                                        cACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGa                                                                                         >  4:393671/1‑149 (MQ=255)
                                                         aCCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGAt                                                                                        >  4:231337/1‑149 (MQ=255)
                                                         aCCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGAt                                                                                        >  1:650895/1‑149 (MQ=255)
                                                           cGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATaa                                                                                      <  1:170403/149‑1 (MQ=255)
                                                            gACTTGGTGGAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGACAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAAt                                                                                     >  2:97104/1‑149 (MQ=255)
                                                               gtggtgTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATATGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGtt                                                                                  <  1:587707/149‑1 (MQ=255)
                                                                tggtgTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGttt                                                                                 <  4:792662/149‑1 (MQ=255)
                                                                tggtgTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGttt                                                                                 >  2:590342/1‑149 (MQ=255)
                                                                  gtgtAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTc                                                                               <  3:164901/149‑1 (MQ=255)
                                                                       ggCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGGTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCAt                                                                          >  2:2346/1‑149 (MQ=255)
                                                                             gCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATcgacga                                                                    >  3:646012/1‑149 (MQ=255)
                                                                                gtgGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACa                                                                 <  2:656204/149‑1 (MQ=255)
                                                                                      aaaGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGa                                                           >  1:28746/1‑149 (MQ=255)
                                                                                       aaGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGaa                                                          >  2:546316/1‑149 (MQ=255)
                                                                                       aaGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGaa                                                          <  3:374038/149‑1 (MQ=255)
                                                                                           cGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTgg                                                      <  4:660839/149‑1 (MQ=255)
                                                                                           cGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTgg                                                      <  2:403964/149‑1 (MQ=255)
                                                                                            gCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATcgac                                                                      <  1:211307/132‑1 (MQ=255)
                                                                                            gCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATcgac                                                                      >  2:211307/1‑132 (MQ=255)
                                                                                            gCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATcga                                                                       <  1:169356/131‑1 (MQ=255)
                                                                                            gCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATcga                                                                       >  2:169356/1‑131 (MQ=255)
                                                                                            gCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGc                                                     <  1:724145/149‑1 (MQ=255)
                                                                                            gCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGc                                                     >  4:847693/1‑149 (MQ=255)
                                                                                                aTCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCgg                                                 >  2:384099/1‑149 (MQ=255)
                                                                                                       acacacGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGTTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTtcataa                                          >  2:544263/1‑147 (MQ=255)
                                                                                                         acacGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGa                                        >  1:381507/1‑149 (MQ=255)
                                                                                                          cacGGTCAGAGAGGAAGGTCACTTCGTTATACCCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGaa                                       <  4:481603/149‑1 (MQ=255)
                                                                                                          cacGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCGATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGaa                                       <  3:522613/149‑1 (MQ=255)
                                                                                                          cacGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTATCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGaa                                       <  1:650172/149‑1 (MQ=255)
                                                                                                            cGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGt                                     <  1:66873/149‑1 (MQ=255)
                                                                                                             ggTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTg                                    >  3:524964/1‑149 (MQ=255)
                                                                                                              gTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGa                                   <  3:55644/149‑1 (MQ=255)
                                                                                                               tCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGAc                                  <  4:229405/149‑1 (MQ=255)
                                                                                                                   agagGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGc                              >  1:161542/1‑149 (MQ=255)
                                                                                                                     agGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAg                            >  1:54272/1‑149 (MQ=255)
                                                                                                                         aGGTCACTTCGTTATAACCCATGCGGGTTTCCGCCAGTTTTCCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACAAGCAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGa                        <  2:711426/149‑1 (MQ=255)
                                                                                                                          ggTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGAc                       >  1:658872/1‑149 (MQ=255)
                                                                                                                            tCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCa                     <  1:590342/149‑1 (MQ=255)
                                                                                                                             cACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCaa                    >  2:25731/1‑149 (MQ=255)
                                                                                                                              aCTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAAt                   <  2:54272/149‑1 (MQ=255)
                                                                                                                               cTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATc                  >  1:656003/1‑149 (MQ=255)
                                                                                                                                ttCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACGGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCg                 >  2:373979/1‑149 (MQ=255)
                                                                                                                                ttCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCg                 >  1:836277/1‑149 (MQ=255)
                                                                                                                                ttCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCg                 >  4:619812/1‑149 (MQ=255)
                                                                                                                                   gTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTtc                                              >  2:220872/1‑117 (MQ=255)
                                                                                                                                   gTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTtc                                              <  1:220872/117‑1 (MQ=255)
                                                                                                                                     tataACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTc                                                                               >  2:76782/1‑82 (MQ=255)
                                                                                                                                     tataACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTc                                                                               <  1:76782/82‑1 (MQ=255)
                                                                                                                                     ggaaaCCCGTATGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCacg            >  1:726658/4‑149 (MQ=255)
                                                                                                                                     ggaaaCCCGTATGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAAAAGTGACGTGCAGCGGACCAATCGGCacg            >  4:553226/4‑149 (MQ=255)
                                                                                                                                      ataACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCacga           <  4:646012/149‑1 (MQ=255)
                                                                                                                                           ccATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACgttg      >  1:112555/1‑149 (MQ=255)
                                                                                                                                           ccATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACgttg      >  3:701492/1‑149 (MQ=255)
                                                                                                                                           ccATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCACCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACgttg      >  2:418204/1‑149 (MQ=255)
                                                                                                                                             aTACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACgttgtt    >  3:743392/1‑149 (MQ=255)
                                                                                                                                              tACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGTTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACGTTGttt   <  2:75542/149‑1 (MQ=255)
                                                                                                                                              tACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACGTTGttt   >  2:832490/1‑149 (MQ=255)
                                                                                                                                               aCGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAAACGGCACGACGTTGtttt  >  2:356849/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AGAATGGCGCGGATCATCTGCGCACGTTCTGGCACCTGAATACCCATCGCGTTTTCCACCGACGTGGTGTAGGCGGTGCTGTGGGCAAAGCCGCAGATCCCGCACACACGGTCAGAGAGGAAGGTCACTTCGTTATAACCCATACGGGTTTCCGCCAGTTTTTCCATGCCGCGATGGACGTAGAACAGACGGTAGTCGGCGTCGATAATGTTTTCGCCATCGACGAACAGACGGAAGTGGCCCGGTTCATCAGAAGTGACGTGCAGCGGACCAATCGGCACGACGTTGTTTT  >  NC_000913/2845648‑2845939

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