Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,085,675 Δ1 bp intergenic (+38/+13) yiiG → / ← frvR hypothetical protein/putative frv operon regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,085,6750G.99.4% 2148.6 / 2.2 483intergenic (+38/+13)yiiG/frvRhypothetical protein/putative frv operon regulatory protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base . (291/188);  minor base C (0/3);  total (292/191)
Fisher's exact test for biased strand distribution p-value = 6.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 526 total aligned reads displayed.

GCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACA  >  NC_000913/4085587‑4085762
                                                                                        |                                                                                       
 cATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGa                                                                                   >  1:2827645/1‑93 (MQ=255)
   ttACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                 <  1:639687/93‑1 (MQ=255)
    tACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGacac                                                                                >  1:7660321/1‑93 (MQ=255)
     aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGAcaca                                                                               <  1:7592457/93‑1 (MQ=255)
     aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGAcaca                                                                               <  1:4716805/93‑1 (MQ=255)
      cgcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAg                                                                              >  2:272018/1‑93 (MQ=255)
       gcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGt                                                                             >  1:7387469/1‑93 (MQ=255)
        cgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGtt                                                                            >  2:8846437/1‑93 (MQ=255)
         gAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCTGGTGGGGCGGATGCGACACAGGTc                                                                           >  2:5429576/1‑93 (MQ=255)
         gAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTc                                                                           >  1:2450025/1‑93 (MQ=255)
          aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCa                                                                          >  1:9993396/1‑93 (MQ=255)
          aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCa                                                                          <  1:6694557/93‑1 (MQ=255)
          aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCa                                                                          >  1:1117540/1‑93 (MQ=255)
            tACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcc                                                                        <  2:9685492/93‑1 (MQ=255)
            tACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcc                                                                        >  2:8073614/1‑93 (MQ=255)
             aCAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccc                                                                       >  1:6291509/1‑93 (MQ=255)
              cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcccc                                                                      <  1:5700389/93‑1 (MQ=255)
               aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccccc                                                                     <  2:1563835/93‑2 (MQ=255)
                  gAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCtgt                                                                  >  2:8308317/1‑93 (MQ=255)
                  aaaaTGTTTAATGACTATAATAGCTTGCTTTAACCTTTTTTGATGATAAATGCGGGGTCGGAGGGGGCCCACCGGACACAGTTCACCCCCtgt                                                                  <  2:9045758/92‑1 (MQ=255)
                   aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCtgtg                                                                 >  2:4656830/1‑93 (MQ=255)
                   aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCtgtg                                                                 >  2:2551649/1‑93 (MQ=255)
                     aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGaa                                                               <  1:6596186/93‑1 (MQ=255)
                     aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGaa                                                               <  1:8239039/93‑1 (MQ=255)
                     aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGaa                                                               >  2:3185424/1‑93 (MQ=255)
                      tGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAg                                                              >  1:10031584/1‑93 (MQ=255)
                       ggTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAga                                                             >  1:6138015/1‑92 (MQ=255)
                       ggTTGATGACTATAATAGCCTGCGTAAAACGTTTTTGATGATAAATGCAGTGTCGGAGGCGGCGCACCCGACACAGTTCACCCCCTGTGAAGt                                                             <  2:3984726/93‑1 (MQ=255)
                        gTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTa                                                            >  1:9338285/1‑93 (MQ=255)
                        cttGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTa                                                            <  1:1754949/92‑1 (MQ=255)
                         ttGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTaa                                                           >  2:6197378/1‑93 (MQ=255)
                         tggatgaCTATAAGAGCCTGCGTTAAACGTTTTTGATGAGAAATGCAGGGTCGGATGCGGCGCACCGGACACAGTTCACCCCCTGTGAAGTaa                                                           <  2:954267/91‑1 (MQ=255)
                           gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACt                                                         <  1:8146406/93‑1 (MQ=255)
                            atgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTc                                                        <  2:1772460/93‑1 (MQ=255)
                              gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCa                                                      >  1:3518590/1‑93 (MQ=255)
                              gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCa                                                      <  2:1315101/93‑1 (MQ=255)
                               aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCAt                                                     <  1:1533038/93‑1 (MQ=255)
                               aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCAt                                                     <  1:2626913/93‑1 (MQ=255)
                                cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATa                                                    >  1:45833/1‑93 (MQ=255)
                                 tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaa                                                   >  1:6877925/1‑93 (MQ=255)
                                 tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaa                                                   >  2:3410499/1‑93 (MQ=255)
                                  ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                  >  2:8683259/1‑93 (MQ=255)
                                  ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                  >  2:4780851/1‑93 (MQ=255)
                                   taataGCCTGCTTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCACCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAAt                                                 <  2:783798/93‑1 (MQ=255)
                                    aataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATg                                                >  2:6445962/1‑93 (MQ=255)
                                     ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACGGTTCACCCCCTGTGAAGTAACTCCATAAATGc                                               <  2:981697/93‑1 (MQ=255)
                                     ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGc                                               >  1:8500668/1‑93 (MQ=255)
                                     ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGc                                               <  1:4292824/93‑1 (MQ=255)
                                      tggCCTGCGTTAACCTTTTTTGATGATAAATGCAGTGTCGTATGCGGCGCATCCGCCACAGTTCACCCCCTGTGAAGTAACGCCATAAATGcc                                              <  1:2129703/91‑1 (MQ=255)
                                      taGCTTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGcc                                              <  2:7731994/93‑1 (MQ=255)
                                         ccTGCTTTAAACGTTTTTGATGATAAATGCATTGTCGGATGCGGCCCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCa                                           <  1:2045006/93‑1 (MQ=255)
                                           tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                         <  1:8340996/93‑1 (MQ=255)
                                           tGCGTTAAACGTTTTTGATAATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                         <  1:7910262/93‑1 (MQ=255)
                                            gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCatat                                        >  2:6371589/1‑93 (MQ=255)
                                            gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCatat                                        >  1:9645873/1‑93 (MQ=255)
                                            gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCatat                                        >  1:5522629/1‑93 (MQ=255)
                                             cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATc                                       <  1:1075484/93‑1 (MQ=255)
                                             cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATc                                       >  1:9254380/1‑93 (MQ=255)
                                              gTTAAACGTTTTTGATGATAAATGCAGTGTCTGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATcc                                      <  2:4417000/93‑1 (MQ=255)
                                               ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                     >  2:5769299/1‑93 (MQ=255)
                                               ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                     <  1:7020909/93‑1 (MQ=255)
                                                tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                    >  2:5140124/1‑93 (MQ=255)
                                                tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                    >  1:5044015/1‑93 (MQ=255)
                                                tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTAACCCCCTGTGAAGAAACTACATAAATGCCTCATATCCTg                                    >  2:8224619/1‑93 (MQ=255)
                                                  aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                  <  2:8556666/93‑1 (MQ=255)
                                                  aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                  >  2:5015223/1‑93 (MQ=255)
                                                   aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                 >  2:352310/1‑93 (MQ=255)
                                                     gTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc                               <  2:4328566/93‑1 (MQ=255)
                                                       ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAccc                             <  1:3617573/93‑1 (MQ=255)
                                                       ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAccc                             <  2:8464029/93‑1 (MQ=255)
                                                        tttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCg                            >  1:9842485/1‑93 (MQ=255)
                                                         ttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg                           >  2:9268417/1‑93 (MQ=255)
                                                          tgatgaTAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                          <  1:4895261/93‑1 (MQ=255)
                                                          tgatgaTAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                          >  1:8746402/1‑93 (MQ=255)
                                                          tgatgaTAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                          >  1:8934513/1‑93 (MQ=255)
                                                            atgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaa                        >  2:2326325/1‑93 (MQ=255)
                                                             tgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaat                       >  2:4205885/1‑93 (MQ=255)
                                                              gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaata                      >  2:5890371/1‑93 (MQ=255)
                                                              gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaata                      >  2:121034/1‑93 (MQ=255)
                                                               atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                     >  1:2889900/1‑93 (MQ=255)
                                                               atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                     >  2:2223195/1‑93 (MQ=255)
                                                                tAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat                    >  2:6841759/1‑93 (MQ=255)
                                                                   aTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtctc                 <  2:6048234/93‑1 (MQ=255)
                                                                    tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTCACCCGGCAATAATCTcc                >  2:4498643/1‑93 (MQ=255)
                                                                    tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                >  2:555625/1‑93 (MQ=255)
                                                                    tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                >  2:6612993/1‑93 (MQ=255)
                                                                     gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccc               >  2:7502305/1‑92 (MQ=255)
                                                                     gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTacg               >  2:3672148/1‑93 (MQ=255)
                                                                      cAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgc              >  1:4059758/1‑93 (MQ=255)
                                                                       aGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc             >  2:2993524/1‑93 (MQ=255)
                                                                            cGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTACGCCGGATg        >  2:6725551/1‑93 (MQ=255)
                                                                             ggATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc       <  1:3301181/93‑1 (MQ=255)
                                                                             ggATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc       <  1:3242958/93‑1 (MQ=255)
                                                                              gATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt      >  1:4228604/1‑93 (MQ=255)
                                                                              gATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt      >  2:3885917/1‑93 (MQ=255)
                                                                               aggcggcgCACCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg     <  2:7235057/91‑1 (MQ=255)
                                                                               aTGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg     >  1:3186555/1‑93 (MQ=255)
                                                                                tGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa    <  1:3997158/93‑1 (MQ=255)
                                                                                tGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa    <  1:4448023/93‑1 (MQ=255)
                                                                                 gcggcgCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac   >  1:222330/1‑93 (MQ=255)
                                                                                        |                                                                                       
GCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACA  >  NC_000913/4085587‑4085762

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