Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,088,091 G→C 15.6% intergenic (+34/+16) yiiG → / ← frvR DUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,088,0910GC15.6% 377.4 / 58.7 134intergenic (+34/+16)yiiG/frvRDUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain
Reads supporting (aligned to +/- strand):  ref base G (61/52);  new base C (11/10);  total (72/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 150 total aligned reads displayed.

GACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGC  >  NC_000913/4087992‑4088189
                                                                                                   |                                                                                                  
gACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcgg                                                                                                    <  1:1787268/100‑1 (MQ=255)
gACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcgg                                                                                                    <  1:1655950/100‑1 (MQ=255)
  ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGCTgcggcg                                                                                                  >  1:855626/1‑100 (MQ=255)
  ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggcg                                                                                                  <  1:3205965/100‑1 (MQ=255)
  ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATAGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcgccg                                                                                                  <  1:1398382/100‑1 (MQ=255)
   tctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGgcgc                                                                                                 >  1:1016420/1‑100 (MQ=255)
     tttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCAt                                                                                               >  2:576417/1‑100 (MQ=255)
      ttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATc                                                                                              <  2:126338/100‑1 (MQ=255)
      ttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATc                                                                                              <  2:1501051/100‑1 (MQ=255)
       tCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATcc                                                                                             <  1:2875485/100‑1 (MQ=255)
        cccGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTATACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCg                                                                                            >  1:2405637/1‑100 (MQ=255)
        cccGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCg                                                                                            >  1:2931226/1‑100 (MQ=255)
        cccGAGAGCATTACGCGAGTACAACGAAATGGTTCATGACTATAATAGCCTGCGTTAAACGTTTTTGATGAGAAATGCAGTGTCGGATCAGGCGCATcc                                                                                             >  1:290826/1‑99 (MQ=255)
         ccGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGa                                                                                           <  2:2007164/100‑1 (MQ=255)
          cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGac                                                                                          >  2:1309023/1‑100 (MQ=255)
          cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                          <  2:79204/100‑1 (MQ=255)
           gagaGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                         <  2:2139623/100‑1 (MQ=255)
            agagCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGacac                                                                                        <  2:2036484/100‑1 (MQ=255)
            agagCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacac                                                                                        <  2:2954464/100‑1 (MQ=255)
            agagCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacac                                                                                        <  2:2597109/100‑1 (MQ=255)
              agCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAg                                                                                      <  1:3319662/100‑1 (MQ=255)
               gCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacactgc                                                                                     <  1:1935348/100‑4 (MQ=255)
                cATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGtt                                                                                    >  2:1156421/1‑100 (MQ=255)
                cATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGtt                                                                                    >  2:102516/1‑100 (MQ=255)
                  ttACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCa                                                                                  >  2:2097980/1‑100 (MQ=255)
                   tACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAc                                                                                 <  1:2459781/100‑1 (MQ=255)
                     cgcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGGCACAGTTCAcgc                                                                               >  1:1814198/1‑98 (MQ=255)
                     cgcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccc                                                                               >  2:3340765/1‑100 (MQ=255)
                     cgcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccc                                                                               >  2:2620166/1‑100 (MQ=255)
                      gcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcccc                                                                              <  2:2561326/100‑1 (MQ=255)
                          gTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtgt                                                                          >  2:740728/1‑100 (MQ=255)
                            aCAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGa                                                                        >  2:1452203/1‑100 (MQ=255)
                             cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                       >  1:254012/1‑100 (MQ=255)
                             cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                       >  2:3330987/1‑100 (MQ=255)
                             cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                       >  1:2090371/1‑100 (MQ=255)
                             cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                       >  1:226921/1‑100 (MQ=255)
                              aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAg                                                                      >  1:3051647/1‑100 (MQ=255)
                              aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAg                                                                      <  2:2525855/100‑1 (MQ=255)
                              aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAg                                                                      >  1:516695/1‑100 (MQ=255)
                              aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAg                                                                      <  2:748599/100‑1 (MQ=255)
                               aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGt                                                                     <  1:3102268/100‑1 (MQ=255)
                               aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGt                                                                     >  1:2899923/1‑100 (MQ=255)
                               aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGt                                                                     <  2:1678873/100‑1 (MQ=255)
                               aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGt                                                                     >  2:1171711/1‑100 (MQ=255)
                                  aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGa                                                                        <  1:1770833/94‑1 (MQ=255)
                                  aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAAc                                                                  <  1:3363769/100‑1 (MQ=255)
                                  aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAAc                                                                  <  2:523056/100‑1 (MQ=255)
                                    aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTc                                                                <  1:1045107/100‑1 (MQ=255)
                                     tGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTcc                                                               <  1:985151/100‑1 (MQ=255)
                                        ttGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATa                                                            <  2:3055568/100‑1 (MQ=255)
                                         tgatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaa                                                           <  2:1928554/100‑1 (MQ=255)
                                          gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaaa                                                          >  1:2122433/1‑100 (MQ=255)
                                          gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaaa                                                          >  1:1843988/1‑100 (MQ=255)
                                          gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaaa                                                          >  2:263936/1‑100 (MQ=255)
                                           atgaCTATAATAGCCTGCGTTAAATGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAAt                                                         <  2:2513458/100‑1 (MQ=255)
                                            tgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATg                                                        >  2:37209/1‑100 (MQ=255)
                                              atcttaataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                      <  1:2566244/96‑1 (MQ=255)
                                                tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTc                                                    <  2:3005599/100‑1 (MQ=255)
                                                     taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATc                                               >  1:1624130/1‑100 (MQ=255)
                                                       gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCt                                             <  2:357472/100‑1 (MQ=255)
                                                         cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGt                                           >  1:816538/1‑100 (MQ=255)
                                                          tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGt                                           <  1:2389374/99‑1 (MQ=255)
                                                            cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCACTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTaa                                        >  1:847922/1‑100 (MQ=255)
                                                            cGTAAAACGATTTTGATGTGAAATGCAGTGTCGGATGCGCCGCATCTGACACAGATGACCCATTGTCAAGTAACGCGATAAATGCCTCGTATCCTGTTaa                                        >  1:326609/1‑100 (MQ=255)
                                                             gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc                                       >  2:247315/1‑100 (MQ=255)
                                                              ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAcc                                      <  1:1211046/100‑1 (MQ=255)
                                                               tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAccc                                     >  1:2471304/1‑100 (MQ=255)
                                                               tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAccc                                     >  2:3168844/1‑100 (MQ=255)
                                                                aaaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCg                                    <  1:875947/100‑1 (MQ=255)
                                                                aaaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCg                                    >  2:2559052/1‑100 (MQ=255)
                                                                  aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                  >  2:3070694/1‑100 (MQ=255)
                                                                    gTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaa                                <  1:1460671/100‑1 (MQ=255)
                                                                     tttttGATGATAAATGCAGTGTCGGATGCGGCGCCTCCGACACAGTTCACCCCTTGTGAAGTTACTTCATAAATGCCTCATATCCTGTTAAcccn                                    <  1:3182997/95‑2 (MQ=255)
                                                                        ttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat                            >  1:2205503/1‑100 (MQ=255)
                                                                         tgatgaTAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtc                           >  2:410091/1‑100 (MQ=255)
                                                                         tgatgaTAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtc                           <  2:641404/100‑1 (MQ=255)
                                                                         tgatgaTAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtc                           >  1:1172568/1‑100 (MQ=255)
                                                                             gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg                       >  2:1758255/1‑100 (MQ=255)
                                                                                      tgtgtCGGATGCGCCGAATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATAACCTGTTAACCCGGCAATAATCTCCGCCGGATg                >  1:2762733/2‑98 (MQ=255)
                                                                                       gtgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg             >  2:2458048/1‑100 (MQ=255)
                                                                                       gtgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg             >  2:2785729/1‑100 (MQ=255)
                                                                                        tgtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa            <  2:2090382/100‑1 (MQ=255)
                                                                                        tgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa            >  2:2584543/1‑100 (MQ=255)
                                                                                         gtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac           <  1:423081/100‑1 (MQ=255)
                                                                                         gtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCCGTTAACCCGGCAATAATCTCCGCCGGATGCTGac           >  1:1537205/1‑100 (MQ=255)
                                                                                         gtCGGATGCGCCGCATCCGACACTGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac           >  1:2731259/1‑100 (MQ=255)
                                                                                         gtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac           >  1:585238/1‑100 (MQ=255)
                                                                                           cGGATGCGCCGCATCCGACACTGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGacac         >  2:547086/1‑100 (MQ=255)
                                                                                             gATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACaa       >  2:3175317/1‑100 (MQ=255)
                                                                                             gATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACaa       >  1:2005516/1‑100 (MQ=255)
                                                                                              aTGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAat      <  2:348449/100‑1 (MQ=255)
                                                                                               tGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAata     >  2:28046/1‑100 (MQ=255)
                                                                                                gcggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAatat    <  2:876124/100‑1 (MQ=255)
                                                                                                gcggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAatat    <  1:1309018/100‑1 (MQ=255)
                                                                                                 cggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg   >  1:1270970/1‑100 (MQ=255)
                                                                                                 cggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg   >  2:1762908/1‑100 (MQ=255)
                                                                                                 cggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg   >  1:873635/1‑100 (MQ=255)
                                                                                                 cggcgCATCCGACACAGTTCACAACTTGTGAAGTAACTCTATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg   >  1:1272592/1‑100 (MQ=255)
                                                                                                  ggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGc  <  2:2261383/100‑1 (MQ=255)
                                                                                                  ggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGc  >  1:3176390/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGC  >  NC_000913/4087992‑4088189

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