Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,128,352 T→A 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,3520TA99.4% 584.2 / ‑6.0 174H15L (CAC→CTC) metJtranscriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (86/87);  minor base C (1/0);  total (87/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 174 total aligned reads displayed.

GTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAT  >  NC_000913/4128255‑4128447
                                                                                                 |                                                                                               
  tCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTc                                                                                             >  1:3943655/1‑100 (MQ=255)
   cACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCa                                                                                            >  2:731021/1‑100 (MQ=255)
   cACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCa                                                                                            >  2:668381/1‑100 (MQ=255)
   cACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCa                                                                                            >  2:1548384/1‑100 (MQ=255)
   cACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCa                                                                                            >  1:3404878/1‑100 (MQ=255)
   cACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCa                                                                                            >  1:3180352/1‑100 (MQ=255)
    aCCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAg                                                                                           >  1:113254/1‑100 (MQ=255)
     ccTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGc                                                                                          >  2:3581011/1‑100 (MQ=255)
     ccTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGc                                                                                          <  2:2790069/100‑1 (MQ=255)
         acgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTAt                                                                                      <  2:88747/100‑1 (MQ=255)
         acgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTAt                                                                                      <  1:3543778/100‑1 (MQ=255)
             cgcgTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGc                                                                                  >  1:828577/1‑100 (MQ=255)
              tccgatcttCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCt                                                                                 <  2:2807051/93‑1 (MQ=255)
              gcgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCt                                                                                 <  1:1650650/100‑1 (MQ=255)
                gtgcgtTCATCGGTGAGGATTTTTTACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGa                                                                               <  1:3076960/100‑1 (MQ=255)
                gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGa                                                                               >  1:4165655/1‑100 (MQ=255)
                gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGa                                                                               >  1:2121506/1‑100 (MQ=255)
                gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGa                                                                               <  2:2949087/100‑1 (MQ=255)
                 tgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGat                                                                              <  2:4127090/100‑1 (MQ=255)
                  gcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGata                                                                             <  2:4108264/100‑1 (MQ=255)
                   cgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATgagatcggaa                                                                            >  2:1936688/1‑91 (MQ=255)
                      tCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCtagatcgga                                                                         >  1:2807083/1‑92 (MQ=255)
                       cATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTc                                                                        >  2:1315803/1‑100 (MQ=255)
                        aTCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcg                                                                       >  1:3069826/1‑100 (MQ=255)
                          cGGTGAGGATTTTTAACACCTTATGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgc                                                                      <  2:1563044/99‑1 (MQ=255)
                          cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgcc                                                                     <  1:1986908/100‑1 (MQ=255)
                           ggTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgccg                                                                    >  1:2967690/1‑100 (MQ=255)
                            gTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgccgc                                                                   <  2:1957068/100‑1 (MQ=255)
                            gTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgc                                                                      >  2:3451174/1‑97 (MQ=255)
                             tGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCt                                                                  <  1:267127/100‑1 (MQ=255)
                               aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTcc                                                                <  2:3895971/100‑1 (MQ=255)
                                ggATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCa                                                               >  1:1228226/1‑100 (MQ=255)
                                ggATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCa                                                               <  1:763703/100‑1 (MQ=255)
                                ggATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCa                                                               >  1:593787/1‑100 (MQ=255)
                                  aTTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTTTTTTTTTTTTTATTCTAACCGAGCTCAGCGTATTGGTTAATATATTCGCCTCTCCAtt                                                             >  2:816907/1‑100 (MQ=255)
                                  aTTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGTCTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAtt                                                             >  2:3726120/1‑100 (MQ=255)
                                  aTTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAtt                                                             <  2:3114485/100‑1 (MQ=255)
                                     tttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAg                                                          <  1:2101651/100‑1 (MQ=255)
                                     tttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAg                                                          <  1:334563/100‑1 (MQ=255)
                                      ttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGc                                                         >  1:234635/1‑100 (MQ=255)
                                       tAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCGTTCAGcc                                                        <  1:2199686/100‑1 (MQ=255)
                                       tAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGcc                                                        >  1:3582718/1‑100 (MQ=255)
                                       tAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGcc                                                        >  1:1834020/1‑100 (MQ=255)
                                       gatctccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGcn                                                        <  2:228112/95‑2 (MQ=255)
                                        aaCACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGcc                                                        >  2:2096043/1‑99 (MQ=255)
                                         acacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGCGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCAt                                                      >  2:2668668/1‑100 (MQ=255)
                                         acacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCAt                                                      <  1:4126779/100‑1 (MQ=255)
                                          cacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATg                                                     >  2:2356615/1‑100 (MQ=255)
                                          cacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATg                                                     >  2:3783868/1‑100 (MQ=255)
                                           acCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATga                                                    <  1:1386018/100‑1 (MQ=255)
                                            ccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCagatc                                                   >  1:228117/1‑96 (MQ=255)
                                             cTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgaga                                                  <  1:2459431/100‑1 (MQ=255)
                                             cTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgaga                                                  <  2:3690686/100‑1 (MQ=255)
                                             cTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgaga                                                  >  1:2084822/1‑100 (MQ=255)
                                              ttAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGAt                                                 >  1:41452/1‑100 (MQ=255)
                                               tAAGAGGAATGGAAACCGTAATCTTTTTTTTTTTTTTACTCTTCTTGCCGAGCTCAACGTATGGGCTGATATATTTGCCGCTCCATTCAGCCATGAGAt                                                 >  2:1851359/1‑99 (MQ=255)
                                                 agagGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACt                                              <  1:2515319/100‑1 (MQ=255)
                                                 agagGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACt                                              >  2:1747789/1‑100 (MQ=255)
                                                 agagGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACt                                              <  2:3296548/100‑1 (MQ=255)
                                                 agagGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACt                                              <  2:1039935/100‑1 (MQ=255)
                                                  gagGAATGTAAACCGTAATCTTTTTGACGTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACtt                                             <  2:2510939/100‑1 (MQ=255)
                                                   agGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTa                                            >  1:2228672/1‑100 (MQ=255)
                                                    ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                           >  1:1838631/1‑100 (MQ=255)
                                                    ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                           >  2:2730351/1‑100 (MQ=255)
                                                    ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                           >  2:2442425/1‑100 (MQ=255)
                                                    ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                           >  2:1028277/1‑100 (MQ=255)
                                                    ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTAGGGGCGGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                           >  2:345140/1‑100 (MQ=255)
                                                     gAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAAt                                          >  1:699310/1‑100 (MQ=255)
                                                       aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                        <  2:3213173/100‑1 (MQ=255)
                                                       aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                        <  1:222710/100‑1 (MQ=255)
                                                       aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                        >  2:2599677/1‑100 (MQ=255)
                                                          gAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctct                                     >  2:2266576/1‑100 (MQ=255)
                                                          gAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctct                                     <  2:306625/100‑1 (MQ=255)
                                                            aaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTc                                   <  2:1228216/100‑1 (MQ=255)
                                                             aCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCg                                  <  2:2580586/100‑1 (MQ=255)
                                                               cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                <  2:1740292/100‑1 (MQ=255)
                                                               cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                <  1:1417001/100‑1 (MQ=255)
                                                                 tAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCaa                              >  1:1118436/1‑100 (MQ=255)
                                                                 tAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCaa                              <  2:1273549/100‑1 (MQ=255)
                                                                  aaTCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAAt                             <  2:113254/100‑1 (MQ=255)
                                                                      tttttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAt                         <  1:4042004/100‑1 (MQ=255)
                                                                       ttttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                        <  1:2709867/100‑1 (MQ=255)
                                                                       ttttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                        >  1:438592/1‑100 (MQ=255)
                                                                       ttttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                        <  1:634889/100‑1 (MQ=255)
                                                                          tGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgag                     <  1:3726161/100‑1 (MQ=255)
                                                                          tGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgag                     <  2:2967658/100‑1 (MQ=255)
                                                                            aCTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagac                   <  1:2075324/100‑1 (MQ=255)
                                                                            aCTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagac                   <  2:4165610/100‑1 (MQ=255)
                                                                            aCTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagac                   <  2:234630/100‑1 (MQ=255)
                                                                            aCTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagac                   >  1:1498444/1‑100 (MQ=255)
                                                                             cTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                  <  2:438588/100‑1 (MQ=255)
                                                                                 ttCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac              >  1:3951354/1‑100 (MQ=255)
                                                                                 ttCACTCTTCCTGCCGAGCTCAGCGTATGGGCTGATATATTCGCTGCTTCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccaga               <  2:3929255/99‑1 (MQ=255)
                                                                                      tcttctTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTTAATAAATTGAGACCAGACcacac         <  2:3404840/100‑1 (MQ=255)
                                                                                      tcttctTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacac         <  2:1407547/100‑1 (MQ=255)
                                                                                      tctgctTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagacc             <  1:96347/96‑1 (MQ=255)
                                                                                        ttcttGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAg       >  2:3371744/1‑100 (MQ=255)
                                                                                         tcttGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGt      <  2:1083241/100‑1 (MQ=255)
                                                                                            tGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGa   <  2:2917250/100‑1 (MQ=255)
                                                                                             gCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAt  <  1:2442450/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAT  >  NC_000913/4128255‑4128447

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