Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,249,979 C→A 6.4% G12V (GGC→GTC)  aruH ← arginine‑‑pyruvate transaminase AruH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,249,9790CA 6.4% 1162.3 / 47.7 346G12V (GGC→GTC) aruHarginine‑‑pyruvate transaminase AruH
Reads supporting (aligned to +/- strand):  ref base C (168/156);  new base A (14/8);  total (182/164)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 351 total aligned reads displayed.

GCCAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCGGCCTGG  >  NC_002947/4249831‑4250128
                                                                                                                                                    |                                                                                                                                                     
gCCAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGcc                                                                                                                                                      >  1:6716396/1‑150 (MQ=255)
 ccAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCg                                                                                                                                                     >  2:380799/1‑150 (MQ=255)
  cAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCgg                                                                                                                                                    >  2:4567108/1‑150 (MQ=255)
   aGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGc                                                                                                                                                   >  1:4941195/1‑150 (MQ=255)
           ggCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGtt                                                                                                                                           <  2:1380537/150‑1 (MQ=255)
               gTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGt                                                                                                                                       <  1:2348567/150‑1 (MQ=255)
                tCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTg                                                                                                                                      >  2:6754298/1‑150 (MQ=255)
                tCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTg                                                                                                                                      >  1:6287328/1‑150 (MQ=255)
                  gATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAg                                                                                                                                    >  2:2665186/1‑150 (MQ=255)
                      gccgccTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTg                                                                                                                                <  2:6160336/150‑1 (MQ=255)
                       ccgccTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTgg                                                                                                                               >  1:3273899/1‑150 (MQ=255)
                          ccTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGt                                                                                                                            >  1:2807997/1‑150 (MQ=255)
                           cTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTa                                                                                                                           <  2:838113/150‑1 (MQ=255)
                            tCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGTGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAg                                                                                                                          >  2:6345609/1‑150 (MQ=255)
                            tCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCCGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAg                                                                                                                          >  2:1211094/1‑150 (MQ=255)
                            tCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAg                                                                                                                          >  2:1440077/1‑150 (MQ=255)
                            tCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAg                                                                                                                          <  1:6792484/150‑1 (MQ=255)
                                   tAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATgg                                                                                                                   <  1:4771659/150‑1 (MQ=255)
                                   tAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATgg                                                                                                                   <  2:5790895/150‑1 (MQ=255)
                                      gcgcGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCaa                                                                                                                >  2:5135155/1‑150 (MQ=255)
                                          gggAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAAtctc                                                                                                            >  1:5847044/1‑150 (MQ=255)
                                          gggAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGACGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAAtctc                                                                                                            <  1:5361957/150‑1 (MQ=255)
                                            gAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCt                                                                                                          >  2:4404961/1‑150 (MQ=255)
                                             aGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCtt                                                                                                         <  2:6266672/150‑1 (MQ=255)
                                               tgtCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTc                                                                                                       <  1:2283574/150‑1 (MQ=255)
                                                gtCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCa                                                                                                      <  2:3188639/150‑1 (MQ=255)
                                                  cGAAGTCCGGGTCACCCACCGACATCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCTATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAAt                                                                                                    <  1:2299072/150‑1 (MQ=255)
                                                   gAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATg                                                                                                   >  1:6850752/1‑150 (MQ=255)
                                                   gAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTGTCAATg                                                                                                   >  2:2827222/1‑150 (MQ=255)
                                                      gTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCCGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCa                                                                                                >  2:1965861/1‑150 (MQ=255)
                                                       tCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCaa                                                                                               >  2:2314542/1‑150 (MQ=255)
                                                       tCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCaa                                                                                               >  2:1792962/1‑150 (MQ=255)
                                                         cGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAAtt                                                                                             >  2:645815/1‑150 (MQ=255)
                                                          gggTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAAttt                                                                                            <  1:3986846/150‑1 (MQ=255)
                                                          gggTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAAttt                                                                                            <  1:4825394/150‑1 (MQ=255)
                                                           ggTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTc                                                                                           >  2:5252553/1‑150 (MQ=255)
                                                            gaccccccccGACTGCGGGCGGCTGTCTTTGCCCGCCGCCGGCTGCGCCGGGGCGGGTTTGGGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCa                                                                                          <  1:370329/146‑1 (MQ=255)
                                                            gTCACCCACCGACAGCAGCAGGATGTCCGTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCa                                                                                          <  1:5428574/150‑1 (MQ=255)
                                                               acccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTc                                                                                       >  1:1257189/1‑150 (MQ=255)
                                                                cccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcg                                                                                      >  1:1734482/1‑150 (MQ=255)
                                                                 ccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcga                                                                                     >  1:3305558/1‑150 (MQ=255)
                                                                 ccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcga                                                                                     >  1:825843/1‑150 (MQ=255)
                                                                  caccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcgac                                                                                    >  1:5573267/1‑150 (MQ=255)
                                                                    ccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcgacga                                                                                  >  2:4260606/1‑150 (MQ=255)
                                                                    ccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGAATGGCAATCTCCTTTCAATGGCAATTTCAGTcgacga                                                                                  <  2:5635190/150‑1 (MQ=255)
                                                                         agcagcagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCCATCTCCTTTCAATGGCACTTTCAGTCGACGAGTATc                                                                             >  2:4392411/1‑150 (MQ=255)
                                                                          gcagcagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCt                                                                            >  1:2091213/1‑150 (MQ=255)
                                                                            agcagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTc                                                                          <  1:6070888/150‑1 (MQ=255)
                                                                             gcagGATGTCCTTGTCCTCAGACTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTaacgcgc                                                                         >  2:5314143/1‑144 (MQ=255)
                                                                              cagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCgg                                                                        >  2:1559682/1‑150 (MQ=255)
                                                                              cagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCgg                                                                        <  1:4470222/150‑1 (MQ=255)
                                                                                ggATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGt                                                                      <  2:6838555/150‑1 (MQ=255)
                                                                                 gATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGtt                                                                     <  2:4831610/150‑1 (MQ=255)
                                                                                        ttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTg                                                              >  1:5802728/1‑150 (MQ=255)
                                                                                          gtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGcc                                                            <  1:380799/150‑1 (MQ=255)
                                                                                           tcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGccc                                                           <  1:5045561/150‑1 (MQ=255)
                                                                                             ctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCa                                                         >  1:112121/1‑150 (MQ=255)
                                                                                             ctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCa                                                         >  1:3525344/1‑150 (MQ=255)
                                                                                             ctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCa                                                         >  2:1114005/1‑150 (MQ=255)
                                                                                             ctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTCTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCa                                                         >  1:4422260/1‑150 (MQ=255)
                                                                                              tCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCGGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCACTTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCAt                                                        >  2:4154067/1‑150 (MQ=255)
                                                                                              tCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCTCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGTTTGGAATTGCCCCAt                                                        >  1:1568612/1‑150 (MQ=255)
                                                                                              tCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCAt                                                        >  1:1950619/1‑150 (MQ=255)
                                                                                               cAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCg                                                                                                                                             >  1:6495461/1‑64 (MQ=255)
                                                                                                aGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATAc                                                      >  2:2927400/1‑150 (MQ=255)
                                                                                                aGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATAc                                                      >  2:2547585/1‑150 (MQ=255)
                                                                                                aGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATAc                                                      >  2:2259179/1‑150 (MQ=255)
                                                                                                  ccTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACgg                                                    >  2:312013/1‑150 (MQ=255)
                                                                                                  ccTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACgg                                                    >  1:4338310/1‑150 (MQ=255)
                                                                                                   cTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGgt                                                   >  2:5587767/1‑150 (MQ=255)
                                                                                                   cTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGgt                                                   >  2:4199112/1‑150 (MQ=255)
                                                                                                       aGCGCCAGGGCGCGGTAGTGGACGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTAc                                               >  1:4601810/1‑150 (MQ=255)
                                                                                                        gcgcCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACa                                              >  1:3193652/1‑150 (MQ=255)
                                                                                                           ccAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCTCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAAt                                           >  2:5769214/1‑150 (MQ=255)
                                                                                                           ccAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAAt                                           >  1:5530018/1‑150 (MQ=255)
                                                                                                             aGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCTCATGGCAATCTCCTTTCAATGTCAATTTCAGTCGACGAGTATCTTCGGGTTTGAATTGCCCCATACTGTGTACAAATCg                                         <  2:668926/150‑1 (MQ=255)
                                                                                                             aGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCg                                         <  2:5823601/150‑1 (MQ=255)
                                                                                                               ggCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGc                                       <  1:2772045/150‑1 (MQ=255)
                                                                                                                 cgcgGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACACATCGGCGa                                     <  2:5847044/150‑1 (MQ=255)
                                                                                                                        gTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAAcgc                              >  1:4135210/1‑150 (MQ=255)
                                                                                                                        gTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAAcgc                              >  2:4735803/1‑150 (MQ=255)
                                                                                                                         tGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAAcgcg                             >  2:2156351/1‑150 (MQ=255)
                                                                                                                           gATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGcc                           >  2:5386268/1‑150 (MQ=255)
                                                                                                                              gTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGtt                        <  2:3474025/150‑1 (MQ=255)
                                                                                                                              gTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGtt                        <  2:3288775/150‑1 (MQ=255)
                                                                                                                               tCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTg                       >  2:2958053/1‑150 (MQ=255)
                                                                                                                                 ccAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGAt                     <  2:6647655/150‑1 (MQ=255)
                                                                                                                                  cAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCAACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGAtt                    <  2:5451517/150‑1 (MQ=255)
                                                                                                                                   ccaacGCATCGCCGGCGACGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGAttt                   >  1:3127248/5‑150 (MQ=255)
                                                                                                                                    ggCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATttc                  <  2:360589/150‑1 (MQ=255)
                                                                                                                                     gCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATttct                 <  2:3412885/150‑1 (MQ=255)
                                                                                                                                     aacGCATCGCCGGCGACGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATttct                 >  1:2957527/3‑150 (MQ=255)
                                                                                                                                      ccGCAGCGCCGCCGCCGGCGATGCGTTGGGTGAGCTTGGCGGAGCGCATGGCAATCTCCTTTCACTGGCAATTTCAGTCGACGAGTATCTTCGGGCTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATttctt                <  2:5520051/150‑1 (MQ=255)
                                                                                                                                       cGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATttcttc               <  2:2025656/150‑1 (MQ=255)
                                                                                                                                       cGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATttcttc               <  1:6662596/150‑1 (MQ=255)
                                                                                                                                        gCATCGCCGGCGACGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTcc              >  2:5040058/1‑150 (MQ=255)
                                                                                                                                         cAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCa             >  2:322181/1‑150 (MQ=255)
                                                                                                                                             gcCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCAtctc         <  1:6555117/150‑1 (MQ=255)
                                                                                                                                              ccGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCg        <  1:2614966/150‑1 (MQ=255)
                                                                                                                                                  cGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCGGCCt    <  2:30647/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCGGCCTGG  >  NC_002947/4249831‑4250128

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