Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,584,031 A→G 75.3% intergenic (‑72/+176) safA ← / ← ydeP two‑component system connector membrane protein, EvgSA to PhoQP/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,584,0310AG75.3% 250.2 / 78.6 141intergenic (‑72/+176)safA/ydePtwo‑component system connector membrane protein, EvgSA to PhoQP/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (20/15);  new base G (50/56);  total (70/71)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 147 total aligned reads displayed.

TGTTTTTCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTTTAATTTAATCCA  >  NC_000913/1583935‑1584127
                                                                                                |                                                                                                
   ttttCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAgg                                                                                            >  1:3746790/1‑100 (MQ=255)
     ttCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAAGGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGaa                                                                                          >  2:2306731/1‑100 (MQ=255)
      tCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGaaa                                                                                         <  1:2619687/100‑1 (MQ=255)
       cACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAAt                                                                                        >  2:1641127/1‑100 (MQ=255)
       cACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAAt                                                                                        >  2:3282467/1‑100 (MQ=255)
         ntgtgGTCGCCTGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGATGGAGAGAAGAGGAAATAg                                                                                      >  2:3334051/2‑100 (MQ=255)
          tgtgGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAgg                                                                                     <  1:3257724/100‑1 (MQ=255)
           gtgGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGa                                                                                    >  1:3559924/1‑100 (MQ=255)
           gtgGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGa                                                                                    >  1:3583854/1‑100 (MQ=255)
             ggTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACt                                                                                  >  1:3177688/1‑100 (MQ=255)
                  catgcatTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGc                                                                             >  1:27464/1‑100 (MQ=255)
                   atgcatTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCt                                                                            <  1:848134/100‑1 (MQ=255)
                    tgcatTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCtt                                                                           >  1:2896067/1‑100 (MQ=255)
                     gcatTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTc                                                                          <  1:3090199/100‑1 (MQ=255)
                       atTTCAAATATGTTTATTTATCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAATAGGAAATAGGAATGGTGCTTCa                                                                         >  2:2037785/1‑99 (MQ=255)
                       atTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAg                                                                        >  2:3130714/1‑100 (MQ=255)
                       atTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAg                                                                        <  1:2900982/100‑1 (MQ=255)
                       atTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAg                                                                        <  1:1033568/100‑1 (MQ=255)
                         ttCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAg                                                                        <  2:425920/98‑1 (MQ=255)
                             aaTATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGa                                                                  >  1:3068304/1‑100 (MQ=255)
                              atatGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGTAAATAGGACTGGTGCTTCAGTACTGa                                                                  <  2:70137/99‑1 (MQ=255)
                               tatGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGc                                                                >  1:1969575/1‑100 (MQ=255)
                               tatGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGc                                                                <  1:178496/100‑1 (MQ=255)
                                atgatcaattaGCGGATAACGTTAAAAATCGTTAAGGAGTATGCGTAATATAACAGGATGGAGAGAAGGGGAAATAGGACTGGTGCTTCAGTACTGAGc                                                                <  2:1177286/91‑1 (MQ=255)
                                atGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGCAGAGGAAATAGGACTGGTGCTTCAGTACTGa                                                                  >  2:2500239/1‑97 (MQ=255)
                                 tGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGAAAATAGGACTAGTGCTTCAGTACTGAGCgg                                                              >  2:3619517/1‑100 (MQ=255)
                                 tGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCgg                                                              <  1:3809117/100‑1 (MQ=255)
                                 tGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCgg                                                              <  2:89819/100‑1 (MQ=255)
                                  gTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGGACTGAGCGGa                                                             >  2:3898334/1‑100 (MQ=255)
                                  gTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTGAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTAGTGCTTCAGTACTGAGCGGa                                                             >  2:643888/1‑100 (MQ=255)
                                   tttatttaGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAg                                                            <  1:2603665/100‑1 (MQ=255)
                                    ttatttaGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGt                                                           <  1:2542293/100‑1 (MQ=255)
                                     tatttaGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGtt                                                          <  1:2241639/100‑1 (MQ=255)
                                        ttaGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCt                                                       <  1:2444523/100‑1 (MQ=255)
                                         taGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCtt                                                      <  1:3883606/100‑1 (MQ=255)
                                         taGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCtt                                                      <  2:3787442/100‑1 (MQ=255)
                                           gCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTa                                                     <  1:3282490/99‑1 (MQ=255)
                                              gATAACGTTAAAAATAGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGc                                                 <  2:3415440/100‑1 (MQ=255)
                                                 aaCGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGt                                              <  2:1049831/100‑1 (MQ=255)
                                                  aCGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTa                                             <  1:1655482/100‑1 (MQ=255)
                                                  aCGTTAAAAATCGTTAATAAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTa                                             <  2:961223/100‑1 (MQ=255)
                                                   cGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAg                                            >  2:3178540/1‑100 (MQ=255)
                                                   cGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAg                                            >  1:2780345/1‑100 (MQ=255)
                                                   cGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAg                                            >  2:3488251/1‑100 (MQ=255)
                                                   cGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAg                                            <  1:1741776/100‑1 (MQ=255)
                                                    gTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAgg                                           >  1:3698783/1‑100 (MQ=255)
                                                    gTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAgg                                           >  2:430431/1‑100 (MQ=255)
                                                     ttAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGc                                          >  1:3684891/1‑100 (MQ=255)
                                                     ttAAAAATCGTTAATAAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACTGCTGTa                                             >  1:388948/1‑97 (MQ=255)
                                                      tAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCCGTACTGAGCGGAGTTTCTTACAGCTGTAGGCa                                         >  2:3507432/1‑100 (MQ=255)
                                                       aaaaaTCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAggaa                                         <  2:1941013/99‑3 (MQ=255)
                                                       aaaaaTCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAg                                        <  1:1097510/100‑1 (MQ=255)
                                                           aTCGTTAAGCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGt                                    <  1:147344/100‑1 (MQ=255)
                                                              gTTAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTta                                 <  2:3741628/100‑1 (MQ=255)
                                                                tAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTtata                               <  1:738075/100‑1 (MQ=255)
                                                                tAATCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTtata                               <  2:2351471/100‑1 (MQ=255)
                                                                 aaTCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTatat                              <  2:1233589/100‑1 (MQ=255)
                                                                  aTCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATc                             >  1:823693/1‑100 (MQ=255)
                                                                   tCAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCt                            <  1:14719/100‑1 (MQ=255)
                                                                   tCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATc                             <  2:1912086/99‑1 (MQ=255)
                                                                   tCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCt                            >  2:2064290/1‑100 (MQ=255)
                                                                    cAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGCCTGGTGCTTCAGTATTGAACGGAGATTATTACAGCTGTAGACAGAAGTTTATATc                             >  2:2151977/1‑98 (MQ=255)
                                                                    cAGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCtt                           >  1:1874708/1‑100 (MQ=255)
                                                                     aGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTta                          <  1:2836698/100‑1 (MQ=255)
                                                                     aGTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTta                          >  2:4076810/1‑100 (MQ=255)
                                                                      gTATGCTTAATATAACAGGCTGGAGAGAAGATGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTgat                         <  2:176570/100‑1 (MQ=255)
                                                                      gTATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTtat                         <  1:3178560/100‑1 (MQ=255)
                                                                      gTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTtat                         >  2:3608337/1‑100 (MQ=255)
                                                                       tATGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTtata                        <  1:3501895/100‑1 (MQ=255)
                                                                        aTGCTTAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAg                       <  1:3297949/100‑1 (MQ=255)
                                                                          gCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTa                     >  2:1418852/1‑100 (MQ=255)
                                                                            ttAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTTGTGCTTCAGTACTTAGCGTAGTTTCTTACAGCTTTAGGCAGAAGTTTATATCTTATAGTATg                   <  2:2742270/100‑1 (MQ=255)
                                                                             tAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGt                  <  2:151949/100‑1 (MQ=255)
                                                                             tAATATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGt                  >  2:2012616/1‑100 (MQ=255)
                                                                              aaTATAACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTa                 >  1:1011376/1‑100 (MQ=255)
                                                                              aaTATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGt                  >  1:2737660/1‑99 (MQ=255)
                                                                                 ataACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtt              >  1:1249986/1‑100 (MQ=255)
                                                                                 ataACAGGCTGGAGAAAAGAGGAAATAGTACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAg                <  1:1156705/98‑1 (MQ=255)
                                                                                  taACAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTGCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGt               >  2:336667/1‑98 (MQ=255)
                                                                                  taACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGttt             >  2:1275965/1‑100 (MQ=255)
                                                                                   aaCAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGttta            >  2:2786121/1‑100 (MQ=255)
                                                                                   aaCAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGttta            <  2:1402494/100‑1 (MQ=255)
                                                                                    aCAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaa           >  1:1367837/1‑100 (MQ=255)
                                                                                     cAGGCTGGAGGGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaat          <  2:1175765/100‑1 (MQ=255)
                                                                                     cAGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaat          <  1:100810/100‑1 (MQ=255)
                                                                                      aGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaatt         <  2:1046159/100‑1 (MQ=255)
                                                                                      aGGCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaatt         <  1:2838237/100‑1 (MQ=255)
                                                                                       ggCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaattt        <  1:196543/100‑1 (MQ=255)
                                                                                       ggCTGGAGAGAAAAGGAAATAGGACTATTGCGTCAGTACTGAGTGGAGTATCTTACAGATTTAGGCAGAAGTTTATATCTTATAGTTTGTAGtttaattt        <  2:1865573/100‑1 (MQ=255)
                                                                                       ggCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACGGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaattt        >  1:2297330/1‑100 (MQ=255)
                                                                                       ggCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATCGTATGTAGtttaattt        >  2:2437865/1‑100 (MQ=255)
                                                                                       ggCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaattt        >  1:604544/1‑100 (MQ=255)
                                                                                        gCTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaattta       >  2:293989/1‑100 (MQ=255)
                                                                                         cTGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGtttaatttaa      >  1:4058123/1‑100 (MQ=255)
                                                                                          tGGAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTttaatttaat     <  1:3517582/100‑1 (MQ=255)
                                                                                           ggAGAGAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTTTAATTTAATc    <  2:4163482/100‑1 (MQ=255)
                                                                                            gagaAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTTTAATTTAATcc   <  1:836145/100‑1 (MQ=255)
                                                                                             agagAAGAGGAAATAGGACTGGTGCTTCAGTACTTAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTTTAATTTAATCCa  <  1:3983579/100‑1 (MQ=255)
                                                                                             agagAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTTTAATTTAATCCa  >  2:2230499/1‑100 (MQ=255)
                                                                                             agagAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTATTACAGCTGTATGCAGAAGTTTATATCTTAGAGTATGTAGTttaatttaat     <  2:2552240/97‑1 (MQ=255)
                                                                                                |                                                                                                
TGTTTTTCACTGTGGTCGCATGCATTTCAAATATGTTTATTTAGCGGATAACGTTAAAAATCGTTAATCAGTATGCTTAATATAACAGGCTGGAGAAAAGAGGAAATAGGACTGGTGCTTCAGTACTGAGCGGAGTTTCTTACAGCTGTAGGCAGAAGTTTATATCTTATAGTATGTAGTTTAATTTAATCCA  >  NC_000913/1583935‑1584127

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