Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,295,027 T→A 6.7% intergenic (+37/+33) slyB → / ← nhaA‑I outer membrane lipoprotein/Na(+)/H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,295,0270TA 6.7% 222.8 / 4.5 72intergenic (+37/+33)slyB/nhaA‑Iouter membrane lipoprotein/Na(+)/H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base A (3/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGCT  >  NC_002947/1294882‑1295173
                                                                                                                                                 |                                                                                                                                                   
tgGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCa                                                                                                                                                 >  1:8049/1‑150 (MQ=255)
   cagcaTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTc                                                                                                                                             >  2:196684/1‑150 (MQ=255)
    agcaTGCGCGCCTATGTGCAGGCCGTATAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCg                                                                                                                                            <  1:479373/150‑1 (MQ=255)
     gcaTGCGAGCCTGTGTGCAGGCCGCCCGCGGTAATGATAGCTGCCGTGGTGGCTACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCATGTCgg                                                                                                                                           <  2:1646231/150‑1 (MQ=255)
      caTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCggg                                                                                                                                          <  1:1883002/150‑1 (MQ=255)
      caTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACTGGCAGCAGCCATGTCggg                                                                                                                                          <  1:1118570/150‑1 (MQ=255)
       aTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCgggg                                                                                                                                         <  1:1823272/150‑1 (MQ=255)
        tGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGt                                                                                                                                        <  1:1505217/150‑1 (MQ=255)
             gcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGtttttt                                                                                                                                   <  2:280502/150‑1 (MQ=255)
               cTATGTGCAGGCTGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCa                                                                                                                                 <  2:985638/150‑1 (MQ=255)
               cTATGTGCAGGCCGTGTAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCa                                                                                                                                 >  2:1279882/1‑150 (MQ=255)
               cTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCa                                                                                                                                 >  1:1428272/1‑150 (MQ=255)
               cTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCa                                                                                                                                 >  2:1761065/1‑150 (MQ=255)
                    tgCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTc                                                                                                                                  <  2:167462/144‑1 (MQ=255)
                    tgCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTc                                                                                                                                  >  1:167462/1‑144 (MQ=255)
                      cAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCa                                                                                                                          <  2:473356/150‑1 (MQ=255)
                       aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAg                                                                                                                         <  2:1428272/150‑1 (MQ=255)
                       aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAg                                                                                                                         <  1:1026386/150‑1 (MQ=255)
                            gTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATGATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAg                                                                                                                    <  2:208683/150‑1 (MQ=255)
                                aGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATAATGAATGTCACCGGGACCAGTCGCGTTTCGCACTGTCCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAg                                                                                                                <  1:1376332/150‑1 (MQ=255)
                                 ggAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCggg                                                                                                                                          >  2:598110/1‑123 (MQ=255)
                                 ggAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCggg                                                                                                                                          <  1:598110/123‑1 (MQ=255)
                                  gagaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGc                                                                                                              >  2:664889/1‑150 (MQ=255)
                                      aCGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATg                                                                                                          <  1:1200812/150‑1 (MQ=255)
                                      aCGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATg                                                                                                          <  1:664889/150‑1 (MQ=255)
                                           aTCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGt                                                                                                     <  2:102095/150‑1 (MQ=255)
                                              ttCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGtttt                                                                                                  >  1:282795/1‑150 (MQ=255)
                                              ttCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGtttt                                                                                                  >  1:733307/1‑150 (MQ=255)
                                                       ggCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGc                                                                                         >  1:624162/1‑150 (MQ=255)
                                                       ggCGACCGAGTACGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGc                                                                                         >  2:551259/1‑150 (MQ=255)
                                                          gACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCAt                                                                                      >  1:1253573/1‑150 (MQ=255)
                                                              gAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGCGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCg                                                                                  <  1:1023789/150‑1 (MQ=255)
                                                                gTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGt                                                                                <  1:765825/150‑1 (MQ=255)
                                                                   cgcATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGAc                                                                             <  1:282154/150‑1 (MQ=255)
                                                                     catcatGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACgg                                                                           >  2:638614/1‑150 (MQ=255)
                                                                     catcatGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGGCGGGGTTTTTTCAGGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACgg                                                                           >  2:1035335/1‑150 (MQ=255)
                                                                      atcatGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGc                                                                          <  1:1140161/150‑1 (MQ=255)
                                                                        catGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGt                                                                        <  2:1238678/150‑1 (MQ=255)
                                                                          tGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGt                                                                                                                                        >  1:1101595/1‑84 (MQ=255)
                                                                          tGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGt                                                                                                                                        <  2:1101595/84‑1 (MQ=255)
                                                                          tGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAg                                                                      >  1:1033966/1‑150 (MQ=255)
                                                                             cTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGTGGTGTGCGGTTACGGTGACGGCGGAGCCg                                                                   >  2:904195/1‑150 (MQ=255)
                                                                              tGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGa                                                                  >  1:527007/1‑150 (MQ=255)
                                                                              tGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGa                                                                  <  2:8049/150‑1 (MQ=255)
                                                                                tAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGa                                                                  <  1:1280894/148‑1 (MQ=255)
                                                                                tAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTATTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGa                                                                  >  2:1280894/1‑148 (MQ=255)
                                                                                            gTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAg                                                    >  2:1231355/1‑150 (MQ=255)
                                                                                            gTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAg                                                    >  2:1503883/1‑150 (MQ=255)
                                                                                              cgcgTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaa                                                  <  1:904195/150‑1 (MQ=255)
                                                                                                  tttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTg                                                                                           >  1:1111098/1‑105 (MQ=255)
                                                                                                  tttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTg                                                                                           <  2:1111098/105‑1 (MQ=255)
                                                                                                  tttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaangaa                                              <  2:1202404/150‑1 (MQ=255)
                                                                                                     cGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCa                                           <  2:1594623/150‑1 (MQ=255)
                                                                                                         cTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGCACGAACCAGTaa                                       >  1:374883/1‑150 (MQ=255)
                                                                                                                    gCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATg                            <  1:551259/150‑1 (MQ=255)
                                                                                                                       gcatgaaaaaaCCCCGAGTGGCAGCAGCCATGGCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGCTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCgg                         >  2:118635/6‑150 (MQ=255)
                                                                                                                        taataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGt                        <  1:196684/150‑1 (MQ=255)
                                                                                                                        taataaAAAACCCCGACATGGCTGCTGCCAGTAGGGGTCTTTTCCTCCATCACCCAAGGCAGACAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGt                        >  1:515342/1‑150 (MQ=255)
                                                                                                                         aataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCAt                             <  2:1315616/144‑1 (MQ=255)
                                                                                                                         aataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCAt                             >  1:1315616/1‑144 (MQ=255)
                                                                                                                         aataaAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTc                       >  1:48153/1‑150 (MQ=255)
                                                                                                                                 aCCCCGACTGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGTCGGCGTAGCCGATCACCGCAGCAAAGAACGAACCGGTAAGAATGCCCATGCGGTCCCTGCCgg               >  2:1613271/1‑150 (MQ=255)
                                                                                                                                  ccccGACTGGCAGCAGCCATGTCGGGGTTTTTTCATTCTTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGc              >  2:1520529/1‑150 (MQ=255)
                                                                                                                                   cccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACg                                                                            >  2:1831859/1‑87 (MQ=255)
                                                                                                                                   cccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACg                                                                            <  1:1831859/87‑1 (MQ=255)
                                                                                                                                   cccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCg             >  1:1258979/1‑150 (MQ=255)
                                                                                                                                   cccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCg             >  1:872152/1‑150 (MQ=255)
                                                                                                                                   cccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCg             >  2:1409489/1‑150 (MQ=255)
                                                                                                                                   cccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGGAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCg             >  2:1409146/1‑150 (MQ=255)
                                                                                                                                    ccGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGGAAGAATGCCCATGCGGTCCATGCCGGCgg            >  1:1713726/1‑149 (MQ=255)
                                                                                                                                     cGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTa           >  1:443832/1‑150 (MQ=255)
                                                                                                                                      gACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTAc          <  2:48153/150‑1 (MQ=255)
                                                                                                                                            gctgctGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGAt                                                                 >  2:52337/1‑89 (MQ=255)
                                                                                                                                            gctgctGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGAt                                                                 <  1:52337/89‑1 (MQ=255)
                                                                                                                                              tgctgcCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGGCGGCGGAGGCGCTTGCCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGcc  >  2:689252/1‑149 (MQ=255)
                                                                                                                                              tgctgcCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGcc  >  2:1393768/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGCT  >  NC_002947/1294882‑1295173

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