Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,245,027 T→G 5.6% D13A (GAC→GCC)  aruH ← arginine‑‑pyruvate transaminase AruH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,245,0270TG 5.6% 367.7 / 11.2 125D13A (GAC→GCC) aruHarginine‑‑pyruvate transaminase AruH
Reads supporting (aligned to +/- strand):  ref base T (53/64);  new base G (4/3);  total (58/67)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 125 total aligned reads displayed.

ATGGCCAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCG  >  NC_002947/4244879‑4245173
                                                                                                                                                     |                                                                                                                                                   
aTGGCCAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCAACGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTc                                                                                                                                                    <  2:2118326/150‑1 (MQ=255)
     cAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCgg                                                                                                                                               <  2:2058831/150‑1 (MQ=255)
     cAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCgg                                                                                                                                               >  1:257424/1‑150 (MQ=255)
        cgcgcAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGa                                                                                                                                           <  2:2846585/150‑1 (MQ=255)
               gCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTg                                                                                                                                    <  1:1815106/150‑1 (MQ=255)
                cTGTCGATGGCCGCCGCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTgg                                                                                                                                   <  2:2148284/150‑1 (MQ=255)
                 tGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTggg                                                                                                                                  >  2:532485/1‑150 (MQ=255)
                  gTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGt                                                                                                                                 <  2:847096/150‑1 (MQ=255)
                    cGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGa                                                                                                                               <  1:570796/150‑1 (MQ=255)
                     gATGGCCGCCTCGACGATAGGCGCGCGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGc                                                                                                                                        >  2:2803135/1‑140 (MQ=255)
                     gATGGCCGCCTCGACGATAGGCGCGCGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGc                                                                                                                                        <  1:2803135/140‑1 (MQ=255)
                      gttgccgccTGGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGc                                                                                                                             <  2:473354/147‑1 (MQ=255)
                      aTGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGc                                                                                                                             <  1:879045/150‑1 (MQ=255)
                         gccgccTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTg                                                                                                                          <  1:428309/150‑1 (MQ=255)
                           cgccTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGGGGTGGTGTGTCTTCCAGGCAGCAGCTGCGTCTTTGGCCATGCTTTGGGTGAGCTTGGc                                                                                                                        >  2:939375/1‑150 (MQ=255)
                           cgccTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACACACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGc                                                                                                                        <  2:2647147/150‑1 (MQ=255)
                             ccTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGACGCATCGCCGGCGACGGCGCTGCggcctg                                                                                                                                  >  1:2514066/1‑133 (MQ=255)
                                     aTAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATg                                                                                                              <  2:2380323/150‑1 (MQ=255)
                                      tAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCACATgg                                                                                                             >  1:2597064/1‑150 (MQ=255)
                                      tAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGACCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATgg                                                                                                             <  1:78150/150‑1 (MQ=255)
                                       aGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGc                                                                                                            <  1:745482/150‑1 (MQ=255)
                                         gcgcGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCaa                                                                                                          >  2:2275787/1‑150 (MQ=255)
                                          cgcgGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAAt                                                                                                         <  2:2620062/150‑1 (MQ=255)
                                           gcgGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAAtc                                                                                                        <  1:2197513/150‑1 (MQ=255)
                                                aGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCtt                                                                                                   <  2:392703/150‑1 (MQ=255)
                                                    tCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCaa                                                                                               >  2:1058246/1‑150 (MQ=255)
                                                     cGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAAt                                                                                              <  1:371744/150‑1 (MQ=255)
                                                      gAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATg                                                                                             <  2:1686997/150‑1 (MQ=255)
                                                         gTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCa                                                                                          >  2:718392/1‑150 (MQ=255)
                                                                   cccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGACGCATCGCCGGCGACGGCGCTGCggc                                                                                                                                     >  2:2136167/1‑95 (MQ=37)
                                                                   cccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGACGCATCGCCGGCGACGGCGCTGCggc                                                                                                                                     <  1:2136167/97‑3 (MQ=37)
                                                                    ccaccGACAGCAGCAGGATGTCCTTGTCGTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcga                                                                               >  1:2716941/1‑150 (MQ=255)
                                                                    ccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcga                                                                               >  1:2054778/1‑150 (MQ=255)
                                                                    ccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcga                                                                               >  2:601204/1‑150 (MQ=255)
                                                                    ccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTc                                                                                 >  2:1341860/1‑148 (MQ=255)
                                                                    ccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTc                                                                                 <  1:1341860/148‑1 (MQ=255)
                                                                    ccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATATAAGTcga                                                                               >  2:2520948/1‑150 (MQ=255)
                                                                    ccaccGACAGCAGCAGGATGTCCTTGTCCTAAGCCTGCAGCACCAGGGCGCGGTTGTGCCTGTCCCAGGCGGCCGCCGGACACTGCGGCAGTTGCCCT‑‑GTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTGCAGTcga                                                                               >  1:1465743/1‑150 (MQ=255)
                                                                     caccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcgac                                                                              >  2:181369/1‑150 (MQ=255)
                                                                       ccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcgacga                                                                            <  2:1757891/150‑1 (MQ=255)
                                                                               agcagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTc                                                                    <  2:257424/150‑1 (MQ=255)
                                                                               agcagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTc                                                                    >  1:1115651/1‑150 (MQ=255)
                                                                                    gATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTcc                                                                                                     <  1:1036825/112‑1 (MQ=255)
                                                                                    gATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTcc                                                                                                     >  2:1036825/1‑112 (MQ=255)
                                                                                    gATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGtt                                                               >  1:1031069/1‑150 (MQ=255)
                                                                                     aTGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTg                                                              <  2:1085676/150‑1 (MQ=255)
                                                                                        tccttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGaa                                                           <  2:1544630/150‑1 (MQ=255)
                                                                                        tccttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGaa                                                           <  1:2275787/150‑1 (MQ=255)
                                                                                          cttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCtt                                                                                                   <  2:2188631/108‑1 (MQ=255)
                                                                                          cttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAAtt                                                         >  1:704140/1‑150 (MQ=255)
                                                                                           ttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGt                                                                          >  1:1711652/1‑132 (MQ=255)
                                                                                           ttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGt                                                                          <  2:1711652/132‑1 (MQ=255)
                                                                                            tgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGc                                                       <  1:181369/150‑1 (MQ=255)
                                                                                             gtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGcc                                                      <  2:1115651/150‑1 (MQ=255)
                                                                                              tcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCACCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGccc                                                     <  2:2054778/150‑1 (MQ=255)
                                                                                                ctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCa                                                   >  1:287921/1‑150 (MQ=255)
                                                                                                 tCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGc                                                       <  2:2354314/145‑1 (MQ=255)
                                                                                                  cAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATa                                                 <  1:601204/150‑1 (MQ=255)
                                                                                                   aGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAAtt                                                         >  2:2691945/1‑141 (MQ=255)
                                                                                                   aGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAAtt                                                         <  1:2691566/141‑1 (MQ=255)
                                                                                                   aGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAAtt                                                         >  2:2691566/1‑141 (MQ=255)
                                                                                                      cTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGgt                                             >  1:945598/1‑150 (MQ=255)
                                                                                                      cTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGgt                                             >  1:3045766/1‑150 (MQ=255)
                                                                                                       tGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGgtg                                            <  2:945598/150‑1 (MQ=255)
                                                                                                        gCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGgtgt                                           >  2:1396080/1‑150 (MQ=255)
                                                                                                         cAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGCCGGTCTTTTTGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTa                                          <  1:532485/150‑1 (MQ=255)
                                                                                                         cAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTa                                          >  1:2929108/1‑150 (MQ=255)
                                                                                                         cAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTa                                          <  1:1058246/150‑1 (MQ=255)
                                                                                                          aGCGCCAGGGCGCGGTAGTGGATGTCACAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTAc                                         >  1:618723/1‑150 (MQ=255)
                                                                                                           gcgcCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACa                                        >  1:1859886/1‑150 (MQ=255)
                                                                                                              ccAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCggg                                                                 >  2:2521215/1‑122 (MQ=255)
                                                                                                               ctgggCGCCGTAGTATATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATc                                    <  2:2184185/148‑1 (MQ=255)
                                                                                                               cAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGgtg                                            <  1:61877/142‑1 (MQ=255)
                                                                                                               cAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGgtg                                            >  2:61877/1‑142 (MQ=255)
                                                                                                                 gggCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACg                                               <  2:1068876/137‑1 (MQ=255)
                                                                                                                  ggCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCCATTTCAGGCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGc                                 <  1:2520948/150‑1 (MQ=255)
                                                                                                                  ggCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCgag                                 >  1:2184185/1‑148 (MQ=255)
                                                                                                                  ggCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGc                                 <  2:1859886/150‑1 (MQ=255)
                                                                                                                   gcgcGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCg                                >  2:599555/1‑150 (MQ=255)
                                                                                                                   gcgcGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCg                                >  2:3048175/1‑150 (MQ=255)
                                                                                                                    cgcgGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGa                               <  1:1396080/150‑1 (MQ=255)
                                                                                                                    cgcgGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGa                               <  1:1029597/150‑1 (MQ=255)
                                                                                                                        gTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACg                                               <  2:231047/130‑1 (MQ=255)
                                                                                                                        gTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACg                                               >  1:231047/1‑130 (MQ=255)
                                                                                                                              gATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGa                               >  2:2357867/1‑140 (MQ=255)
                                                                                                                              gATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGa                               <  1:2357867/140‑1 (MQ=255)
                                                                                                                                  tCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTg                 >  1:2907154/1‑150 (MQ=255)
                                                                                                                                    ccAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGAt               >  1:2781239/1‑150 (MQ=255)
                                                                                                                                    ccAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGAt               <  2:287921/150‑1 (MQ=255)
                                                                                                                                     cAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAg                                                                                   >  2:409087/1‑81 (MQ=255)
                                                                                                                                     cAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGAtt              <  2:1465743/150‑1 (MQ=255)
                                                                                                                                      aGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGAttt             >  1:1905532/1‑150 (MQ=255)
                                                                                                                                           gCATCGCCGGCGACGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTcc        >  2:2487452/1‑150 (MQ=255)
                                                                                                                                              gcgcCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCAtc     <  2:2781239/150‑1 (MQ=255)
                                                                                                                                              atcgcCGGCGACGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGgt                                             <  2:2571724/109‑1 (MQ=255)
                                                                                                                                              atcgcCGGCGACGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGgt                                             >  1:2571724/3‑111 (MQ=255)
                                                                                                                                                 gcCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCAtctc   >  2:2085374/1‑150 (MQ=255)
                                                                                                                                                 gcCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGCGCTGGAATCGGCGAGCAACGCGCCGTTGATTTCTTCCAtctc   >  1:2337003/1‑150 (MQ=255)
                                                                                                                                                  ccGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCg  >  1:2449844/1‑150 (MQ=255)
                                                                                                                                                   cGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGccc                                                     >  1:656925/1‑98 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATGGCCAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCG  >  NC_002947/4244879‑4245173

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