Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,011,749 A→G A369A (GCA→GCG metE → tetrahydropteroyltriglutamate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,011,7490AG99.1% 651.6 / ‑3.3 211A369A (GCA→GCGmetEtetrahydropteroyltriglutamate methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (114/95);  total (115/96)
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.41e-01

Only 100 of 229 total aligned reads displayed.

TGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCAGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAT  >  NC_000913/4011659‑4011838
                                                                                          |                                                                                         
tGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACggcggc                                                                                         >  2:1622837/1‑93 (MQ=255)
tGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACggcggc                                                                                         >  2:2860833/1‑93 (MQ=255)
agCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGGCACACggcggc                                                                                         >  1:1171480/2‑93 (MQ=255)
  cAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGctc                                                                                       <  2:2591723/93‑1 (MQ=255)
   aGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGctct                                                                                      <  1:3489478/93‑1 (MQ=255)
    gAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTg                                                                                     <  1:2748098/93‑1 (MQ=255)
    gAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTg                                                                                     >  1:3139216/1‑93 (MQ=255)
      aGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGc                                                                                   >  1:2962751/1‑93 (MQ=255)
       gTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCa                                                                                  >  1:3812670/1‑93 (MQ=255)
       gTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCa                                                                                  >  2:2050496/1‑93 (MQ=255)
         gAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGATACGCCGGCTCTGGCaga                                                                                >  2:3740314/1‑93 (MQ=255)
           aaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGt                                                                              <  1:4441308/93‑1 (MQ=255)
            aaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTg                                                                             >  1:4360488/1‑93 (MQ=255)
              gCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGa                                                                           >  2:1907424/1‑93 (MQ=255)
               cTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAg                                                                          <  1:1823752/93‑1 (MQ=255)
                tGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgc                                                                         <  2:2827478/93‑1 (MQ=255)
                tGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgc                                                                         >  2:2249192/1‑93 (MQ=255)
                tGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgc                                                                         >  1:2639108/1‑93 (MQ=255)
                tGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgc                                                                         <  2:1455303/93‑1 (MQ=255)
                tGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGAACTTCTGCGCGATGCGCTGCACCGTGGTGACACGGCGGCGCTGTCAGAGTGGAgg                                                                         >  2:4572784/1‑92 (MQ=255)
                  gTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgcgc                                                                       >  1:2796455/1‑93 (MQ=255)
                   tttGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGcc                                                                      >  1:4085700/1‑93 (MQ=255)
                   tttGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGcc                                                                      >  1:4596378/1‑93 (MQ=255)
                     tGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGcccc                                                                    >  1:319298/1‑93 (MQ=255)
                       ccTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGGTCTGGCAGAGGGGAGCGCCCCGa                                                                  >  1:3669014/1‑93 (MQ=255)
                       ccTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCgc                                                                  >  1:186492/1‑92 (MQ=255)
                       ccTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGa                                                                  >  2:3507852/1‑93 (MQ=255)
                        cTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGAGGAGCGCCCCGAt                                                                 >  1:2679549/1‑93 (MQ=255)
                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGAtt                                                                >  1:1121152/1‑93 (MQ=255)
                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGAtt                                                                >  2:1961166/1‑93 (MQ=255)
                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGAtt                                                                <  2:3637012/93‑1 (MQ=255)
                          tCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTc                                                               <  1:1686522/93‑1 (MQ=255)
                             cccTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAgg                                                            >  2:1388496/1‑93 (MQ=255)
                             cccTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAgg                                                            <  1:3080304/93‑1 (MQ=255)
                              ccTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGc                                                           <  2:3246095/93‑1 (MQ=255)
                              ccTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGc                                                           >  2:1931274/1‑93 (MQ=255)
                               cTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCa                                                          >  2:1451931/1‑93 (MQ=255)
                               cTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCa                                                          >  1:419668/1‑93 (MQ=255)
                               cTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCa                                                          <  1:46700/93‑1 (MQ=255)
                                  aaaaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgt                                                       <  2:1736546/92‑1 (MQ=255)
                                   aaaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtc                                                      >  1:157015/1‑93 (MQ=255)
                                   aaaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtc                                                      <  2:2521943/93‑1 (MQ=255)
                                   aaaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtc                                                      <  1:2554169/93‑1 (MQ=255)
                                     aaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtcgt                                                    >  2:309116/1‑93 (MQ=255)
                                      aaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACgtcgtc                                                   >  2:4336410/1‑93 (MQ=255)
                                      aaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACgtcgtc                                                   <  2:4103713/93‑1 (MQ=255)
                                       aTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCa                                                  <  2:447297/93‑1 (MQ=255)
                                       aTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGGGGAGCGCCCCGATTCAGGCACgtcgtcc                                                  >  2:3385102/1‑92 (MQ=255)
                                        tGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCAc                                                 <  1:4498306/93‑1 (MQ=255)
                                         gCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGAGGAGCGCCCCGATTCAGGCACGTCGTCAct                                                >  1:633061/1‑93 (MQ=255)
                                          ccATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCActc                                               <  1:2850747/93‑1 (MQ=255)
                                          ccATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCActc                                               <  2:116133/93‑1 (MQ=255)
                                          ccATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCActc                                               <  2:1322625/93‑1 (MQ=255)
                                           cATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCActct                                              >  1:4517518/1‑93 (MQ=255)
                                           cATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCActct                                              >  2:3987271/1‑93 (MQ=255)
                                              gAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTAcc                                           >  1:4682784/1‑93 (MQ=255)
                                              gAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTAcc                                           >  1:1031644/1‑93 (MQ=255)
                                              gAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCACGATTCAGGCACGTCGTCACTCTAcc                                           <  2:943350/93‑1 (MQ=255)
                                               aaCTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTAccc                                          >  1:4000142/1‑93 (MQ=255)
                                                 cTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCcgc                                        <  2:710698/93‑1 (MQ=255)
                                                  tGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCcgcg                                       >  2:1562828/1‑93 (MQ=255)
                                                  cggCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCcgcg                                       <  1:1182510/92‑1 (MQ=255)
                                                   ggCACTGCTGCGCGATGCGCTGAACAGTGGTGACGCGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTCCCCGCGt                                      <  1:1391527/93‑1 (MQ=255)
                                                   ggCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGt                                      >  1:3220436/1‑93 (MQ=255)
                                                      aCTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACa                                   >  1:391983/1‑93 (MQ=255)
                                                      aCTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACa                                   <  1:3973428/93‑1 (MQ=255)
                                                      aCTGCGGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACa                                   <  2:3843109/93‑1 (MQ=255)
                                                       ctgctgCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACAt                                  >  1:511357/1‑93 (MQ=255)
                                                        tgctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATa                                 <  1:3241470/93‑1 (MQ=255)
                                                        tgctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATa                                 <  1:2011790/93‑1 (MQ=255)
                                                        tgctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATa                                 >  2:2063596/1‑93 (MQ=255)
                                                        tgctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCAGATTCAGGCACGTCGTCACTCTACCCGCGTAAATa                                 >  1:2563939/1‑93 (MQ=255)
                                                         gctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATaa                                >  1:1753343/1‑93 (MQ=255)
                                                          ctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAAt                               <  2:2340105/93‑1 (MQ=255)
                                                           tgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATc                              <  1:4134934/93‑1 (MQ=255)
                                                              gcgATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcgg                           <  2:2683788/93‑1 (MQ=255)
                                                                gATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggcg                         <  1:74651/93‑1 (MQ=255)
                                                                 aTGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggcgg                        >  1:543879/1‑93 (MQ=255)
                                                                  tGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCCCTCTACCCGCGTACATAATTCGGCTGt                       >  2:2871591/1‑93 (MQ=255)
                                                                  tGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGt                       >  1:4014969/1‑93 (MQ=255)
                                                                  tGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGt                       >  1:2086075/1‑93 (MQ=255)
                                                                    cgcTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAg                     <  2:2618064/93‑1 (MQ=255)
                                                                       tGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACTTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaa                  <  2:3658975/93‑1 (MQ=255)
                                                                       tGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaa                  >  2:1337503/1‑93 (MQ=255)
                                                                       tGAACAGTGGTGACACGGCGCCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaa                  >  2:4551571/1‑93 (MQ=255)
                                                                        gAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaaa                 <  2:2734710/93‑1 (MQ=255)
                                                                        gAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaaa                 >  1:2388437/1‑93 (MQ=255)
                                                                         aaCAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAg                <  2:1492274/93‑1 (MQ=255)
                                                                           cAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCg              >  2:157538/1‑93 (MQ=255)
                                                                            aGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGt             >  1:2995642/1‑93 (MQ=255)
                                                                             gtggtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTc            >  2:582152/1‑93 (MQ=255)
                                                                              tggtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCt           <  1:850375/93‑1 (MQ=255)
                                                                               ggtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTg          <  1:3522998/93‑1 (MQ=255)
                                                                               ggtgACACGGCGGCTCTGGCAGAGTGGAGCCCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTg          <  2:4169362/93‑1 (MQ=255)
                                                                                gtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTgg         >  2:2760779/1‑93 (MQ=255)
                                                                                gtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTgg         <  1:1310635/93‑1 (MQ=255)
                                                                                  gACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggcg       >  2:4140626/1‑93 (MQ=255)
                                                                                  gACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggcg       >  1:3319616/1‑93 (MQ=255)
                                                                                     acGGCAGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGgcggcg    <  2:1055197/93‑1 (MQ=255)
                                                                                       ggcggcTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAt  >  2:3557466/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCAGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAT  >  NC_000913/4011659‑4011838

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