Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,128,359 C→T A13T (GCT→ACT)  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,3590CT99.2% 406.8 / NA 120A13T (GCT→ACT) metJtranscriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (57/62);  total (58/62)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

TGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTAC  >  NC_000913/4128262‑4128455
                                                                                                 |                                                                                                
tGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAgtgt                                                                                                <  1:103131/100‑1 (MQ=255)
 gacgacGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTa                                                                                               >  2:3185490/1‑100 (MQ=255)
     acgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATggg                                                                                           >  1:2619631/1‑100 (MQ=255)
       gcgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCt                                                                                         <  2:2958427/100‑1 (MQ=255)
        cgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGTCTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTg                                                                                        >  1:1693820/1‑100 (MQ=255)
         gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTgn                                                                                       <  2:1652269/100‑2 (MQ=255)
         gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGa                                                                                       <  2:1689785/100‑1 (MQ=255)
         gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGa                                                                                       <  2:1075708/100‑1 (MQ=255)
          tgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGat                                                                                      >  2:844577/1‑100 (MQ=255)
           gcgtTCATCGTTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCCCTCTTCTTGCCGTGCTCAGTGTATGGGCTGata                                                                                     <  2:2147550/100‑1 (MQ=255)
            cgtTCATCGGTGAGTATTTTTAACACCTTAAGAGGAATGTAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGat                                                                                      <  2:3249181/98‑1 (MQ=255)
            cgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGatat                                                                                    >  2:1393942/1‑100 (MQ=255)
              ccgatctgTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGatatat                                                                                  <  2:3260060/93‑1 (MQ=255)
                cATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTTTTTTTTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTc                                                                                >  2:650738/1‑100 (MQ=255)
                cATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTc                                                                                >  1:607293/1‑100 (MQ=255)
                 aTCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcg                                                                               <  2:1335139/100‑1 (MQ=255)
                 aTCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcg                                                                               >  1:1850198/1‑100 (MQ=255)
                 aTCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcg                                                                               <  2:79124/100‑1 (MQ=255)
                  tCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcgc                                                                              <  2:515531/100‑1 (MQ=255)
                   cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcgcc                                                                             <  2:1474556/100‑1 (MQ=255)
                   cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcgcc                                                                             <  2:3135454/100‑1 (MQ=255)
                   cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcgcc                                                                             <  1:845286/100‑1 (MQ=255)
                   cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcgcc                                                                             <  2:2307426/100‑1 (MQ=255)
                    ggTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcgccg                                                                            >  1:3042961/1‑100 (MQ=255)
                     gTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGatatatagatcgg                                                                           >  1:3260099/1‑93 (MQ=255)
                     gTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcgccgc                                                                           >  1:3066873/1‑100 (MQ=255)
                       gAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTc                                                                         >  1:932995/1‑100 (MQ=255)
                        aGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCGTTTTGACTTGTTCACTCTTCGTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTcc                                                                        <  2:1944598/100‑1 (MQ=255)
                        aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTcc                                                                        >  2:2062215/1‑100 (MQ=255)
                        aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTcc                                                                        >  1:212411/1‑100 (MQ=255)
                           aTTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCAtt                                                                     >  2:2521872/1‑100 (MQ=255)
                            tttttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTc                                                                    <  2:2631708/100‑1 (MQ=255)
                            tttttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTc                                                                    <  1:798184/100‑1 (MQ=255)
                              tttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAg                                                                  <  2:2497450/100‑1 (MQ=255)
                              tttAACACCTTAAGAGGAATGGAAACCGTAATATTTTTGACTTGTCCACTCTTCTTGCCGTGCTCAGTGTATGGGATGATATATTCGCCGCTCCATTCAg                                                                  >  2:2759751/1‑100 (MQ=255)
                                 aaCACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCa                                                               >  1:2928343/1‑100 (MQ=255)
                                  acacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCAt                                                              >  2:797412/1‑100 (MQ=255)
                                   cacCTTAAGAGGAATGGAAACCGTAATCTTTTTTTTTTTTTTTCTCTTCTTACCGTGCTCAGTGTATGGTCTGATATATTCGCCGCTCCATTCAGCCAt                                                              >  2:3031522/1‑99 (MQ=255)
                                   cacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATg                                                             <  2:2619598/100‑1 (MQ=255)
                                   cacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATg                                                             <  1:1910134/100‑1 (MQ=255)
                                    acCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATga                                                            <  1:2284580/100‑1 (MQ=255)
                                    acCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATga                                                            <  2:3193651/100‑1 (MQ=255)
                                     ccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgag                                                           >  1:2716310/1‑100 (MQ=255)
                                     ccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgag                                                           >  2:389172/1‑100 (MQ=255)
                                     ccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgag                                                           <  1:1936882/100‑1 (MQ=255)
                                       ttAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGAt                                                         >  1:213094/1‑100 (MQ=255)
                                         aaGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATAc                                                       <  1:559785/100‑1 (MQ=255)
                                           gagGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCCGATATATTCGCCGCTCCATTCAGCCATGAGATACtt                                                     <  1:2299975/100‑1 (MQ=255)
                                             ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                                   >  2:1704243/1‑100 (MQ=255)
                                             ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                                   >  1:3242166/1‑100 (MQ=255)
                                             ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                                   >  2:1855089/1‑100 (MQ=255)
                                              gAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAAt                                                  <  1:1368868/100‑1 (MQ=255)
                                                aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                                <  2:1645456/100‑1 (MQ=255)
                                                aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCCCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                                <  2:2490741/100‑1 (MQ=255)
                                                  ggAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctc                                              >  2:375780/1‑100 (MQ=255)
                                                  ggAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctc                                              >  2:112183/1‑100 (MQ=255)
                                                    aaaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCtt                                            <  1:1440767/100‑1 (MQ=255)
                                                    aaaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCtt                                            <  1:389179/100‑1 (MQ=255)
                                                      aCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCg                                          <  2:1669173/100‑1 (MQ=255)
                                                       ccGTAATCTTTGTGACTTGTTCACTCTTCTTGCCGTTCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGt                                         <  2:2858856/100‑1 (MQ=255)
                                                        cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                        <  1:1304528/100‑1 (MQ=255)
                                                        cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                        <  2:2860501/100‑1 (MQ=255)
                                                        cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                        <  1:1393956/100‑1 (MQ=255)
                                                         gTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCa                                       <  1:687732/100‑1 (MQ=255)
                                                          tAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCaa                                      <  1:844587/100‑1 (MQ=255)
                                                             tCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATaa                                   <  2:686974/100‑1 (MQ=255)
                                                               tttttGACTTGGTCACTCGTCTTGCGGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCGTCTTCGTCAATAAAt                                 >  2:1467997/1‑100 (MQ=255)
                                                                ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                                <  1:1240401/100‑1 (MQ=255)
                                                                ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                                <  2:2038247/100‑1 (MQ=255)
                                                                ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                                <  1:1240402/100‑1 (MQ=255)
                                                                ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                                <  1:1498853/100‑1 (MQ=255)
                                                                ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                                <  1:668130/100‑1 (MQ=255)
                                                                  ttGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTga                              >  2:1047087/1‑100 (MQ=255)
                                                                  ttGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTga                              <  1:367227/100‑1 (MQ=255)
                                                                    gACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCTGCTCCATTCAGCCATGACATACTTAATCCTCTTCGTCAATAAATTgaga                            >  2:1430965/1‑100 (MQ=255)
                                                                    gACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgaga                            >  1:288462/1‑100 (MQ=255)
                                                                      cTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                          <  2:172591/100‑1 (MQ=255)
                                                                      cTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                          >  1:2906320/1‑100 (MQ=255)
                                                                      cTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                          >  2:2941898/1‑100 (MQ=255)
                                                                         gTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccaga                       >  2:583261/1‑100 (MQ=255)
                                                                           tCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagacc                     >  2:704809/1‑100 (MQ=255)
                                                                           tCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagacc                     <  1:3126115/100‑1 (MQ=255)
                                                                            cACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAgaccagacca                    >  2:1595711/1‑100 (MQ=255)
                                                                               tcttctTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacac                 <  2:2102308/100‑1 (MQ=255)
                                                                               tcttctTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacac                 <  2:160803/100‑1 (MQ=255)
                                                                               tcttcgTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcaca                  <  2:1176580/99‑1 (MQ=255)
                                                                                cttcttGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacaca                <  1:2112589/100‑1 (MQ=255)
                                                                                   cttGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGtt             >  2:109997/1‑100 (MQ=255)
                                                                                   cttGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGtt             <  2:533545/100‑1 (MQ=255)
                                                                                     tGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGa           >  1:968455/1‑100 (MQ=255)
                                                                                     tGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGa           <  1:3185527/100‑1 (MQ=255)
                                                                                      gCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAt          >  2:1762384/1‑100 (MQ=255)
                                                                                       ccGTGCTCAGTGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtg         <  2:3066836/100‑1 (MQ=255)
                                                                                       ccGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCGCTTCGTCAATAAATTGAGACCAGAGCACACAGTTGAtg         <  2:2015390/100‑1 (MQ=255)
                                                                                        cGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtgt        >  1:1735813/1‑100 (MQ=255)
                                                                                         gTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtgtg       >  1:1979551/1‑100 (MQ=255)
                                                                                           gCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCCTCAATAAATTGAGACCAGACCACACAGTTGATGTggg     <  2:1746676/100‑1 (MQ=255)
                                                                                            cTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAGATTGAGACCAGACCACACAGTTGATGTGGGt    >  2:1233236/1‑100 (MQ=255)
                                                                                             tCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTa   >  1:833741/1‑100 (MQ=255)
                                                                                              cAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTAc  >  2:2665273/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTAC  >  NC_000913/4128262‑4128455

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: