Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,124 +G intergenic (+58/+26) ilvC → / ← ppiC ketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,957,1231.G97.9% 1138.1 / 5.9 382intergenic (+57/+27)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (8/0);  new base G (106/267);  total (115/267)
Fisher's exact test for biased strand distribution p-value = 5.38e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 413 total aligned reads displayed.

GACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGA‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTC  >  NC_000913/3957033‑3957210
                                                                                                   |                                                                                       
gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                          <  2:6203091/93‑1 (MQ=21)
gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                          <  1:4340448/93‑1 (MQ=21)
gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                          <  1:4273713/93‑1 (MQ=21)
gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                          <  1:3610067/93‑1 (MQ=21)
     atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGATAACACTCAAggcct                                                                                    <  2:6454326/93‑5 (MQ=21)
     atatGAAACGTATTGCTGTTGCGAGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcct                                                                                    <  1:1522566/93‑5 (MQ=18)
      tatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCCCCC‑‑CTCTCCCCCAGGGGGAGGGGGGAAACACTCAAAGCCtt                                                                                    >  2:718416/1‑93 (MQ=18)
       atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                                                                   <  2:2198483/93‑1 (MQ=18)
       atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                                                                   <  1:4654929/93‑1 (MQ=18)
        tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                                                                  <  2:2343567/93‑1 (MQ=18)
        tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                                                                  <  2:4309308/93‑1 (MQ=18)
         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                                                                 <  1:6644371/93‑1 (MQ=18)
         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                                                                 <  2:5572490/93‑1 (MQ=18)
            aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                              <  1:2532422/93‑1 (MQ=18)
            aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                              <  1:7073373/93‑1 (MQ=18)
            aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                              <  2:1786916/93‑1 (MQ=18)
            aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCCCA‑‑CCCCCACCCCCTCCCCCCGCGAGAGGGAGAAAACACCCAACGCCTTCTCCTg                                                                              >  2:1020134/1‑93 (MQ=16)
              gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACCGGGGGAGGGGGAAAACACTCAAAGCCCTCTCCTCga                                                                            >  2:4839186/1‑93 (MQ=18)
              gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                            <  2:3229776/93‑1 (MQ=18)
                 ttGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                         <  1:7419598/93‑1 (MQ=18)
                 ttGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                         <  2:4064336/93‑1 (MQ=18)
                 ttGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                         <  2:1385545/93‑1 (MQ=21)
                 ttGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                         <  1:1042997/93‑1 (MQ=18)
                   gCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                                                       <  2:5822223/93‑1 (MQ=18)
                    cTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCC‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                      <  1:1601308/93‑1 (MQ=18)
                    cTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                      <  1:5983613/93‑1 (MQ=18)
                    cTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                      <  2:2890828/93‑1 (MQ=18)
                     tGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                                     <  2:591435/93‑1 (MQ=18)
                      gTTGCGGGTTAAGTGCGCGCTGATGCCCCCCCC‑‑CCCCCC‑‑CCCTCCCCCAGGGGGGGGGGGAAAAAACTCAAAGCCCTCTCCCTGGGAAgggcc                                                                    >  1:3721389/1‑90 (MQ=9)
                      gTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                    <  1:5504766/93‑1 (MQ=18)
                      gTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                    <  2:4762407/93‑1 (MQ=18)
                        tGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                                  <  1:2098444/93‑1 (MQ=18)
                        tGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                                  <  1:641766/93‑1 (MQ=18)
                         gCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                                                                 <  1:7229567/93‑1 (MQ=18)
                         gCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                                                                 <  2:5393144/93‑1 (MQ=18)
                           gggTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACCCTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                                                               <  2:2540807/93‑1 (MQ=18)
                            ggTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                                                              <  1:3090853/93‑1 (MQ=18)
                              ttAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                            <  2:2143146/93‑1 (MQ=18)
                               tAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                                                           <  1:1856043/93‑1 (MQ=18)
                               tAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                                                           <  1:1985038/93‑1 (MQ=18)
                                aaGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                                                          <  2:7498528/93‑1 (MQ=18)
                                  gTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGATAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                                                        <  1:5610471/93‑1 (MQ=18)
                                   tGCGCGCTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                       <  2:7322482/93‑1 (MQ=21)
                                       cgcTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTa                                                   <  1:4078766/93‑1 (MQ=21)
                                        gcTGATG‑‑CCCTCACCCCCCCC‑‑CTCTCCCCCAGGGGGGGGGGGAAAACACCCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAc                                                  >  1:6079129/1‑93 (MQ=16)
                                        gcTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAc                                                  <  1:668962/93‑1 (MQ=21)
                                         cTGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                                 <  2:301339/93‑1 (MQ=21)
                                          tGATG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                                                <  2:7450691/93‑1 (MQ=25)
                                            aTG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCa                                              <  1:1168425/93‑1 (MQ=25)
                                             tG‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAc                                             <  2:2769643/93‑1 (MQ=25)
                                              g‑‑CCCTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                                            <  1:5131444/93‑1 (MQ=25)
                                                 cccTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt                                           <  2:1860029/93‑1 (MQ=25)
                                                 cccTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt                                           <  2:6784023/93‑1 (MQ=25)
                                                  ccTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt                                          <  1:442348/93‑1 (MQ=25)
                                                   cTCA‑‑CCCCGACCCTCTCCCCCAGGGGGAGGGGGAAAAAACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTa                                         >  2:2480432/1‑93 (MQ=255)
                                                   cTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTa                                         <  1:4889566/93‑1 (MQ=37)
                                                   cTCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTa                                         >  2:3077036/1‑93 (MQ=37)
                                                    tCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa                                        >  2:1989154/1‑93 (MQ=37)
                                                    tCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa                                        <  2:5250433/93‑1 (MQ=37)
                                                    tCA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa                                        <  1:979458/93‑1 (MQ=37)
                                                     cA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                       <  1:2281938/93‑1 (MQ=37)
                                                     cA‑‑CCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                       >  1:1637760/1‑93 (MQ=37)
                                                              accctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGa                                >  2:1273191/4‑93 (MQ=255)
                                                              accctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGa                                <  2:7094255/90‑1 (MQ=255)
                                                               ccctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat                               >  2:3320665/3‑93 (MQ=255)
                                                               ccctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat                               <  2:3776040/91‑1 (MQ=255)
                                                                cctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGata                              <  1:684264/92‑1 (MQ=255)
                                                                 ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGatat                             >  2:5077752/1‑93 (MQ=255)
                                                                  cctcCCACAGGGAGAGGCGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGatat                             <  1:6172927/92‑1 (MQ=255)
                                                                   ctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc                           <  1:1009669/93‑1 (MQ=255)
                                                                   ctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc                           >  2:4204355/1‑93 (MQ=255)
                                                                   ctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc                           <  2:3451972/93‑1 (MQ=255)
                                                                    tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCg                          >  2:6923014/1‑93 (MQ=255)
                                                                    tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCg                          >  1:5220962/1‑93 (MQ=255)
                                                                     cccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGa                         <  1:754230/93‑1 (MQ=255)
                                                                      ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGaa                        >  2:3831131/1‑93 (MQ=255)
                                                                      ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGaa                        <  1:3524667/93‑1 (MQ=255)
                                                                       cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                       <  1:2351008/93‑1 (MQ=255)
                                                                       cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                       <  1:5290517/93‑1 (MQ=255)
                                                                       cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                       <  1:499455/93‑1 (MQ=255)
                                                                       cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                       >  2:4631926/1‑93 (MQ=255)
                                                                           gggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTggg                   >  2:5983903/1‑93 (MQ=255)
                                                                            ggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                  <  1:6560707/93‑1 (MQ=255)
                                                                            ggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                  <  1:584520/93‑1 (MQ=255)
                                                                              agagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtgt                <  2:3679227/93‑1 (MQ=255)
                                                                              agagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtgt                <  1:5734925/93‑1 (MQ=255)
                                                                              agagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtgt                <  1:3323768/93‑1 (MQ=255)
                                                                               gagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGtgtg               <  1:2635583/93‑1 (MQ=255)
                                                                                      gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg         <  2:6258644/93‑1 (MQ=255)
                                                                                       aaaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCggg        <  1:3423524/93‑1 (MQ=255)
                                                                                       aaaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCggg        >  1:734785/1‑93 (MQ=255)
                                                                                        aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc       <  1:1358480/93‑1 (MQ=255)
                                                                                        aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc       <  1:1303850/93‑1 (MQ=255)
                                                                                          acacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCg     <  2:4416793/93‑1 (MQ=255)
                                                                                           cacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCgg    <  2:2732323/93‑1 (MQ=255)
                                                                                           cacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCgg    <  2:3959475/93‑1 (MQ=255)
                                                                                           cacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCgg    <  1:6935897/93‑1 (MQ=255)
                                                                                           cacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCgg    <  1:2870514/93‑1 (MQ=255)
                                                                                            acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGt   >  2:1735684/1‑93 (MQ=255)
                                                                                            acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGt   <  2:6373364/93‑1 (MQ=255)
                                                                                                   |                                                                                       
GACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCGA‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTC  >  NC_000913/3957033‑3957210

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