Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,225,139 T→C intergenic (+49/‑171) metH → / → yjbB B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase/putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,225,1390TC99.4% 527.4 / NA 168intergenic (+49/‑171)metH/yjbBB12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase/putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (97/70);  minor base A (1/0);  total (98/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 213 total aligned reads displayed.

GAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAG  >  NC_000913/4225049‑4225226
                                                                                          |                                                                                        
gggcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCCATGAgcgc                                                                                        >  2:1559241/3‑93 (MQ=255)
gggcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCGTACAACAAACAGGGCGCTCAATGAgcgc                                                                                        >  2:6445527/3‑93 (MQ=255)
     ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGt                                                                                   >  2:3298523/1‑93 (MQ=255)
      tggctggcACCGAATCTGGGGTAGGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtc                                                                                  <  2:5091990/93‑1 (MQ=255)
       ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtct                                                                                 >  1:2877446/1‑93 (MQ=255)
          tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCtt                                                                              >  1:4160024/1‑93 (MQ=255)
           ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGTGTGCTCAATGAGCTACCTCTCTCttt                                                                             >  1:4857454/1‑93 (MQ=255)
              cccGAATCTGGGGTATGTCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTttatt                                                                          <  2:2091556/92‑1 (MQ=255)
              aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCtttttt                                                                          >  1:1676451/1‑93 (MQ=255)
              aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTttatt                                                                          >  1:6427279/1‑93 (MQ=255)
              aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTttatt                                                                          >  1:3369096/1‑93 (MQ=255)
               ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTttatta                                                                         <  1:861534/93‑1 (MQ=255)
                 gAATCGGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAat                                                                       <  1:6367356/93‑1 (MQ=255)
                      tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaa                                                                  >  1:5901182/1‑93 (MQ=255)
                      tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaa                                                                  >  1:6574036/1‑93 (MQ=255)
                      tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaa                                                                  >  1:621526/1‑93 (MQ=255)
                       ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaaa                                                                 >  1:2049119/1‑93 (MQ=255)
                        gggTATGACGCGGACTGATTCACAAATCTGTCATTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAac                                                                <  1:3115023/93‑1 (MQ=255)
                        gggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAac                                                                <  1:4483038/93‑1 (MQ=255)
                          gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAacac                                                              >  2:6664101/1‑93 (MQ=255)
                           taccaCGCGGCGTGATTCACAAATCTGTCACTTTTCATTACAACACACAGGCCGCTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACt                                                            <  2:307836/89‑1 (MQ=255)
                             tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCCGTCTCCTTATTAATATGAAAAAATta                                                          >  1:2455944/1‑93 (MQ=255)
                             tGACGCGGACTGATTCACAAATCTGTCACTTTCCCTTACAACAAACAGGGCGCTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTta                                                          <  2:3858635/93‑1 (MQ=255)
                               aCGCGGACTGATCCACAAATCTGTCACTTTTCCTTACAACAAACGGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTtata                                                        <  1:6092098/93‑1 (MQ=255)
                                cgcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                       >  1:5101192/1‑93 (MQ=255)
                                 gcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                      >  2:2649726/1‑93 (MQ=255)
                                  gggACTTATTCACAACTTTTTCTTTTTTCCTAAAAAAAAACGGGCCCCTAAAGGCCCCCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                     <  1:3691856/92‑1 (MQ=255)
                                  cgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGGGCCATGAGCGCCCCCTTTCTTTTTTAATATTAAAAAAATATACTg                                                     >  2:3581482/1‑93 (MQ=255)
                                   ggACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                    >  2:5303180/1‑93 (MQ=255)
                                    gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCCATGAGCGCCCCGTCTTTTTTTTTATATGAAAAACTTATTCTTGa                                                   >  1:3574857/1‑93 (MQ=255)
                                    gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATTAAACACTTATACTGGa                                                   >  2:2549414/1‑93 (MQ=255)
                                     aCTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCCGTCTCTTTATTAATATGAAAAACTTATACTTGaa                                                  >  1:1335191/1‑93 (MQ=255)
                                      cTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa                                                 <  2:3991885/93‑1 (MQ=255)
                                       tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                >  2:4425307/1‑93 (MQ=255)
                                       tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                >  2:4299732/1‑93 (MQ=255)
                                       tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                >  1:6717826/1‑93 (MQ=255)
                                        gATTCACAAATCTGTCTCTTTTCCTAACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                               <  2:5750599/93‑1 (MQ=255)
                                        gATTCACAAATCGTTCACTTTTCCTTACAACAAACAGGGCGTTCATTGACCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                               <  2:6277845/93‑1 (MQ=255)
                                         aTTCACAAATCTGTCATTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAg                                              <  2:4162160/93‑1 (MQ=255)
                                         aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAg                                              <  1:2964095/93‑1 (MQ=255)
                                         aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAg                                              >  1:490197/1‑93 (MQ=255)
                                           tCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGc                                            >  2:6123211/1‑93 (MQ=255)
                                             acaAGTCTGTCACTTTTCCTTACACAAAACAGGGCTCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTg                                          <  1:4108107/93‑1 (MQ=255)
                                             acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAAAAGCTg                                          >  1:640000/1‑93 (MQ=255)
                                             acaAATCTGTCACTTGCCCTTACAACAAGCAGGGCGTCCAAGGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTg                                          <  2:1944979/93‑1 (MQ=255)
                                              caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                         >  1:460917/1‑93 (MQ=255)
                                              caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                         <  1:1919086/93‑1 (MQ=255)
                                              caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                         >  2:2203663/1‑93 (MQ=255)
                                               aaaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                        >  1:5318925/1‑93 (MQ=255)
                                               aaaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGAGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                        >  2:2979137/1‑93 (MQ=255)
                                                aaTCTGTCACTTTTCCTTACAACAAACGGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                       <  2:1357200/93‑1 (MQ=255)
                                                   ctctCATTTTTCCTCACACAAACCGGGGCGCTATATGAGCACCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                    <  2:579619/90‑1 (MQ=255)
                                                    tGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGgataat                                   >  2:6273139/1‑93 (MQ=255)
                                                    tGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGACCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataat                                   <  2:4860584/93‑1 (MQ=255)
                                                     gTCTCTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                  <  1:5077758/93‑1 (MQ=255)
                                                      tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                 <  1:1883635/93‑1 (MQ=255)
                                                       taCTTTTCCGTACAACAAACAGGGCGCTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCtt                                <  1:2012722/92‑1 (MQ=255)
                                                       cACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCtt                                <  1:2219181/93‑1 (MQ=255)
                                                        aCTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTTGAATAATCTTa                               >  2:380999/1‑93 (MQ=255)
                                                        aCTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                               >  2:4474294/1‑93 (MQ=255)
                                                        aCTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                               <  2:6729268/93‑1 (MQ=255)
                                                        aCTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                               <  1:898615/93‑1 (MQ=255)
                                                         cTTTTCCTTACAACAAACAGGGGGGTCCATGAAGGCCCCGTCTCTTTTTTTATATTAAAAACATTTAATTGAAAAAAGGCGGGAATAATTTTa                               >  2:4660900/1‑93 (MQ=255)
                                                          ttttCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                             >  1:2516993/1‑93 (MQ=255)
                                                          ttttCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                             >  2:5549863/1‑93 (MQ=255)
                                                          ttttCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                             >  2:3385880/1‑93 (MQ=255)
                                                           tttCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGcc                            >  1:2907420/1‑93 (MQ=255)
                                                            ttCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                           >  1:3044218/1‑93 (MQ=255)
                                                            taccTTACAACAAACAGGGCCCTAAATGAGACCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                           <  2:3925058/91‑1 (MQ=255)
                                                               cTTACAACCAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                        >  1:6543179/1‑93 (MQ=255)
                                                               cTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                        <  2:1773778/93‑1 (MQ=255)
                                                               cTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                        >  2:5465961/1‑93 (MQ=255)
                                                               cTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                        >  1:5530649/1‑93 (MQ=255)
                                                                 tACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaa                      >  2:4089833/1‑93 (MQ=255)
                                                                  acaacaAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                     <  2:2222320/93‑1 (MQ=255)
                                                                    aacaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGGCGGAAAAACg                   >  2:6082023/1‑93 (MQ=255)
                                                                    aacaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                   >  1:958850/1‑93 (MQ=255)
                                                                    aacaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                   >  1:5421107/1‑93 (MQ=255)
                                                                     acaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGGATAATCTTAGACGGAAAAAAGa                  >  1:3829208/1‑93 (MQ=255)
                                                                     acaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGa                  >  2:5587519/1‑93 (MQ=255)
                                                                      caaacaGGGCGTTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                 <  1:275681/93‑1 (MQ=255)
                                                                      caaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGGATAATCTTAGCCGGAAAAACGaa                 >  2:3319587/1‑93 (MQ=255)
                                                                      caaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                 <  1:706869/93‑1 (MQ=255)
                                                                       aaacgGGGCGCACAAGGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                <  1:4711553/93‑1 (MQ=255)
                                                                        aacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg               >  2:2605806/1‑93 (MQ=255)
                                                                        aacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg               >  1:5571170/1‑93 (MQ=255)
                                                                         acggggCGCTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt              <  2:3654731/90‑1 (MQ=255)
                                                                         acaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGGATAATCTTAGCCCGAAAAACGAAAGt              >  2:201673/1‑93 (MQ=255)
                                                                         acaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt              >  2:1105991/1‑93 (MQ=255)
                                                                              gcgcTCAAGGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa         <  2:600311/93‑1 (MQ=255)
                                                                               cgcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg        >  2:6223546/1‑93 (MQ=255)
                                                                               cgcTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg        <  2:6251019/93‑1 (MQ=255)
                                                                                gcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg       >  2:150168/1‑93 (MQ=255)
                                                                                 cTCAATGAGCTCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc      <  1:6408148/93‑1 (MQ=255)
                                                                                 cTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc      >  1:1004756/1‑93 (MQ=255)
                                                                                  taaaTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGca     <  2:4917644/91‑1 (MQ=255)
                                                                                  tCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGca     >  1:4023941/1‑93 (MQ=255)
                                                                                  tCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGca     >  1:3152986/1‑93 (MQ=255)
                                                                                    aaTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca   >  2:3163508/1‑93 (MQ=255)
                                                                                     aTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg  <  1:4826804/93‑1 (MQ=255)
                                                                                          |                                                                                        
GAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAG  >  NC_000913/4225049‑4225226

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