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,1390TC92.6% 253.4 / 9.0 95intergenic (+49/‑171)metH/yjbBB12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase/putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base T (5/2);  new base C (46/42);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

Only 100 of 113 total aligned reads displayed.

GAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAG  >  NC_000913/4225049‑4225226
                                                                                          |                                                                                        
gggcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAgcgc                                                                                        >  1:1036598/3‑93 (MQ=255)
gggcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAgcgc                                                                                        >  1:1619385/3‑93 (MQ=255)
 ggcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGcc                                                                                       <  2:69503/92‑1 (MQ=255)
  gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGccc                                                                                      <  1:3501655/93‑1 (MQ=255)
  gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGccc                                                                                      <  1:2591310/93‑1 (MQ=255)
  gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGccc                                                                                      <  2:4554733/93‑1 (MQ=255)
      tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtc                                                                                  <  1:4694870/93‑1 (MQ=255)
      tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtc                                                                                  <  1:4005787/93‑1 (MQ=255)
      tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtc                                                                                  <  2:1617809/93‑1 (MQ=255)
       ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtct                                                                                 >  1:2311024/1‑93 (MQ=255)
        gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCTATGAGCGCCCTGtctc                                                                                >  1:4219409/1‑93 (MQ=255)
        gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCCCCCTGtctc                                                                                <  2:1637749/93‑1 (MQ=255)
         ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTctct                                                                               >  2:1105845/1‑93 (MQ=255)
         ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTctct                                                                               >  2:1228678/1‑93 (MQ=255)
          tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCtt                                                                              <  2:1673729/93‑1 (MQ=255)
          tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCtt                                                                              <  1:2499506/93‑1 (MQ=255)
           ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCttt                                                                             >  2:2103401/1‑93 (MQ=255)
           ggcACCGAAGCTGGGGTGTGACGCGGACTGATTCACAAATCTGTCATTTTGCATTACAACAAACAGGGCGTTCAATGAGCGCCCTGTCTCttt                                                                             <  2:1700663/93‑1 (MQ=255)
            gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                            <  1:740660/93‑1 (MQ=255)
            gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCtttt                                                                            >  2:2747760/1‑92 (MQ=255)
              aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTttatt                                                                          <  1:1539735/93‑1 (MQ=255)
                  aaTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCCATGAGCGCTCTGTCTCTTTATTtata                                                                      >  1:4220022/1‑93 (MQ=255)
                      tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaa                                                                  >  1:1245056/1‑93 (MQ=255)
                      tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaa                                                                  >  2:3284300/1‑93 (MQ=255)
                      tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaa                                                                  >  2:1528494/1‑93 (MQ=255)
                       ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaaa                                                                 <  2:3762880/93‑1 (MQ=255)
                        gggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAac                                                                <  1:3311950/93‑1 (MQ=255)
                          gTATGGCGGGAACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCTCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAacac                                                              <  2:2374154/93‑1 (MQ=255)
                           tATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACt                                                             <  1:2302347/93‑1 (MQ=255)
                            aTGACGAGCACTTATTCAAAACTCTAACACTTTACTTAACAACAACCGGGGCCCGCAAGGGGCGCCCTGTCTCTTTATTAATATGAAACACtt                                                            <  2:3591747/93‑1 (MQ=255)
                             tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGACCTGTCTCTTTATTAATATGAAACACTta                                                           >  1:1574009/1‑93 (MQ=255)
                               aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTtata                                                         <  1:3301162/93‑1 (MQ=255)
                                 gcgGACTGATCCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                       <  2:390043/93‑1 (MQ=255)
                                  cgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                      <  1:549841/93‑1 (MQ=255)
                                   ggACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                     >  2:4460082/1‑93 (MQ=255)
                                    gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                    >  1:3911562/1‑93 (MQ=255)
                                    gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                    >  1:3553357/1‑93 (MQ=255)
                                    gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCCCAAAGAGCGCCCTGGCCCCCTTTTAATAAGAAAAACATacacaaga                                                    >  1:4204617/1‑89 (MQ=255)
                                      cTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCCGTCTCTTTTTTAATATGAAACACTTATACTGGaaa                                                  >  2:142943/1‑93 (MQ=255)
                                       tGATTCACAAATCTGTCACTTTCCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                 <  2:3000715/93‑1 (MQ=255)
                                        gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                >  1:4572993/1‑93 (MQ=255)
                                          ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAgg                                              >  1:2198596/1‑93 (MQ=255)
                                          ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAgg                                              >  1:2623059/1‑93 (MQ=255)
                                           tCAAAAATCTGTCACTTTTCCTTACAACAAACGGGGCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGc                                             <  1:1154871/93‑1 (MQ=255)
                                             acaAATCTGTCACTTTTCCTTACAACAAACGGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTg                                          <  1:1336055/93‑1 (MQ=255)
                                             acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTg                                          <  2:4554091/93‑1 (MQ=255)
                                              caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                         <  1:636111/93‑1 (MQ=255)
                                              caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAAAGGCTgg                                         >  2:1271649/1‑93 (MQ=255)
                                               aaaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTCATATGAAACACTTATACTGGAAACAGGCTGGa                                        >  1:2967095/1‑93 (MQ=255)
                                                aaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                       >  1:314942/1‑93 (MQ=255)
                                                aaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                       <  2:4047678/93‑1 (MQ=255)
                                                aaTCCGCCACTTTTCCCTACAAAAAACGGGGCGCTCAACAATCCCCCTGTCTCTTTATTAATATGAAACCCTTATACTGGAAACAGGCTGGaa                                       <  1:1279865/93‑1 (MQ=255)
                                                 agcTGTCCTTTTTCCTAACAACAAAAGGGGCATGAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                      <  2:1303997/91‑1 (MQ=255)
                                                  tCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata                                     >  2:2459906/1‑93 (MQ=255)
                                                  tCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata                                     >  1:1156548/1‑93 (MQ=255)
                                                    tGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataat                                   >  1:489683/1‑93 (MQ=255)
                                                     gTCACTTTTCCTTACAACAAAGAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                  <  2:441335/93‑1 (MQ=255)
                                                     gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                  >  1:1661322/1‑93 (MQ=255)
                                                     gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                  >  2:2816590/1‑93 (MQ=255)
                                                      tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                 <  2:1562804/93‑1 (MQ=255)
                                                      tCACTGTTCCTTACAACAAACAGGGCGCTCAATGAGCTCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                 <  1:746581/93‑1 (MQ=255)
                                                       cACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGGATAATCtt                                >  2:2768917/1‑93 (MQ=255)
                                                        ccTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                               <  2:1609870/92‑1 (MQ=255)
                                                          ttttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                             <  2:2004668/93‑1 (MQ=255)
                                                          ttttCCTAACAACAAACAGGGCGCTAAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                             <  1:2035410/93‑1 (MQ=255)
                                                             tCCTTACAACAAACAGGGGGCTCCAAGAGGGCCCCGTCTCTTTTTTTATATTTAAAACTTATAATCGGAACAGGCTGGGAAAAACTTTAGCCg                           >  2:1876739/1‑93 (MQ=255)
                                                             tCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCgg                          <  1:3971618/93‑1 (MQ=255)
                                                             tCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCgg                          <  2:3851146/93‑1 (MQ=255)
                                                              ccTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAAGCTGGAATAATATTAGCCGGa                         >  2:458164/1‑93 (MQ=255)
                                                               cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                        >  1:1975538/1‑93 (MQ=255)
                                                                ttACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGGATAATCTTAGCCGTaaa                       >  1:527562/1‑93 (MQ=255)
                                                                ttACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                       >  2:3142556/1‑93 (MQ=255)
                                                                ttACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                       <  2:2788553/93‑1 (MQ=255)
                                                                  acaacaAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                     >  1:4085365/1‑93 (MQ=255)
                                                                    aacaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGCAAAACg                   >  2:235014/1‑93 (MQ=255)
                                                                    aacaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                   <  1:4424948/93‑1 (MQ=255)
                                                                    aacaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                   >  1:4330243/1‑93 (MQ=255)
                                                                    aaaaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                   <  1:3012164/90‑1 (MQ=255)
                                                                     aaaaacgGGGCGCTAAATGAGCTCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGa                  <  2:3447524/91‑1 (MQ=255)
                                                                     aaaaacgGGGCGCTAAATGAGCTCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGa                  <  1:4172137/91‑1 (MQ=255)
                                                                      caaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                 >  1:2389164/1‑93 (MQ=255)
                                                                      aaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                 <  1:2233497/92‑1 (MQ=255)
                                                                       aaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                >  1:3123141/1‑93 (MQ=255)
                                                                        agcggggCACTCAGGGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGAAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg               <  2:4215167/89‑1 (MQ=255)
                                                                        aacaGGGCGCTCAATGAGCGCCCTGTCTCTTTTTTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAAGAATg               >  1:1959440/1‑93 (MQ=255)
                                                                         acggggCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt              <  1:309447/90‑1 (MQ=255)
                                                                          caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATATc             >  2:4006359/1‑93 (MQ=255)
                                                                          caGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc             <  1:2694478/93‑1 (MQ=255)
                                                                          caGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAAGAATGTc             >  1:4545064/1‑93 (MQ=255)
                                                                           aGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCa            >  1:2909787/1‑93 (MQ=255)
                                                                            gggCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa           >  2:1052547/1‑93 (MQ=255)
                                                                              gcgcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa         <  2:2224376/93‑1 (MQ=255)
                                                                                gcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg       >  1:1038196/1‑93 (MQ=255)
                                                                                gcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg       >  2:3801448/1‑93 (MQ=255)
                                                                                gcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg       >  1:1320567/1‑93 (MQ=255)
                                                                                gcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg       >  2:1592270/1‑93 (MQ=255)
                                                                                gcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg       <  1:1512010/93‑1 (MQ=255)
                                                                                  tCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAggga     >  1:945388/1‑91 (MQ=255)
                                                                                     aTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg  >  2:1115223/1‑93 (MQ=255)
                                                                                     aTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg  >  2:2916311/1‑93 (MQ=255)
                                                                                          |                                                                                        
GAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAG  >  NC_000913/4225049‑4225226

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