Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,342 A→G 6.8% intergenic (+412/‑108) PP_t18 → / → PP_1703 tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,899,3420AG 6.8% 168.5 / 2.5 59intergenic (+412/‑108)PP_t18/PP_1703tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (26/29);  new base G (2/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGCCCCTACTGCGG  >  NC_002947/1899197‑1899490
                                                                                                                                                 |                                                                                                                                                      
tgCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAAc                                                                                                                                                      <  1:318152/148‑1 (MQ=255)
         ctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGt                                                                                                                                                    <  3:171572/141‑1 (MQ=255)
         ctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATATACCTAATACCGGGCTGTCAACGt                                                                                                                                                    >  4:171572/1‑141 (MQ=255)
                ggCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTcacca                                                                                                                                    <  8:571/149‑1 (MQ=255)
                  cACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACATCATCCCGGGCTGTCAACGTCGa                                                                                                                                                 >  6:288684/1‑135 (MQ=255)
                  cACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACATCATCCCGGGCTGTCAACGTCGa                                                                                                                                                 <  5:288684/135‑1 (MQ=255)
                     gCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATGTCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCa                                                                                                                               >  2:210283/1‑149 (MQ=255)
                     gCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGGATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCa                                                                                                                               >  4:333822/1‑149 (MQ=255)
                             ccACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAg                                                                                                                       >  2:100273/1‑149 (MQ=255)
                             ccACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAg                                                                                                                       >  4:77155/1‑149 (MQ=255)
                               acaGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCGACGTTGACAGcccg                                                                                                                                          >  8:323915/1‑128 (MQ=255)
                               acaGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCGACGTTGACAGcccg                                                                                                                                          <  7:323915/129‑2 (MQ=255)
                                   gTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc                                                                                                                                            >  8:213093/1‑123 (MQ=255)
                                   gTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc                                                                                                                                            <  7:213093/123‑1 (MQ=255)
                                    tACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCgt                                                                                                                <  3:77155/149‑1 (MQ=255)
                                     aCTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCAc                                                                                                                              <  5:138605/134‑1 (MQ=255)
                                     aCTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCAc                                                                                                                              >  6:138605/1‑134 (MQ=255)
                                      cTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCgtgt                                                                                                              <  4:172582/149‑1 (MQ=255)
                                         cgccgcCCTCAGGCCTTGCTATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACaaa                                                                                                                        <  7:245896/136‑1 (MQ=255)
                                         cgccgcCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCa                                                                                                                               <  7:302600/129‑1 (MQ=255)
                                         cgccccCCTCAGGCCTTGCGATATCACTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGATGGCCCATAAATTGCTGATCTACCTCATCCAGGGCTGTCAACGTCGAAAGACCTACCCACTTCACGACaaa                                                                                                                        >  8:245896/1‑136 (MQ=38)
                                                cTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCtt                                                                                                   <  2:55599/149‑1 (MQ=255)
                                                     gCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCCCCCCTTCACGACAAAGGCGCCGTGTTGCCCCGCGTGCAt                                                                                               >  8:193600/1‑148 (MQ=255)
                                                          gCGATATCCCTGTGGTTGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGc                                                                                           <  2:82510/147‑1 (MQ=255)
                                                          gCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGc                                                                                           <  5:266960/147‑1 (MQ=255)
                                                               atCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGgc                                                                                                                     <  2:41552/117‑1 (MQ=255)
                                                               atCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGgc                                                                                                                     >  1:41552/1‑117 (MQ=255)
                                                                tCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCggg                                                                                   >  3:342296/1‑149 (MQ=255)
                                                                  gcTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCa                                                                                             <  8:323145/136‑1 (MQ=255)
                                                                  ccTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGc                                                                                 <  1:210283/149‑1 (MQ=255)
                                                                      tgGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGc                                                                              <  3:333822/148‑1 (MQ=255)
                                                                               gCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGtt                                                                                                             >  3:293374/1‑109 (MQ=255)
                                                                               gCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGtt                                                                                                             <  4:293374/109‑1 (MQ=255)
                                                                               gCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCttt                                                                    >  2:246342/1‑149 (MQ=255)
                                                                                    cccGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCGGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGc                                                                <  7:193600/148‑1 (MQ=255)
                                                                                         gAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTc                                                          <  4:84758/149‑1 (MQ=255)
                                                                                                ccGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCACCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGa                                                   <  1:248578/149‑1 (MQ=255)
                                                                                                            gCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTGACCACTTCACGACAAAGGCGCCGTGTTGCCCCGACTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATg                                        <  1:100273/149‑1 (MQ=255)
                                                                                                            gCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATg                                        <  4:342296/148‑1 (MQ=255)
                                                                                                              ccATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGc                                                                                           >  4:78799/1‑95 (MQ=255)
                                                                                                              ccATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGc                                      >  5:256531/1‑148 (MQ=255)
                                                                                                              ccATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGc                                      >  3:113949/1‑148 (MQ=255)
                                                                                                              ccATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGc                                      >  3:144660/1‑148 (MQ=255)
                                                                                                                aTAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCaa                                   >  2:194003/1‑149 (MQ=255)
                                                                                                                aTAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCa                                    >  8:334386/1‑148 (MQ=255)
                                                                                                                        cTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGaa                            <  1:194003/148‑1 (MQ=255)
                                                                                                                          gATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTa                         >  4:66417/1‑149 (MQ=255)
                                                                                                                           aTCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTAc                        >  2:110383/1‑149 (MQ=255)
                                                                                                                                ccTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCa                     <  7:334386/147‑1 (MQ=255)
                                                                                                                                ccTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAg                    <  6:256531/148‑1 (MQ=255)
                                                                                                                                 cTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGa                             >  4:343831/1‑138 (MQ=255)
                                                                                                                                 cTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGa                             <  3:343831/138‑1 (MQ=255)
                                                                                                                                   cATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTc                >  4:217294/1‑149 (MQ=255)
                                                                                                                                     tCCCGGGCTGTCAACGTCGACAGCCCTTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGGCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCt               >  1:191464/1‑149 (MQ=255)
                                                                                                                                      tgagggCTGTCGACGTTGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCa                     <  6:96630/138‑1 (MQ=255)
                                                                                                                                      tgagggCTGTCGACGTTGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCa                     >  5:96630/4‑141 (MQ=255)
                                                                                                                                       ccGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCTTGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGc             <  1:110383/148‑1 (MQ=255)
                                                                                                                                               tCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGCCCCTACTgc    <  6:112246/149‑1 (MQ=255)
                                                                                                                                               tCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGCCCCTACTgc    >  5:112246/1‑149 (MQ=255)
                                                                                                                                                 aaCGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATg                                        <  5:241938/111‑1 (MQ=255)
                                                                                                                                                 aaCGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATg                                        >  6:241938/1‑111 (MQ=255)
                                                                                                                                                 aaCGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGc             <  3:204383/138‑1 (MQ=255)
                                                                                                                                                 aaCGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGc             >  4:204383/1‑138 (MQ=255)
                                                                                                                                                 aaCGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGCCCCTACTgcgg  >  2:275134/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
TGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTCACGACAAAGGCGCCGTGTTGCCCCGCTTGCATCAGCATGCCGGGGCAGCCGGCGCCTTTTTGCGTTTTCCACAGGAGCCCGCCCATGGCCAACAGCGAAGTACGCAGCGTCTGCCCCTACTGCGG  >  NC_002947/1899197‑1899490

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