Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,128,352 T→A 100% H15L (CAC→CTC)  metJ ← transcriptional repressor, S‑adenosylmethionine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,128,3520TA97.8% 452.4 / ‑0.2 151H15L (CAC→CTC) metJtranscriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base A (70/76);  total (73/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

GTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATG  >  NC_000913/4128255‑4128448
                                                                                                 |                                                                                                
gttCACCTGATGACTCGTGCGTTCATCGGTGATGATTTTTAACACTTTAAGAGGAATGGAAACCGTAATCTTTCTGACTTGTTCACTCTTCTTGCCGAGc                                                                                                <  1:2375873/100‑1 (MQ=255)
 ttCACCTGACGACGCGTGCGTTCATCGGTTAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGc                                                                                                <  1:2463328/99‑1 (MQ=255)
 ttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCt                                                                                               <  2:1814713/100‑1 (MQ=255)
   cACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCa                                                                                             >  1:1989353/1‑100 (MQ=255)
   cACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAAAAGAAATATTTTTGACTTGTTCACTCTTCTTGCCGAGCTCa                                                                                             >  1:1451429/1‑100 (MQ=255)
     ccTTACTACGCTTGCTTTCATCGTTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAg                                                                                            <  1:1732721/99‑1 (MQ=255)
     ccTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCCGc                                                                                           >  2:635715/1‑100 (MQ=255)
        gacgacGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTa                                                                                        >  1:159677/1‑100 (MQ=255)
        gacgacGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTa                                                                                        >  2:2478392/1‑100 (MQ=255)
        gacgacGCGTGCGTTCATCGGAGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCg                                                                                          >  2:612771/1‑98 (MQ=255)
          cgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTAt                                                                                       >  1:3564908/1‑99 (MQ=255)
          cgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATg                                                                                      >  1:2584795/1‑100 (MQ=255)
           gacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATgg                                                                                     <  2:1145953/100‑1 (MQ=255)
            acgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTATTCTTTTTGAGTTGTTCACTCTTCTTGCCGAGCTCAGCGTATggg                                                                                    <  2:1382567/100‑1 (MQ=255)
              gcgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCt                                                                                  <  1:3809236/100‑1 (MQ=255)
               cgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTg                                                                                 >  1:2790465/1‑100 (MQ=255)
               cgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTg                                                                                 <  2:230496/100‑1 (MQ=255)
                gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATATTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGa                                                                                >  1:503715/1‑100 (MQ=255)
                 tgtgtTCATCGGTGAGGATTTTTAAAACCTTAAGAAGAATGGAAACCGTAATCTTTTTGACTTGATAAATGTTCTTGCTGAGCTCAGCGTATGGTCTTat                                                                               >  2:3401090/4‑100 (MQ=255)
                  gcgtTCATCGGTGAGGATTTTTAACTCCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGata                                                                              <  1:2193212/100‑1 (MQ=255)
                  gcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGata                                                                              <  2:3845738/100‑1 (MQ=255)
                  gcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAAGGGAAGCCGTAGTCTTTTTGACTTGTTCACTCTTCTTGCCGAGGTCAGCGTATGGGCTGat                                                                               <  2:1887957/99‑1 (MQ=255)
                   cgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGGAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGatat                                                                             >  2:2171383/1‑100 (MQ=255)
                    gtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGatata                                                                            <  2:2020539/100‑1 (MQ=255)
                       cATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTc                                                                         >  2:2420522/1‑100 (MQ=255)
                       cATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTc                                                                         <  2:1275062/100‑1 (MQ=255)
                        aTCGGTGAGGATTTTTAACACCTTAAGAGGAATGGATACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcg                                                                        <  1:168759/100‑1 (MQ=255)
                         tCGTTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgc                                                                       <  2:2561480/100‑1 (MQ=255)
                         tCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgc                                                                       <  1:1430994/100‑1 (MQ=255)
                           ggTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgccg                                                                     >  1:1884499/1‑100 (MQ=255)
                           ggTGAGGAGTTTTATCATCTTAAGAGGAATGGAAGTCGTAATCTTTTTGGCTTGTTCACTCGTCTTGCCGTGCTCAGCTCATGCGCTGagctcgg                                                                          <  2:159156/95‑7 (MQ=255)
                            gTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgccgc                                                                    >  1:1320573/1‑100 (MQ=255)
                            gTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgccgc                                                                    >  2:763765/1‑100 (MQ=255)
                             tGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGTCTGATATATTCGCCGCt                                                                   <  2:2452867/100‑1 (MQ=255)
                              gAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTc                                                                  >  1:1507965/1‑100 (MQ=255)
                               aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTcc                                                                 <  1:3639420/100‑1 (MQ=255)
                                ggATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCa                                                                >  2:2623197/1‑100 (MQ=255)
                                ggATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGATCAGCGTATGGGCTGATATATTCGCCGCTCCa                                                                >  1:3016998/1‑100 (MQ=255)
                                ggATTTTTAACACCTTAAGAGGAATAAAAACCGTAATCTTTTTAACTTGTTCACTCTTCTTGCCGAGCTCAAAGTATGGGTTGCTATATTCCCCGCTcc                                                                 >  2:2925247/1‑99 (MQ=255)
                                 gATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAt                                                               <  1:2420544/100‑1 (MQ=255)
                                  aTTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAtt                                                              >  1:398293/1‑100 (MQ=255)
                                  aTTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAtt                                                              >  1:3453518/1‑100 (MQ=255)
                                  aTTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAtt                                                              >  1:2102044/1‑100 (MQ=255)
                                    ttttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCa                                                            >  1:2907982/1‑100 (MQ=255)
                                      ttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGc                                                          >  1:1768843/1‑100 (MQ=255)
                                      ttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGc                                                          <  2:379746/100‑1 (MQ=255)
                                       tAACACCTTAAGAGTAATGGAAACCGTATTCTTTTTTACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGcc                                                         <  1:835638/100‑1 (MQ=255)
                                       tAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGcc                                                         <  2:2857400/100‑1 (MQ=255)
                                       tAACACCTTAAGAGGAATAGAAACCGGAATCTTTTTGACTTGTTCACTCTTCTTGCCGGGCTCAGCGTATGGGCTGATGGATTCGCCTCTCCATTCAGcc                                                         >  2:1750794/1‑100 (MQ=255)
                                         acacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCAt                                                       >  1:305586/1‑100 (MQ=255)
                                          cgatctAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATg                                                      <  1:900209/95‑1 (MQ=255)
                                          cacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATAGGCTGATATATTCGCCGCTCCATTCAGCCATg                                                      >  1:679963/1‑100 (MQ=255)
                                            ccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgag                                                    >  2:225447/1‑100 (MQ=255)
                                            ccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGAGCTGATATATTCGCCGCTCCATTCAGCCATgag                                                    >  1:2572943/1‑100 (MQ=255)
                                              ttAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGAt                                                  <  1:1209374/100‑1 (MQ=255)
                                               naaGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATa                                                 >  1:1544216/2‑100 (MQ=255)
                                                aaGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATAc                                                <  1:2447052/100‑1 (MQ=255)
                                                  gagGAATGGAAACCGTAAGCTTTGTGACTTTGTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTGGTCGCTCCATTCAGCCATGAGATACtt                                              <  2:1696090/100‑1 (MQ=255)
                                                    ggAATGGGAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                            >  1:2051744/1‑100 (MQ=255)
                                                    ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                            >  2:1717941/1‑100 (MQ=255)
                                                    ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                            >  2:2251021/1‑100 (MQ=255)
                                                    ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGATCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                            >  2:1192835/1‑100 (MQ=255)
                                                      aaTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATc                                          <  1:2213343/100‑1 (MQ=255)
                                                       aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTTAGCCATGAGATACTTAATcc                                         <  2:262697/100‑1 (MQ=255)
                                                       aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                         <  1:2816387/100‑1 (MQ=255)
                                                       aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                         <  2:3205289/100‑1 (MQ=255)
                                                       aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                         <  1:740709/100‑1 (MQ=255)
                                                        tGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCct                                        >  1:2724640/1‑100 (MQ=255)
                                                             aCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCg                                   <  1:107596/100‑1 (MQ=255)
                                                              ccGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGt                                  >  2:53056/1‑100 (MQ=255)
                                                              ccGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGt                                  >  1:2149865/1‑100 (MQ=255)
                                                              ccGTAATCTATATGACAAGTTCACTCGTCTGGCGGAGCCGAGCGTATGGGCTGTTATATTCGCCGCTCCATTGAGCCATGGGATACTTaagcgtc                                       >  2:320428/1‑90 (MQ=255)
                                                               cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                 <  2:2493101/100‑1 (MQ=255)
                                                                gTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCTATGAGATACTTAATCCTCTTCGTc                                 <  2:2340901/99‑1 (MQ=255)
                                                                gTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCa                                <  2:2877286/100‑1 (MQ=255)
                                                                 tAATCTTTTTGACTTGTTCACTCTTCTTGCCAAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACCTAATCCTCTTCGTCaa                               >  2:2287848/1‑100 (MQ=255)
                                                                      tttttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAt                          >  1:3012212/1‑100 (MQ=255)
                                                                        tttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                        <  1:1717956/100‑1 (MQ=255)
                                                                        tttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                        <  1:3899914/100‑1 (MQ=255)
                                                                         ttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTga                       >  1:3451150/1‑100 (MQ=255)
                                                                           gACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTTCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgaga                     <  2:789574/100‑1 (MQ=255)
                                                                           gACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgaga                     <  1:3689314/100‑1 (MQ=255)
                                                                             cTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                   >  2:494381/1‑100 (MQ=255)
                                                                             cTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                   <  1:430963/100‑1 (MQ=255)
                                                                             cTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagac                    >  2:3546720/1‑99 (MQ=255)
                                                                               tGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTTATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccag                 <  1:1008422/100‑1 (MQ=255)
                                                                               tGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccag                 >  2:1852069/1‑100 (MQ=255)
                                                                                gTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccaga                <  2:1027062/100‑1 (MQ=255)
                                                                                  tCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagacc              >  2:1418069/1‑100 (MQ=255)
                                                                                  tCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagacc              <  2:976579/100‑1 (MQ=255)
                                                                                   cACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAgaccagacca             <  2:1551246/100‑1 (MQ=255)
                                                                                    aCTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcac            >  2:3320676/1‑100 (MQ=255)
                                                                                    aCTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcac            >  2:13977/1‑100 (MQ=255)
                                                                                     ctctTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcaca           >  2:1935307/1‑100 (MQ=255)
                                                                                        ttcttGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCTATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAg        <  1:3038796/100‑1 (MQ=255)
                                                                                           ttGCCGAGCTCAGGGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTg     <  2:1935831/100‑1 (MQ=255)
                                                                                             gCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAt   >  2:3383652/1‑100 (MQ=255)
                                                                                             gCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAt   <  2:2930571/100‑1 (MQ=255)
                                                                                             gCCGAGCTCAGCGTATGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtg  <  1:1836752/100‑1 (MQ=255)
                                                                                             gCCGAGCGCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAt   >  1:2057787/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATG  >  NC_000913/4128255‑4128448

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