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,3520TA98.3% 398.4 / ‑2.3 116H15L (CAC→CTC) metJtranscriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (58/56);  total (60/56)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 116 total aligned reads displayed.

TCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATG  >  NC_000913/4128257‑4128448
                                                                                               |                                                                                                
tCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTc                                                                                              <  2:2090019/100‑1 (MQ=255)
 cACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCa                                                                                             >  1:132118/1‑100 (MQ=255)
 cACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCa                                                                                             >  2:1723976/1‑100 (MQ=255)
  aCCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAg                                                                                            >  1:1309223/1‑100 (MQ=255)
  aCCTGACGAAGCGTGCGTTCATCGGTGAGGATTTTTAACACCATAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAg                                                                                            <  1:2462894/100‑1 (MQ=255)
     tGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGt                                                                                         >  1:994071/1‑100 (MQ=255)
     tGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGt                                                                                         >  2:1266494/1‑100 (MQ=255)
     tGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGt                                                                                         >  2:2571481/1‑100 (MQ=255)
      gacgacGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTa                                                                                        >  1:974388/1‑100 (MQ=255)
      gacgacGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTa                                                                                        <  1:2259434/100‑1 (MQ=255)
       acgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTAt                                                                                       >  1:202083/1‑100 (MQ=255)
        cgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATg                                                                                      <  1:2641907/100‑1 (MQ=255)
           cgcgTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGc                                                                                   >  1:1821166/1‑100 (MQ=255)
            gcgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCt                                                                                  <  2:538937/100‑1 (MQ=255)
            gcgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCt                                                                                  >  1:213321/1‑100 (MQ=255)
              gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGa                                                                                >  1:1175715/1‑100 (MQ=255)
               tgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGat                                                                               <  2:213320/100‑1 (MQ=255)
                tccgatctCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGata                                                                              <  2:1518555/93‑1 (MQ=255)
                gcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGata                                                                              <  2:2194153/100‑1 (MQ=255)
                 cgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGatat                                                                             >  2:1443672/1‑100 (MQ=255)
                  gtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAACCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGatata                                                                            <  2:1199899/100‑1 (MQ=255)
                  gtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGatata                                                                            <  2:1810669/100‑1 (MQ=255)
                      aTCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGatat                                                                             <  2:2972007/95‑1 (MQ=255)
                       tCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgc                                                                       >  2:1085766/1‑100 (MQ=255)
                        cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgcc                                                                      <  1:505963/100‑1 (MQ=255)
                        cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgcc                                                                      <  1:658378/100‑1 (MQ=255)
                        cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgcc                                                                      >  2:959561/1‑100 (MQ=255)
                         ggTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgccg                                                                     >  2:401015/1‑100 (MQ=255)
                          gTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTcgccgc                                                                    <  2:907117/100‑1 (MQ=255)
                           tGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCt                                                                   >  1:1136674/1‑100 (MQ=255)
                            gAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTc                                                                  <  2:1040420/100‑1 (MQ=255)
                            gAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTc                                                                  >  2:1657330/1‑100 (MQ=255)
                             aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTcc                                                                 <  2:1526945/100‑1 (MQ=255)
                             aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTcc                                                                 <  2:778579/100‑1 (MQ=255)
                              ggATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTTTTTTTTTTTCTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCa                                                                >  2:826572/1‑100 (MQ=255)
                              ggATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCa                                                                >  1:2107253/1‑100 (MQ=255)
                               gATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAt                                                               <  1:2300368/100‑1 (MQ=255)
                               gATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAt                                                               >  1:2056823/1‑100 (MQ=255)
                                aTTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAtt                                                              >  1:2736708/1‑100 (MQ=255)
                                aTTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAtt                                                              >  2:942120/1‑100 (MQ=255)
                                   tttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAg                                                           <  1:1468165/100‑1 (MQ=255)
                                    ttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGc                                                          >  2:2064202/1‑100 (MQ=255)
                                    ttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGc                                                          >  1:2410366/1‑100 (MQ=255)
                                      aaCACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCa                                                        <  2:1586287/100‑1 (MQ=255)
                                       acacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCAt                                                       >  1:1456732/1‑100 (MQ=255)
                                        cacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTCCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATg                                                      >  1:1370238/1‑100 (MQ=255)
                                          ccTTAAGAGGAATGGAAACCGTAATCTTTTTTTTTTTTTTTTTCTTCTTGCCGAGCTCAGCGTATGGTCTGATATATTCGCCGCTCCATTCAGCCATgag                                                    >  2:2864038/1‑100 (MQ=255)
                                          ccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgag                                                    >  1:1133307/1‑100 (MQ=255)
                                          ccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgag                                                    >  1:2284081/1‑100 (MQ=255)
                                           cTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgaga                                                   >  1:611508/1‑100 (MQ=255)
                                            ttAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGAt                                                  <  2:2622522/100‑1 (MQ=255)
                                            ttAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGAt                                                  <  2:2688799/100‑1 (MQ=255)
                                             tAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATa                                                 <  1:1555363/100‑1 (MQ=255)
                                             tAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATa                                                 >  1:1409793/1‑100 (MQ=255)
                                                gagGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACtt                                              <  2:11926/100‑1 (MQ=255)
                                                  ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                            <  1:2547782/100‑1 (MQ=255)
                                                  ggAATGGAAACCGTAATCTTGTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                            >  2:1006892/1‑100 (MQ=255)
                                                    aaTGGAAACCGTAATCTTTTTGACTTGTTCCCTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCACCGCTCCATTCAGCCATGAGATACTTAATc                                          >  1:1090172/1‑100 (MQ=255)
                                                    aaTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATc                                          >  1:2560298/1‑100 (MQ=255)
                                                     ctGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                         <  1:2930947/99‑1 (MQ=255)
                                                     aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                         <  1:733050/100‑1 (MQ=255)
                                                     aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                         <  2:1872316/100‑1 (MQ=255)
                                                     aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                         <  2:1514680/100‑1 (MQ=255)
                                                     aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                         <  1:2788655/100‑1 (MQ=255)
                                                      tGGAAGCGGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCct                                        <  1:2416491/100‑1 (MQ=255)
                                                      tGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCct                                        >  1:2310139/1‑100 (MQ=255)
                                                      tGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCct                                        >  2:1830677/1‑100 (MQ=255)
                                                       ggAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctc                                       >  1:1067687/1‑100 (MQ=255)
                                                       ggAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCAAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                         >  2:1210778/1‑98 (MQ=255)
                                                        gAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctct                                      <  1:172151/100‑1 (MQ=255)
                                                         aaaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCtt                                     <  1:3031989/100‑1 (MQ=255)
                                                           aCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCg                                   <  1:2528723/100‑1 (MQ=255)
                                                             cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                 <  1:954894/100‑1 (MQ=255)
                                                              gTAATCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCa                                >  1:2549347/1‑100 (MQ=255)
                                                                  tCTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATaa                            <  2:3081209/100‑1 (MQ=255)
                                                                   cTTTTTGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATaaa                           >  1:1989135/1‑100 (MQ=255)
                                                                    tttttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAt                          >  1:2287749/1‑100 (MQ=255)
                                                                     ttttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                         <  2:354973/100‑1 (MQ=255)
                                                                     ttttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                         <  1:1087356/100‑1 (MQ=255)
                                                                     ttttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                         <  1:1087363/100‑1 (MQ=255)
                                                                      tttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                        >  2:1675750/1‑100 (MQ=255)
                                                                      tttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                        <  1:959572/100‑1 (MQ=255)
                                                                       ttGACTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTga                       <  2:2182007/100‑1 (MQ=255)
                                                                        tGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgag                      >  2:97562/1‑100 (MQ=255)
                                                                         gACTTGTTCACTCTTCTTGCCGAGATCAGCGTATCGGCTGAGCTATTAGCCGCTACATTAAGCCATGAGATATTTAATCCTCTTCGTCAATAAATAgaga                     >  2:2365038/1‑100 (MQ=255)
                                                                           cTTGTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                   >  2:12532/1‑100 (MQ=255)
                                                                              gTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccaga                >  2:2500825/1‑100 (MQ=255)
                                                                              gTTCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGCTACTTACTCCTCTTCGTCACTAAATTGagaccaga                >  1:1800366/1‑100 (MQ=255)
                                                                               ttCACTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac               >  1:2901466/1‑100 (MQ=255)
                                                                                  aCTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcac            <  1:353387/100‑1 (MQ=255)
                                                                                  aCTCTTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcac            >  1:2502875/1‑100 (MQ=255)
                                                                                   ctctTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcaca           >  2:27006/1‑100 (MQ=255)
                                                                                   ctctTCTTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcaca           >  2:1520186/1‑100 (MQ=255)
                                                                                    tcttctTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacac          <  1:2181855/100‑1 (MQ=255)
                                                                                    tcttctTGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacac          <  1:1617487/100‑1 (MQ=255)
                                                                                      ttcttGCTGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTTAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGAACAGACCACACAg        >  2:2586441/1‑100 (MQ=255)
                                                                                         ttGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTg     <  2:1309212/100‑1 (MQ=255)
                                                                                         ttGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTg     <  2:611506/100‑1 (MQ=255)
                                                                                         ttGCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTg     <  1:2311548/100‑1 (MQ=255)
                                                                                           gCCGAGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAt   >  1:541029/1‑100 (MQ=255)
                                                                                               |                                                                                                
TCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATG  >  NC_000913/4128257‑4128448

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