Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,028 G→A 7.1% intergenic (+56/‑464) 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,0280GA 7.1% 248.6 / 7.5 85intergenic (+56/‑464)PP_t18/PP_1703tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base G (43/36);  new base A (3/3);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01

ACCACTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCT  >  NC_002947/1898879‑1899177
                                                                                                                                                     |                                                                                                                                                     
acctggcGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATCCCAAAGAGATACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAg                                                                                                                                                       >  1:1870714/7‑150 (MQ=11)
accacTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAg                                                                                                                                                       >  2:2386654/1‑150 (MQ=32)
accacTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCAGAAGTAAAAAAGAAAGCGACCTTAg                                                                                                                                                       >  1:1059478/1‑150 (MQ=17)
   acTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAgg                                                                                                                                                      <  2:133596/148‑1 (MQ=33)
   acTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAgg                                                                                                                                                      >  1:133596/1‑148 (MQ=33)
    cTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTc                                                                                                                                                   <  1:2399481/150‑1 (MQ=35)
       gATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCt                                                                                                                                                >  1:537431/1‑150 (MQ=35)
          gTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCtttt                                                                                                                                             <  1:1357304/150‑1 (MQ=35)
           tataGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCATGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCttttt                                                                                                                                            <  2:1737390/150‑1 (MQ=25)
             taGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAAAAAGCGACCCTAAGGTCGctttcttt                                                                                                                                          <  1:2035927/150‑1 (MQ=11)
               gcgcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTAt                                                                                                                                        <  1:77569/150‑1 (MQ=37)
                cgcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCATGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTAtt                                                                                                                                       >  2:2126254/1‑150 (MQ=255)
                 gcAGCGGTAGTTCAGTCGGTTCGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCAAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGCCCTGAGGGTCGCTTTTTTTAttt                                                                                                                                      >  1:178853/1‑150 (MQ=255)
                    gCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCATGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTg                                                                                                                                   <  1:2126254/150‑1 (MQ=255)
                     cGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGc                                                                                                                                  <  2:1974372/150‑1 (MQ=255)
                       gTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGcat                                                                                                                                >  1:844790/1‑150 (MQ=255)
                          gTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGcatcat                                                                                                                             <  1:987531/150‑1 (MQ=255)
                           ttCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTAttt                                                                                                                                      >  2:101977/1‑140 (MQ=40)
                           ttCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTAttt                                                                                                                                      <  1:101977/140‑1 (MQ=39)
                             aaGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTCAGCACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCa                                                                                                                          <  1:416353/149‑1 (MQ=255)
                                          tACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTGTCTTTTTTTATTTCTGCATCATGCAATCGTTTACGccc                                                                                                             >  1:129077/1‑150 (MQ=255)
                                               gCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCt                                                                                                            >  1:2193709/1‑146 (MQ=255)
                                               gCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCt                                                                                                            <  2:2193709/146‑1 (MQ=255)
                                                 cTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGcaaca                                                                                                      >  1:2229495/1‑150 (MQ=255)
                                                  tGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAAAAAGCGACCCTAAGGTCGctttcttt                                                                                                                                          >  2:459740/1‑113 (MQ=11)
                                                  tGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAAAAAGCGACCCTAAGGTCGctttcttt                                                                                                                                          <  1:459740/113‑1 (MQ=11)
                                                    tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAgg                                                                                                   >  1:2202454/1‑150 (MQ=255)
                                                      aCGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAgggg                                                                                                 <  1:352812/150‑1 (MQ=255)
                                                         ccGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGca                                                                                                         >  2:123158/1‑139 (MQ=255)
                                                         ccGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGca                                                                                                         <  1:123158/139‑1 (MQ=255)
                                                          cGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAg                                                                                             >  1:362824/1‑150 (MQ=255)
                                                           gggggTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTa                                                                                                                                         >  2:1146079/1‑105 (MQ=38)
                                                           gggggTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTa                                                                                                                                         <  1:1146079/105‑1 (MQ=38)
                                                              ggTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTgc                                                                                         <  1:244514/150‑1 (MQ=255)
                                                              ggTCGCGGGTTCGAGTCCCGTCCGCTGAGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAAAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTgc                                                                                         <  1:1110163/150‑1 (MQ=255)
                                                                tCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTCGGGTCGCCTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTCCAACCAGCGGCTAGATgcgc                                                                                       >  2:232423/1‑150 (MQ=18)
                                                                  gcgGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCa                                                                                     >  2:666430/1‑150 (MQ=255)
                                                                       ttCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACtctt                                                                                >  1:86296/1‑150 (MQ=255)
                                                                        tCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACtcttc                                                                               <  1:1437097/150‑1 (MQ=255)
                                                                        tCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACtcttc                                                                               <  1:1439115/150‑1 (MQ=255)
                                                                          gAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCa                                                                                     >  2:2295706/1‑142 (MQ=255)
                                                                          gAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCa                                                                                     <  1:2295706/142‑1 (MQ=255)
                                                                          gAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTcttctt                                                                             >  2:322450/1‑150 (MQ=255)
                                                                               ccGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAAGACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCAAGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGc                                                                        >  1:2303402/1‑150 (MQ=255)
                                                                               ccGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTAt                                                                                                                                        <  1:1157664/86‑1 (MQ=255)
                                                                               ccGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTAt                                                                                                                                        >  2:1157664/1‑86 (MQ=255)
                                                                                cGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTgcgc                                                                                       >  1:1435984/1‑134 (MQ=255)
                                                                                cGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTgcgc                                                                                       <  2:1435984/134‑1 (MQ=255)
                                                                                     gCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTc                                                                  <  2:844790/150‑1 (MQ=255)
                                                                                     gCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTc                                                                  >  1:185572/1‑150 (MQ=255)
                                                                                          gcCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAAGAAGGGGCTCGCTGCGCCACTCTTCTTCCCGCAGGGTCAGCAc                                                             >  1:526054/1‑150 (MQ=255)
                                                                                          gcCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCAc                                                             <  2:185572/150‑1 (MQ=255)
                                                                                             atatTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTc                                                          >  1:733256/1‑150 (MQ=255)
                                                                                                 tCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAAc                                                      >  1:514367/1‑150 (MQ=255)
                                                                                                 tCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAAc                                                      >  1:512942/1‑150 (MQ=255)
                                                                                                      ttCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAgggg                                                                                                 >  2:99895/1‑102 (MQ=255)
                                                                                                      ttCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAgggg                                                                                                 <  1:99895/102‑1 (MQ=255)
                                                                                                      ttCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACtcttc                                                                               <  2:813779/120‑1 (MQ=255)
                                                                                                      ttCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACtcttc                                                                               >  1:813779/1‑120 (MQ=255)
                                                                                                      ttCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTc                                                 >  2:617503/1‑150 (MQ=255)
                                                                                                       tCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTcc                                                <  2:2229495/150‑1 (MQ=255)
                                                                                                        ccAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCt                                               >  1:2392241/1‑150 (MQ=255)
                                                                                                            ggCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGAt                                           >  1:779662/1‑150 (MQ=255)
                                                                                                               ccACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACAATGCAGTCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCAc                                        >  1:1832550/1‑150 (MQ=38)
                                                                                                                      acacCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTCGCTGCGCCCCTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCTTCCTCGATCACCGTCAtc                                 <  2:2202454/150‑1 (MQ=255)
                                                                                                                      acacCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCAtc                                 >  1:1745274/1‑150 (MQ=255)
                                                                                                                        acCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCAtctc                               >  1:2231369/1‑150 (MQ=255)
                                                                                                                          ccGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGcaac                                                                                                       >  1:129063/1‑76 (MQ=255)
                                                                                                                          ccGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGcaac                                                                                                       <  2:129063/76‑1 (MQ=255)
                                                                                                                          ccGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCt                             >  1:835728/1‑150 (MQ=255)
                                                                                                                          ccGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCt                             <  2:733256/150‑1 (MQ=255)
                                                                                                                          ccGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCt                             >  1:1756068/1‑150 (MQ=255)
                                                                                                                                  aaaaaaGAAAGCGACCTTAGGGTCGTTTTTTTTATTTATGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATg                     <  2:537431/150‑1 (MQ=255)
                                                                                                                                  aaaaaaGAAAGCGACCGTAGGGTCGCTTTTTTTATTTCTGCATAATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATg                     >  2:1677053/1‑150 (MQ=255)
                                                                                                                                   aaaaaGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATgg                    >  2:1111568/1‑150 (MQ=255)
                                                                                                                                   aaaaaGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTAAAGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATgg                    >  1:297746/1‑150 (MQ=255)
                                                                                                                                    aaaaGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCa                                                              >  1:1107632/1‑107 (MQ=255)
                                                                                                                                    aaaaGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCa                                                              <  2:1107632/107‑1 (MQ=255)
                                                                                                                                       agaaagCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCt                             >  1:1362246/1‑137 (MQ=255)
                                                                                                                                       agaaagCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCt                             <  2:1362246/137‑1 (MQ=255)
                                                                                                                                        aaaagCGACCCTAAGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAAc                                                      <  2:231396/110‑1 (MQ=255)
                                                                                                                                        aaaagCGACCCTAAGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAAc                                                      >  1:231396/2‑111 (MQ=255)
                                                                                                                                        aaaagCGACCCTAAGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAggg               >  2:1740199/2‑150 (MQ=255)
                                                                                                                                         aaagCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGc              >  1:849718/1‑150 (MQ=255)
                                                                                                                                          aagCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCa             <  2:2392241/150‑1 (MQ=255)
                                                                                                                                            gCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTTCGCCCTGCAACGAGGGGCTAGCTGCCCCCCTCTTCTTCCCGCAGTGTCAGCACCTCAAACCCGGCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGt           <  2:178853/150‑1 (MQ=255)
                                                                                                                                                 cTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGc                                                                     >  1:1356807/1‑87 (MQ=255)
                                                                                                                                                 cTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGc                                                                     <  2:1356807/87‑1 (MQ=255)
                                                                                                                                                 cTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCa      <  2:512942/150‑1 (MQ=255)
                                                                                                                                                 cTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCa      <  2:514367/150‑1 (MQ=255)
                                                                                                                                                     gggTCGCTTTTTTTATTTTTTTTTTTTGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTACATTATACCAGCCTCCTTGGACATGGGAGGGCAGTGCTCATGCt  >  2:195081/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCACTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCT  >  NC_002947/1898879‑1899177

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: