Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 875,714 T→G 51.3% V127G (GTC→GGC)  yliF → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913875,7140TG51.3% 10.1 / 162.7 115V127G (GTC→GGC) yliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (32/21);  new base G (32/22);  total (64/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

Only 100 of 117 total aligned reads displayed.

TCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGG  >  NC_000913/875615‑875812
                                                                                                   |                                                                                                  
tCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGt                                                                                                    >  1:686813/1‑100 (MQ=255)
tCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGt                                                                                                    <  1:4030441/100‑1 (MQ=255)
   gTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCaa                                                                                                 >  2:3249357/1‑100 (MQ=255)
     ccTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCaaaa                                                                                               >  2:2643843/1‑100 (MQ=255)
        tATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATa                                                                                            >  2:306350/1‑100 (MQ=255)
        tATCTGATTTCAAGTAACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATa                                                                                            <  1:770897/100‑1 (MQ=255)
         aTCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATAc                                                                                           <  2:2664735/100‑1 (MQ=255)
         aTCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATAc                                                                                           >  2:1737002/1‑100 (MQ=255)
          tCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACa                                                                                          >  2:3145083/1‑100 (MQ=255)
            tGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAgg                                                                                        >  2:3725522/1‑100 (MQ=255)
             gatctCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAggg                                                                                       <  2:1735126/96‑1 (MQ=255)
             gATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAggg                                                                                       >  1:3359016/1‑100 (MQ=255)
             gATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAAGTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAggg                                                                                       >  2:1542307/1‑100 (MQ=255)
              aTTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGt                                                                                      <  1:2164843/100‑1 (MQ=255)
               tttCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGtt                                                                                     >  1:3608719/1‑100 (MQ=255)
               tttCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAAGTATTGTGCAGATATTATGCAGGTTGAGAAGTATCGACATGTCAAAAATACAGGGt                                                                                      >  2:3320441/1‑99 (MQ=255)
                  ccAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGCTATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTAt                                                                                  >  2:806492/1‑100 (MQ=255)
                  ccAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAgggagatc                                                                                  >  1:1735138/1‑95 (MQ=255)
                  ccAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTAt                                                                                  <  1:3145102/100‑1 (MQ=255)
                   cAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGACAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTAtc                                                                                 >  2:539656/1‑100 (MQ=255)
                    aGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTAtct                                                                                >  2:2170204/1‑100 (MQ=255)
                      tGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATctct                                                                              <  1:1644435/100‑1 (MQ=255)
                         cGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCtttt                                                                           <  1:2127084/100‑1 (MQ=255)
                         cGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCtttt                                                                           >  1:1159100/1‑100 (MQ=255)
                         cGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCtttt                                                                           >  1:993738/1‑100 (MQ=255)
                         cGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCtttt                                                                           >  2:3967797/1‑100 (MQ=255)
                          gAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTg                                                                          <  1:1232018/100‑1 (MQ=255)
                           aaTGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGa                                                                         >  1:3628394/1‑100 (MQ=255)
                            aTGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGAc                                                                        <  1:1737014/100‑1 (MQ=255)
                            aTGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGAc                                                                        <  1:2012342/100‑1 (MQ=255)
                            aTGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTAAAAAAAAAAAAAAAATACAGGGTTTATCTCTTTTGAc                                                                        <  2:3567354/100‑1 (MQ=255)
                             tGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACg                                                                       <  2:570573/100‑1 (MQ=255)
                             tGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACg                                                                       <  1:3678610/100‑1 (MQ=255)
                              gTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACgg                                                                      <  2:2094603/100‑1 (MQ=255)
                              gTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACgg                                                                      <  1:283536/100‑1 (MQ=255)
                              gTATTCAGGTAAAGGAGTTCGACAAAAATAATTGTGCAGATATCATGCAGGTTGTGAAGTAAAAAAAAAAAAAAAATACAGGGTTTATCTCTTTTGACgg                                                                      <  2:2643352/100‑1 (MQ=255)
                               tATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGt                                                                     >  2:3074710/1‑100 (MQ=255)
                               tATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGt                                                                     >  1:4128510/1‑100 (MQ=255)
                                aTTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTa                                                                    >  1:3172652/1‑100 (MQ=255)
                                aTTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTa                                                                    >  2:591740/1‑100 (MQ=255)
                                aTTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTa                                                                    >  1:615401/1‑100 (MQ=255)
                                     ggTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAAcc                                                               >  1:3564803/1‑100 (MQ=255)
                                     ggTAAAGGAGTTCGACAAAAATTATTGTGAAGATATTATGCAGGTTGTGAAGTAAAAAAAAAAAAAAAATACAGGGTTTATCTCTTTTGACGGTAAAAcc                                                               <  2:3296592/100‑1 (MQ=255)
                                      gTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCt                                                              >  2:3194240/1‑100 (MQ=255)
                                         aaGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCg                                                           <  1:1561914/100‑1 (MQ=255)
                                          aGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGt                                                          <  1:508996/100‑1 (MQ=255)
                                           ggAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTc                                                         <  1:2833729/100‑1 (MQ=255)
                                           ggAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTc                                                         <  1:1994620/100‑1 (MQ=255)
                                            gAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCt                                                        <  1:894788/100‑1 (MQ=255)
                                             aGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCTACATGGCAAAACTACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTa                                                       >  2:343253/1‑100 (MQ=255)
                                              gTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTAt                                                      >  2:1741441/1‑100 (MQ=255)
                                              gTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTAt                                                      >  2:2936501/1‑100 (MQ=255)
                                               ttCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTAtt                                                     >  1:4051603/1‑100 (MQ=255)
                                               ttCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTAtt                                                     <  1:3792101/100‑1 (MQ=255)
                                                tCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTa                                                    <  2:914164/100‑1 (MQ=255)
                                                tCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTa                                                    <  1:3904373/100‑1 (MQ=255)
                                                tCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTa                                                    <  1:1344432/100‑1 (MQ=255)
                                                tCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTa                                                    >  1:323540/1‑100 (MQ=255)
                                                 cGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTAc                                                   >  2:3465119/1‑100 (MQ=255)
                                                 cGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTAc                                                   >  1:931682/1‑100 (MQ=255)
                                                 cGAGAAAAATTATTGTGCAGATATTATACAGGATGTGAAGTATCGACATGTCAAAAATAAAGGATTTATCTCTTTTGACGGTAAAACCTTCGTATATTa                                                    >  2:3490786/1‑99 (MQ=255)
                                                   agaAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACct                                                 >  2:1395866/1‑100 (MQ=255)
                                                    gaAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACctc                                                <  2:1159093/100‑1 (MQ=255)
                                                    gaAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACctc                                                <  1:1344739/100‑1 (MQ=255)
                                                      aaaaTTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTa                                              >  1:1657434/1‑100 (MQ=255)
                                                       aaaTTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTAt                                             >  2:856884/1‑100 (MQ=255)
                                                         attattGTGCAGATATTATGCAGTTTGTGAAGTAAAAAAAAAAAAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATcc                                           <  2:675781/100‑1 (MQ=255)
                                                         attattGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATcc                                           <  1:3330388/100‑1 (MQ=255)
                                                         attattGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATcc                                           <  1:1573494/100‑1 (MQ=255)
                                                         attattGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATcc                                           >  2:1536820/1‑100 (MQ=255)
                                                          ttattGTGCAGATATTAGGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGCGTTTATCTCGTTTGACGGTAAAACCTTCGTCTATTACCTCTATc                                            <  2:1492716/98‑1 (MQ=255)
                                                            attGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGt                                        <  2:2702652/100‑1 (MQ=255)
                                                                tgCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACt                                    <  1:3361919/100‑1 (MQ=255)
                                                                  cAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTca                                  >  2:2617006/1‑100 (MQ=255)
                                                                  cAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTca                                  >  1:4044467/1‑100 (MQ=255)
                                                                   aGATATTATGCAGGTTGTGAAGAAAAAAAAAAAAAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTca                                  <  2:360291/99‑1 (MQ=255)
                                                                    gATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTcaca                                >  1:2400980/1‑100 (MQ=255)
                                                                    gATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTcaca                                >  2:1626801/1‑100 (MQ=255)
                                                                     atatTATGCAGGTTGTGAATTAACAAAAAAAAAAAAATACAGGTTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTcaca                                <  2:589264/99‑1 (MQ=255)
                                                                        ttatGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATc                             >  1:1448725/1‑99 (MQ=255)
                                                                        ttatGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCg                            >  2:1785714/1‑100 (MQ=255)
                                                                          ntGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCAAAATCGTa                          >  2:1946036/2‑100 (MQ=255)
                                                                          atGCAGGTTGTGAAGTATCGTCATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTa                          >  2:2516169/1‑100 (MQ=255)
                                                                          atGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTa                          <  2:3019969/100‑1 (MQ=255)
                                                                           tGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAg                         <  1:3460526/100‑1 (MQ=255)
                                                                            gCATGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACACTCGTAg                         <  2:3029517/99‑1 (MQ=255)
                                                                             cAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTc                       <  1:2170217/100‑1 (MQ=255)
                                                                                  tgtgAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGata                  >  1:1484848/1‑100 (MQ=255)
                                                                                   gtgAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGatat                 >  1:2841436/1‑100 (MQ=255)
                                                                                     gAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATAttt               >  2:1105239/1‑100 (MQ=255)
                                                                                        gTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTg            >  2:3227354/1‑100 (MQ=255)
                                                                                        gTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTg            >  2:91846/1‑100 (MQ=255)
                                                                                        gTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTg            >  2:2423644/1‑100 (MQ=255)
                                                                                          aTCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATAtttttttt          >  2:231209/1‑97 (MQ=255)
                                                                                          aTCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCt          >  2:3931007/1‑100 (MQ=255)
                                                                                            cGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCtt         >  2:1660468/1‑99 (MQ=255)
                                                                                            cGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTg        >  2:2223264/1‑100 (MQ=255)
                                                                                             gACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTgg       >  2:3349069/1‑100 (MQ=255)
                                                                                               cATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTc     <  1:357143/100‑1 (MQ=255)
                                                                                                  gTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTgg  <  1:3297287/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TCGGTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGG  >  NC_000913/875615‑875812

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: