Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,030,463 G→C 9.1% intergenic (+420/‑274) pcaY → / → mgtA aromatic acid chemoreceptor/ATP‑dependent magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,030,4630GC 9.1% 220.1 / 11.8 74intergenic (+420/‑274)pcaY/mgtAaromatic acid chemoreceptor/ATP‑dependent magnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (27/40);  new base C (4/3);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

TGTTGGCCGGTTGGAGGCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAA  >  NC_002947/3030321‑3030607
                                                                                                                                                |                                                                                                                                                
tgtTGGCCGGTTGGAGGCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCggg                                                                                                                                           <  1:2159830/150‑1 (MQ=255)
             gAGGCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCgg                                                                                                                                            >  2:1925956/1‑136 (MQ=255)
             gAGGCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCgg                                                                                                                                            <  1:1925956/136‑1 (MQ=255)
              aGGCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCt                                                                                                                             >  1:2265985/1‑150 (MQ=255)
              aGGCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCt                                                                                                                             <  2:1685804/150‑1 (MQ=255)
                gCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTc                                                                                                                           <  2:1720771/150‑1 (MQ=255)
                  aaGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCAt                                                                                                                         <  1:793134/150‑1 (MQ=255)
                  aaGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCAt                                                                                                                         <  1:1859803/150‑1 (MQ=255)
                   aGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATg                                                                                                                        <  2:800415/150‑1 (MQ=255)
                   aGCTCATGAACCTCACTTCATTCAATTCACATCGCCAGCCGTCGTGCCTTGCGCCGCGGCAGGCCACTCAGCTCGCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGATGGCGCGTCTTCATg                                                                                                                        <  1:2156123/150‑1 (MQ=255)
                      tCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCGCCAGCACCCCGCAACTACTTGGGGGGTGCt                                                                                                                                      >  1:2297743/1‑133 (MQ=39)
                      tCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCGCCAGCACCCCGCAACTACTTGGGGGGTGCt                                                                                                                                      <  2:2297743/133‑1 (MQ=39)
                          gAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCtt                                                                                                                 <  2:836148/150‑1 (MQ=255)
                          gAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCtt                                                                                                                 >  1:1944179/1‑150 (MQ=255)
                            aCCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGa                                                                                                               <  2:1944179/150‑1 (MQ=255)
                            aCCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGa                                                                                                               <  1:549127/150‑1 (MQ=255)
                             ccTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCGCAACTACTTGGGGGGTGCTgg                                                                                                                                    <  1:86832/128‑1 (MQ=39)
                             ccTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCGCAACTACTTGGGGGGTGCTgg                                                                                                                                    >  2:86832/1‑128 (MQ=39)
                             ccTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAg                                                                                                              >  2:1211472/1‑150 (MQ=255)
                              cTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGc                                                                                                             >  2:1418095/1‑150 (MQ=255)
                                cACTTCATTCGTTTCATATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGccc                                                                                                           >  2:256642/1‑150 (MQ=255)
                                  cTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGccccc                                                                                                         >  1:200728/1‑150 (MQ=255)
                                              cacaTCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCgg                                                                                                                                            >  2:780997/1‑103 (MQ=255)
                                              cacaTCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCgg                                                                                                                                            <  1:780997/103‑1 (MQ=255)
                                              cacaTCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACt                                                                                             <  1:305682/150‑1 (MQ=255)
                                               acaTCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACt                                                                                             >  2:692872/1‑149 (MQ=255)
                                               acaTCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACt                                                                                             <  1:692872/149‑1 (MQ=255)
                                                 aTCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGaa                                                                                          >  2:1529534/1‑150 (MQ=255)
                                                    gCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCa                                                                                                                      >  1:1883812/1‑119 (MQ=255)
                                                    gCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCa                                                                                                                      <  2:1883812/119‑1 (MQ=255)
                                                      cAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATaaaa                                                                                     <  2:2011912/150‑1 (MQ=255)
                                                       aGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAg                                                                                    >  1:139797/1‑150 (MQ=255)
                                                           tGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATaaaa                                                                                     <  2:1373663/145‑1 (MQ=255)
                                                           tGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATaaaa                                                                                     >  1:1373663/1‑145 (MQ=255)
                                                               ggCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCTCCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTCAATAAATGATATTGaa                                                                            <  2:580164/149‑1 (MQ=255)
                                                                tgcagacgccgcGGCCAGGCGCGGGCCTGCCGCGGCGCCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAAc                                                                           >  2:1953695/7‑150 (MQ=38)
                                                                gCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAAc                                                                           <  1:256642/150‑1 (MQ=255)
                                                                gCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAAc                                                                           <  2:2170598/150‑1 (MQ=255)
                                                                   ttGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGat                                                                                  >  2:1773550/1‑140 (MQ=255)
                                                                   ttGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGat                                                                                  <  1:1773550/140‑1 (MQ=255)
                                                                     gcgcCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGACTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGACCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTAt                                                                      >  1:1463294/1‑150 (MQ=255)
                                                                         cgcgGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGTCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCttt                                                                  <  1:823644/150‑1 (MQ=255)
                                                                         cgcgGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCttt                                                                  <  2:200728/150‑1 (MQ=255)
                                                                         cgcgGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCttt                                                                  <  1:2076401/150‑1 (MQ=255)
                                                                                  gCCCGCGCCTTGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGa                                                          >  2:1059506/1‑150 (MQ=255)
                                                                                    ccGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAAc                                                        >  1:1720164/1‑150 (MQ=255)
                                                                                          ccttgCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTc                                                   <  1:1059506/146‑1 (MQ=255)
                                                                                          ccTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTc                                                   <  2:1720164/150‑1 (MQ=255)
                                                                                           cctgCCGCGGCGCCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCa                                                  <  1:1953695/147‑1 (MQ=255)
                                                                                           cTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCtt                                                                                                                 >  1:1967755/1‑87 (MQ=255)
                                                                                           cTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCtt                                                                                                                 <  2:1967755/87‑1 (MQ=255)
                                                                                                  cggcgTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTaaa                                           >  2:219830/1‑150 (MQ=255)
                                                                                                  cggcgTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTaaa                                           <  1:219830/150‑1 (MQ=255)
                                                                                                     cgTCTGCACGACTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGACCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGt                                        <  2:1463294/150‑1 (MQ=255)
                                                                                                         tGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTcgc                                    <  1:1368747/150‑1 (MQ=255)
                                                                                                         tGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATACTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTcgc                                    >  1:225666/1‑150 (MQ=255)
                                                                                                         tGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGATGGCGCGTCTTCATGCACGCTTGAGCCCCCGCAATTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTcgc                                    >  1:757910/1‑150 (MQ=255)
                                                                                                             cGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGa                                <  2:1944219/150‑1 (MQ=255)
                                                                                                                 tACCTTGACCACCCCCAGCCCTCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAAc                            <  2:139797/150‑1 (MQ=255)
                                                                                                                 tACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAAc                            >  1:1240871/1‑150 (MQ=255)
                                                                                                                    cTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTAt                         >  1:676758/1‑150 (MQ=255)
                                                                                                                        accaccaccAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTcc                     <  2:1190359/150‑1 (MQ=255)
                                                                                                                         ccaccaccaGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATaaaa                                                                                     >  2:1491914/1‑85 (MQ=255)
                                                                                                                         ccaccaccaGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATaaaa                                                                                     <  1:1491914/85‑1 (MQ=255)
                                                                                                                             cacctGCACCCCAAAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTTAGCCCCCGCACTTCTCACTAAATAAAAGATATTGAACGGTATCTTTGTACTCGAATTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCtgtg                <  1:2231296/150‑1 (MQ=255)
                                                                                                                             caccaGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTAACATGAAAAGTTcga                                    <  2:1738096/130‑2 (MQ=255)
                                                                                                                              gccaGCACCCCGCAACTACTTGGGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTgtgt               >  2:417363/2‑150 (MQ=255)
                                                                                                                               ccaGCACCCCCCAAGTAGTTGCGGGGTACTTTCGCGTCTTCATGCACGATTGCGCCCCCGCAGCTGGCACGGAAGAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGtt              <  1:1211472/150‑1 (MQ=255)
                                                                                                                                caGCACCCCGCAACTACTTGGGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTcc                     >  2:222873/1‑142 (MQ=255)
                                                                                                                                caGCACCCCGCAACTACTTGGGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTcc                     <  1:222873/142‑1 (MQ=255)
                                                                                                                                 aGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTa            <  2:2265985/150‑1 (MQ=255)
                                                                                                                                      cccccAAGTAGTTGCGGGCTGCTGGCGAGTCTTCATGCACGCTCACGCGATGGCATTTGTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAATTCGCGGGATAACTATTTCTATGTGTGTATGaaa       >  2:911040/1‑150 (MQ=255)
                                                                                                                                        cccAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGAACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAACGTTCGCGGGATAACTATTTCCCTGTGTTTCTGAAAAg     <  1:1418095/150‑1 (MQ=255)
                                                                                                                                           aaGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTAATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGaa  >  1:317472/1‑150 (MQ=255)
                                                                                                                                            agtagtTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTaaaa                                          >  2:703672/1‑109 (MQ=255)
                                                                                                                                            agtagtTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTaaaa                                          <  1:703672/109‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TGTTGGCCGGTTGGAGGCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAA  >  NC_002947/3030321‑3030607

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