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,3590CT97.3% 380.4 / 0.0 113A13T (GCT→ACT) metJtranscriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (59/50);  total (60/53)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 113 total aligned reads displayed.

GACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTAC  >  NC_000913/4128263‑4128455
                                                                                                |                                                                                                
gacgacGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTa                                                                                               <  1:867129/100‑1 (MQ=255)
gacgacGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACAGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTa                                                                                               >  1:2424018/1‑100 (MQ=255)
 acgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTAt                                                                                              <  2:1272372/100‑1 (MQ=255)
  cgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATg                                                                                             >  2:2171524/1‑100 (MQ=255)
   gacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATgg                                                                                            >  1:257285/1‑100 (MQ=255)
     cgcgTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGc                                                                                          <  2:633471/100‑1 (MQ=255)
      gcgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCt                                                                                         <  2:557747/100‑1 (MQ=255)
      gcgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCt                                                                                         <  1:666074/100‑1 (MQ=255)
       cgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTg                                                                                        >  2:1563904/1‑100 (MQ=255)
       cgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTg                                                                                        <  2:2563897/100‑1 (MQ=255)
         tgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGat                                                                                      <  2:1512024/100‑1 (MQ=255)
          gcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGata                                                                                     >  1:1228956/1‑100 (MQ=255)
          gcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGata                                                                                     <  2:1561266/100‑1 (MQ=255)
           cgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGatat                                                                                    <  2:874606/100‑1 (MQ=255)
           cgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGatat                                                                                    >  1:820788/1‑100 (MQ=255)
           cgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGatat                                                                                    <  2:530881/100‑1 (MQ=255)
              tCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATAtt                                                                                 <  2:637078/100‑1 (MQ=255)
               cATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTc                                                                                <  1:1388168/100‑1 (MQ=255)
                aTCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcg                                                                               >  1:2614272/1‑100 (MQ=255)
                  cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcgcc                                                                             <  1:2747858/100‑1 (MQ=255)
                  cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcgcc                                                                             <  1:1854493/100‑1 (MQ=255)
                   ggTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcgccg                                                                            >  1:2047458/1‑100 (MQ=255)
                   ggTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTcgccg                                                                            >  1:457378/1‑100 (MQ=255)
                     tGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTTACTTTTTTACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCt                                                                          >  2:2462694/1‑100 (MQ=255)
                     tGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCt                                                                          >  1:1216952/1‑100 (MQ=255)
                      gAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTc                                                                         >  2:2379041/1‑100 (MQ=255)
                       aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTcc                                                                        >  1:1523706/1‑100 (MQ=255)
                       aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTcc                                                                        <  2:829717/100‑1 (MQ=255)
                        ggATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCa                                                                       >  1:1639041/1‑100 (MQ=255)
                        ggATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCa                                                                       >  1:1384598/1‑100 (MQ=255)
                        ggATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCa                                                                       >  2:2032440/1‑100 (MQ=255)
                        ggATTTTTAACACCCTAAGAGGAATGGAAACCGTAATTTTTTTAACTTGTTCACTCTTCTTACCGTGCTCATTTTATGGGCTTCTATATTctcccc                                                                           >  2:1404533/1‑94 (MQ=255)
                           tttttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTc                                                                    >  2:1272603/1‑100 (MQ=255)
                           tttttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTc                                                                    <  2:1486774/100‑1 (MQ=255)
                             tttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAg                                                                  >  1:2225306/1‑100 (MQ=255)
                              ttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGc                                                                 <  2:157251/100‑1 (MQ=255)
                               tAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGcc                                                                <  2:1020321/100‑1 (MQ=255)
                               tAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGcc                                                                <  1:902866/100‑1 (MQ=255)
                                  cacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATg                                                             >  1:1924408/1‑100 (MQ=255)
                                   acCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATga                                                            <  1:204269/100‑1 (MQ=255)
                                   acCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATga                                                            >  2:345931/1‑100 (MQ=255)
                                    ccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgag                                                           <  1:1404541/100‑1 (MQ=255)
                                     cTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgaga                                                          >  1:481221/1‑100 (MQ=255)
                                       tAAGAGGAATGGAAACCGTAATCTTTTTTTTTTTTTTTTTCTTCTTGCCGTGCTCAGTGTATGGTCTGATATATTTGCTGCTCCATTCAGCCATGAGATa                                                        >  2:2054270/1‑100 (MQ=255)
                                       tAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATa                                                        >  2:2390430/1‑100 (MQ=255)
                                         agagGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACt                                                      >  1:2759013/1‑100 (MQ=255)
                                          gagGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACt                                                      >  2:1623586/1‑99 (MQ=255)
                                            ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                                   >  1:1076918/1‑100 (MQ=255)
                                            ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                                   >  2:700930/1‑100 (MQ=255)
                                             gAATGGAAACCGTAATCTTTTTGACTTGTTCTCTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAAt                                                  <  1:124175/100‑1 (MQ=255)
                                             gAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAAt                                                  >  1:2913783/1‑100 (MQ=255)
                                             gAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAAt                                                  <  2:1884309/100‑1 (MQ=255)
                                              aaTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAAt                                                  <  2:1237389/99‑1 (MQ=255)
                                              aaTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATc                                                 <  2:2485741/100‑1 (MQ=255)
                                                 ggAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctc                                              >  1:259631/1‑100 (MQ=255)
                                                 ggAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctc                                              >  1:300222/1‑100 (MQ=255)
                                                  gAAACCGTAATCTTTTTGACTTGTTCAGTCTTCTTGCCATGCTCACTGTATAGGCTGATATATAAGCAGGCCGATTCAGCCATGAGATAGTTAATActc                                              >  2:2861795/1‑99 (MQ=255)
                                                  gAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGTGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctct                                             >  2:2949077/1‑100 (MQ=255)
                                                  gAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctct                                             <  2:2424003/100‑1 (MQ=255)
                                                  gAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctct                                             >  1:785566/1‑100 (MQ=255)
                                                      ccGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCg                                          >  2:1359232/1‑99 (MQ=255)
                                                      ccGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGt                                         >  2:1358368/1‑100 (MQ=255)
                                                      ccGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGt                                         >  1:2577130/1‑100 (MQ=255)
                                                       cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                        <  2:259629/100‑1 (MQ=255)
                                                       cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                        <  1:1572717/100‑1 (MQ=255)
                                                       cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                        <  2:1513409/100‑1 (MQ=255)
                                                       cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGGGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                        <  1:87834/100‑1 (MQ=255)
                                                        gTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCa                                       >  2:1455618/1‑100 (MQ=255)
                                                        gTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCa                                       <  2:1364339/100‑1 (MQ=255)
                                                              tttttGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAt                                 >  1:2077094/1‑100 (MQ=255)
                                                              tttttGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAt                                 >  2:2286078/1‑100 (MQ=255)
                                                                tttGACTTGTACACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCTATCAGATACTTAATCCTCTTCGTCAATAAATTg                               <  2:1962964/100‑1 (MQ=255)
                                                                  tGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgag                             <  2:2776840/100‑1 (MQ=255)
                                                                  tGACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgag                             <  2:1228950/100‑1 (MQ=255)
                                                                   gACTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgaga                            <  1:2116734/100‑1 (MQ=255)
                                                                    aCTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagac                           <  2:1779123/100‑1 (MQ=255)
                                                                     cTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                          <  2:2913767/100‑1 (MQ=255)
                                                                     cTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                          >  1:2110585/1‑100 (MQ=255)
                                                                     cTTGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                          >  2:252678/1‑100 (MQ=255)
                                                                      ttGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacca                         <  1:1198645/100‑1 (MQ=255)
                                                                       tGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccag                        >  2:2012439/1‑100 (MQ=255)
                                                                       tGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccag                        <  2:2614257/100‑1 (MQ=255)
                                                                       tGTTCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccag                        >  2:3613/1‑100 (MQ=255)
                                                                          tCACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagacc                     <  2:743923/100‑1 (MQ=255)
                                                                           cACTCTTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAgaccatacca                    >  1:2682008/1‑100 (MQ=255)
                                                                             ctctTCTTGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcac                   <  1:35089/99‑1 (MQ=255)
                                                                               cttcttGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacaca                <  1:1913375/100‑1 (MQ=255)
                                                                               cttcttGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacaca                <  1:2171536/100‑1 (MQ=255)
                                                                                  cttGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGtt             <  1:337598/100‑1 (MQ=255)
                                                                                   ttGCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTg            <  2:2455083/100‑1 (MQ=255)
                                                                                     gCCGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAt          >  2:874396/1‑100 (MQ=255)
                                                                                      ccGTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtg         <  2:1523697/100‑1 (MQ=255)
                                                                                        gTGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtgtg       <  1:1563913/100‑1 (MQ=255)
                                                                                         tGCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTgg      <  2:2633834/100‑1 (MQ=255)
                                                                                          gCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTggg     <  1:592321/100‑1 (MQ=255)
                                                                                          gCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTggg     >  2:1983478/1‑100 (MQ=255)
                                                                                          gCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTggg     <  2:2536809/100‑1 (MQ=255)
                                                                                          gCTCAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTggg     >  1:1515502/1‑100 (MQ=255)
                                                                                           cTCAGTGTATGGGATGATATATTCGCCGCTCCATTCAGGCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGt    >  2:1839098/1‑100 (MQ=255)
                                                                                             cAGTGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTAc  >  2:2223181/1‑100 (MQ=255)
                                                                                                |                                                                                                
GACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTAC  >  NC_000913/4128263‑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: