Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,295,025:1 +T 9.5% intergenic (+41/+29) 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,0251.T 9.5% 318.8 / 28.9 84intergenic (+41/+29)slyB/nhaA‑Iouter membrane lipoprotein/Na(+)/H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base . (37/39);  new base T (3/5);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGCTACCCGG  >  NC_002947/1294878‑1295173
                                                                                                                                                    |                                                                                                                                                     
gcagcaTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCatg                                                                                                                                                      <  1:286157/149‑3 (MQ=255)
    caTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCgg                                                                                                                                                  <  6:120968/149‑1 (MQ=255)
        cgcgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGt                                                                                                                                               <  4:403810/148‑1 (MQ=255)
        cgcgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCGGGGt                                                                                                                                               <  1:133874/148‑1 (MQ=255)
         gcgcCTATGTGCAGGCCGTGCAGGAGAACGAATTCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGGTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGttt                                                                                                                                             <  6:418563/149‑1 (MQ=255)
         gcgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGttt                                                                                                                                             >  6:353515/1‑149 (MQ=255)
             cTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTc                                                                                                                                         <  8:50858/149‑1 (MQ=255)
             cTATGTGCAGGCCGTGCAGGAGAACGAAATCTGCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGATTCTCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTc                                                                                                                                         <  8:107631/149‑1 (MQ=255)
              tATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCa                                                                                                                                        <  4:360711/149‑1 (MQ=255)
               aTGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCAt                                                                                                                                       >  1:129170/1‑149 (MQ=255)
                  tgCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGcc                                                                                                                                    >  4:73893/1‑149 (MQ=255)
                    cAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTc                                                                                                                                  <  6:372518/149‑1 (MQ=255)
                     aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCt                                                                                                                                   >  6:311138/1‑147 (MQ=255)
                     aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCa                                                                                                                                 <  3:350941/149‑1 (MQ=255)
                      ggCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAg                                                                                                                                <  2:137947/149‑1 (MQ=255)
                      ggCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGGGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAg                                                                                                                                <  5:49730/149‑1 (MQ=255)
                        ccGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGcc                                                                                                                              <  6:132981/149‑1 (MQ=255)
                        ccGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGc                                                                                                                               >  3:427723/1‑148 (MQ=255)
                           tGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAg                                                                                                                           >  4:90776/1‑149 (MQ=255)
                             cAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTc                                                                                                                         <  2:129170/149‑1 (MQ=255)
                                gagaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAgg                                                                                                                      >  5:412799/1‑149 (MQ=255)
                                gagaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAgg                                                                                                                      <  6:334445/149‑1 (MQ=255)
                                gagaACGAAATCTTCCGGGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGGACCAGTCGCGGTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAgg                                                                                                                      <  8:85719/149‑1 (MQ=255)
                                     cGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCGTGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATg                                                                                                                 >  4:132344/1‑149 (MQ=255)
                                     cGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCGTGCCTCAGCCAAGTCAGGCAAt                                                                                                                  >  2:98006/1‑148 (MQ=255)
                                     cGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCggg                                                                                                                                                 >  3:413948/1‑117 (MQ=38)
                                     cGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCggg                                                                                                                                                 <  4:413948/117‑1 (MQ=38)
                                      gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAg                                                                                                                                <  8:60120/133‑1 (MQ=255)
                                      gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAg                                                                                                                                >  7:60120/1‑133 (MQ=255)
                                      gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCGGGGtttttta                                                                                                                                         >  8:390237/1‑123 (MQ=38)
                                      gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCGGGGtttttta                                                                                                                                         <  7:390237/124‑2 (MQ=38)
                                                gtgtTGGCGACCGAGTGCGCATCGTGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCg                                                                                                      <  3:132344/149‑1 (MQ=255)
                                                  gtTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCAt                                                                                                                                       <  6:273508/114‑1 (MQ=255)
                                                  gtTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCAt                                                                                                                                       <  2:382571/114‑1 (MQ=255)
                                                  gtTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCAt                                                                                                                                       >  5:273508/1‑114 (MQ=255)
                                                  gtTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCAt                                                                                                                                       >  1:382571/1‑114 (MQ=255)
                                                     ggCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTgg                                                                                                 >  5:261994/1‑149 (MQ=255)
                                                     ggCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTg                                                                                                  <  3:73893/148‑1 (MQ=255)
                                                        gACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTTGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCa                                                                                              >  8:315037/1‑149 (MQ=255)
                                                          ccGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGc                                                                                                       <  2:395582/138‑1 (MQ=255)
                                                          ccGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGc                                                                                                       >  1:395582/1‑138 (MQ=255)
                                                          ccGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCAtt                                                                                            <  1:148612/149‑1 (MQ=255)
                                                              gTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCgg                                                                                        <  1:32464/149‑1 (MQ=255)
                                                               tGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGt                                                                                       <  3:90776/149‑1 (MQ=255)
                                                               tGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCT‑CCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGt                                                                                       >  7:128643/1‑148 (MQ=255)
                                                                 cgcATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAAGGCTGGTTTTGCGGCTGGCCATTGCGGTGa                                                                                     >  8:405829/1‑149 (MQ=255)
                                                                  gcATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑c                                                                                   >  5:238582/1‑149 (MQ=255)
                                                                     tcatGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGc                                                                                <  5:353515/149‑1 (MQ=255)
                                                                        tGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCg                                                                               >  7:1723/1‑147 (MQ=255)
                                                                            tGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGcc                                                                          >  7:9200/1‑148 (MQ=255)
                                                                                gACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCa                                                                     <  8:1723/149‑1 (MQ=255)
                                                                                 aCGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCAc                                                                    <  4:427723/149‑1 (MQ=255)
                                                                                   ggTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCg                                                                  <  5:311138/149‑1 (MQ=255)
                                                                                   ggTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCg                                                                  <  6:238582/149‑1 (MQ=255)
                                                                                       ccAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGATCGGCTACGCCG‑TCACCgcaat                                                              <  1:127137/149‑3 (MQ=255)
                                                                                           tCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCgcagca                                                             >  6:211415/1‑146 (MQ=255)
                                                                                             gcgTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCT‑CCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAgaa                                                        <  8:128643/148‑1 (MQ=255)
                                                                                                  tCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                                   >  3:178959/1‑149 (MQ=255)
                                                                                                   cGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                                   <  6:412799/148‑1 (MQ=255)
                                                                                                    gCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                                   >  1:355603/1‑147 (MQ=255)
                                                                                                      aCTGACCGTTGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGt                                               <  7:315037/149‑1 (MQ=255)
                                                                                                         gACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGt                                               >  1:384322/1‑146 (MQ=255)
                                                                                                                      taataaAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCgg                               >  2:267338/1‑149 (MQ=255)
                                                                                                                       atgaaaaaaCCCCGAC‑TGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCa                                                 >  3:340332/4‑130 (MQ=37)
                                                                                                                       atgaaaaaaCCCCGAC‑TGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCa                                                 <  4:340332/127‑1 (MQ=37)
                                                                                                                        ataaAAATCCCCGACATGGCTGCTGCCA‑GTCGGGTTTTTTTCATTCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGTCTGGCCATTGCGTTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAATAATGACCATGCGGTc                             <  4:422294/149‑1 (MQ=255)
                                                                                                                             aaaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGc                        >  6:236280/1‑149 (MQ=255)
                                                                                                                               aCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCg                      <  2:355603/149‑1 (MQ=255)
                                                                                                                                 cccGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGc                    >  3:85510/1‑149 (MQ=255)
                                                                                                                                 cccGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGc                    <  5:211415/149‑1 (MQ=255)
                                                                                                                                    gACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTa                 <  6:261994/149‑1 (MQ=255)
                                                                                                                                     aCATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTa                 >  4:266750/1‑148 (MQ=255)
                                                                                                                                      cATGGCTGCTGCCA‑GTCGGGGTTTTTTCGTGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACt               <  1:98006/149‑1 (MQ=255)
                                                                                                                                      cATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                                   <  7:306038/113‑1 (MQ=255)
                                                                                                                                      cATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                                   >  8:306038/1‑113 (MQ=255)
                                                                                                                                       aTGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTACCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTc              <  3:266750/149‑1 (MQ=255)
                                                                                                                                        tGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTct             >  4:134284/1‑148 (MQ=255)
                                                                                                                                         ggCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCAc            >  5:286313/1‑149 (MQ=255)
                                                                                                                                            tgctgcCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGc         >  7:400454/1‑149 (MQ=255)
                                                                                                                                            tgctgcCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGc         >  3:164575/1‑149 (MQ=255)
                                                                                                                                                gcCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGCTAcc     >  5:228857/1‑149 (MQ=255)
                                                                                                                                                gcCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGCTAcc     >  8:393489/1‑149 (MQ=255)
                                                                                                                                                  cA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCg                   <  1:100387/133‑1 (MQ=255)
                                                                                                                                                  cA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCg                   >  2:100387/1‑133 (MQ=255)
                                                                                                                                                  cA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGCTACCCg   >  2:415982/1‑149 (MQ=255)
                                                                                                                                                   a‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGCTACCCgg  <  2:384322/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGA‑CGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGCTACCCGG  >  NC_002947/1294878‑1295173

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