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.7% 1457.1 / ‑3.7 324intergenic (+38/+13)yiiG/frvRhypothetical protein/putative frv operon regulatory protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (196/127);  minor base C (0/1);  total (196/128)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 347 total aligned reads displayed.

GCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACA  >  NC_000913/4085587‑4085762
                                                                                        |                                                                                       
  aTTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGac                                                                                  <  2:1113585/93‑1 (MQ=255)
   ttACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                 <  1:1109137/93‑1 (MQ=255)
   ttACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                 <  1:7868132/93‑1 (MQ=255)
   ttACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                 >  2:1279444/1‑93 (MQ=255)
   ttACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGCTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                 >  1:752919/1‑93 (MQ=255)
         gAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTc                                                                           <  2:6081896/93‑1 (MQ=255)
          aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCa                                                                          >  2:4968471/1‑93 (MQ=255)
          aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCa                                                                          >  1:4529076/1‑93 (MQ=255)
          aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCa                                                                          <  1:3795836/93‑1 (MQ=255)
           gTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAc                                                                         >  1:5124352/1‑93 (MQ=255)
             aCAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccc                                                                       >  2:1804347/1‑93 (MQ=255)
               aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccccc                                                                     <  2:1819741/93‑2 (MQ=255)
                aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGTGGCGCATCCGACACAGTTCAccccct                                                                    <  1:1030990/93‑3 (MQ=255)
                  gAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCtgt                                                                  >  2:4047279/1‑93 (MQ=255)
                  gAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCtgt                                                                  <  1:487624/93‑1 (MQ=255)
                   aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCtgtg                                                                 >  2:7780964/1‑93 (MQ=255)
                    aaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGa                                                                >  2:1935337/1‑93 (MQ=255)
                    aaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGa                                                                >  1:2439035/1‑93 (MQ=255)
                      tGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAg                                                              >  2:5736796/1‑93 (MQ=255)
                       ggTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGt                                                             <  2:7649061/93‑1 (MQ=255)
                        gTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTa                                                            >  1:3694171/1‑93 (MQ=255)
                         ttGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTaa                                                           >  1:405285/1‑93 (MQ=255)
                         ttGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGACAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTaa                                                           >  1:2811800/1‑93 (MQ=255)
                           gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACt                                                         <  1:2498566/93‑1 (MQ=255)
                            atgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTc                                                        <  2:1038177/93‑1 (MQ=255)
                            atgaCTATAATAGCCTGCGTTAAACGCTTTTGATGATAAATGCAGTGTCGGATTCGGCGCATCCGACACAGTTCACCTCCTGTGAAGCAACTc                                                        >  1:7930524/1‑93 (MQ=255)
                             tgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTcc                                                       >  2:6431738/1‑93 (MQ=255)
                                cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTTACCCCCTGTGAAGTAACTCCATa                                                    >  2:5371022/1‑93 (MQ=255)
                                cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATa                                                    <  1:2731155/93‑1 (MQ=255)
                                cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATa                                                    >  1:873788/1‑93 (MQ=255)
                                cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATa                                                    >  1:1439190/1‑93 (MQ=255)
                                 tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaa                                                   >  2:2799475/1‑93 (MQ=255)
                                  ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                  >  1:2359441/1‑93 (MQ=255)
                                  ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                  >  1:3960820/1‑93 (MQ=255)
                                  ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                  <  1:5484589/93‑1 (MQ=255)
                                    aataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATg                                                >  2:1616271/1‑93 (MQ=255)
                                     ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGc                                               >  1:3041308/1‑93 (MQ=255)
                                      taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGcc                                              <  1:2606981/93‑1 (MQ=255)
                                        gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTc                                            <  1:6332845/93‑1 (MQ=255)
                                         ccTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCa                                           >  2:4517312/1‑93 (MQ=255)
                                         ccTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCa                                           <  1:5024694/93‑1 (MQ=255)
                                           tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                         >  2:1987633/1‑93 (MQ=255)
                                           tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                         >  2:6028934/1‑93 (MQ=255)
                                           tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                         <  2:7897143/93‑1 (MQ=255)
                                           tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                         >  1:1826251/1‑93 (MQ=255)
                                            gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCatat                                        >  1:2018702/1‑93 (MQ=255)
                                             cGTTAACCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATc                                       <  1:2188137/93‑1 (MQ=255)
                                             cGTTAACCGTTTTTGATGATAAATGCAGTGTCGGAGGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATc                                       <  1:4235608/93‑1 (MQ=255)
                                             cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATc                                       >  1:5799757/1‑93 (MQ=255)
                                              gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATcc                                      >  1:1929192/1‑93 (MQ=255)
                                              gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATcc                                      >  1:1256700/1‑93 (MQ=255)
                                              gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATcc                                      >  2:3702861/1‑93 (MQ=255)
                                               ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                     >  1:7649613/1‑93 (MQ=255)
                                               ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                     <  2:4167218/93‑1 (MQ=255)
                                               ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                     >  2:7722818/1‑93 (MQ=255)
                                               ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                     <  2:5957135/93‑1 (MQ=255)
                                                tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCGCCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                    >  1:4448401/1‑93 (MQ=255)
                                                tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                    >  1:7708639/1‑93 (MQ=255)
                                                tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                    <  1:7574022/93‑1 (MQ=255)
                                                tAAACGGTTTTGATGATAAATGCAGTGTCGGATTCGGCGCATCCGACACAGTTCACCCCCTGTGAAGGAACTCCATAAATGCCTCATATCCTg                                    >  1:4764097/1‑93 (MQ=255)
                                                  aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                  >  1:434213/1‑93 (MQ=255)
                                                   aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                 >  2:67390/1‑93 (MQ=255)
                                                   aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                 <  1:6553502/93‑1 (MQ=255)
                                                    cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCTCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTaa                                >  1:3207065/1‑93 (MQ=255)
                                                     gTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc                               >  1:6869858/1‑93 (MQ=255)
                                                     gTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc                               >  2:3950045/1‑93 (MQ=255)
                                                      tttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATCTCCTGTTAAcc                              >  1:6773678/1‑93 (MQ=255)
                                                       ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAccc                             >  2:508468/1‑93 (MQ=255)
                                                            atgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaa                        >  2:4748510/1‑93 (MQ=255)
                                                             tgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaat                       >  2:1538536/1‑93 (MQ=255)
                                                             tgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaat                       <  2:2254282/93‑1 (MQ=255)
                                                             tgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaat                       >  1:5998352/1‑93 (MQ=255)
                                                              gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaata                      >  2:58312/1‑93 (MQ=255)
                                                              gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaata                      <  1:2134724/93‑1 (MQ=255)
                                                               atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                     >  2:854612/1‑93 (MQ=255)
                                                               atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                     >  1:6526714/1‑93 (MQ=255)
                                                               atAAAGGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                     <  2:3195648/93‑1 (MQ=255)
                                                                tAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat                    <  1:2537533/93‑1 (MQ=255)
                                                                 aaaTGCAGTGTCGGATGCGGCGCCTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtc                   <  2:4329977/93‑1 (MQ=255)
                                                                 aaaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtc                   >  1:3511172/1‑93 (MQ=255)
                                                                   aTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtctc                 >  1:7207055/1‑93 (MQ=255)
                                                                    tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                >  1:7100274/1‑93 (MQ=255)
                                                                     gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg               >  2:3511701/1‑93 (MQ=255)
                                                                     gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg               >  2:730425/1‑93 (MQ=255)
                                                                       aGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc             >  2:6695950/1‑93 (MQ=255)
                                                                        gtgcCGGAGCCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgccg            <  2:7087806/93‑1 (MQ=255)
                                                                         tgtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCgg           >  2:4613855/1‑93 (MQ=255)
                                                                          gtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa          <  1:4793100/93‑1 (MQ=255)
                                                                          gtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa          >  2:1381239/1‑93 (MQ=255)
                                                                          gtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa          >  2:2403584/1‑93 (MQ=255)
                                                                           tCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt         <  1:3376517/93‑1 (MQ=255)
                                                                             ggATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc       >  2:7458221/1‑93 (MQ=255)
                                                                             ggAGGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc       <  2:3011337/93‑1 (MQ=255)
                                                                              gATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt      >  2:6357723/1‑93 (MQ=255)
                                                                              gATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt      >  2:2696410/1‑93 (MQ=255)
                                                                              gATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt      >  1:5307613/1‑93 (MQ=255)
                                                                                 gcggcgCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac   >  2:5921495/1‑93 (MQ=255)
                                                                                 gcggcgCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac   >  2:509279/1‑93 (MQ=255)
                                                                                 gcggcgCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac   <  1:6229930/93‑1 (MQ=255)
                                                                                  cggcgCATCCCACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGaca  <  2:2535473/93‑1 (MQ=255)
                                                                                        |                                                                                       
GCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACA  >  NC_000913/4085587‑4085762

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