Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,259 C→T 6.3% intergenic (+287/‑233) 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,2590CT 6.3% 270.1 / 7.1 95intergenic (+287/‑233)PP_t18/PP_1703tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base C (48/40);  new base T (2/4);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

CAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTC  >  NC_002947/1899113‑1899407
                                                                                                                                                   |                                                                                                                                                    
cAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGCGTGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCC                                                                                                                                                   <  1:1155202‑M1/150‑1 (MQ=255)
  gCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCG                                                                                                                                                 >  1:1257618‑M1/1‑150 (MQ=255)
    aCCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCT                                                                                                                                               >  2:140732‑M1/1‑150 (MQ=255)
      cTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGTAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGATGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCC                                                                                                                                             >  2:141410‑M1/1‑150 (MQ=255)
      cTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCC                                                                                                                                             >  1:628427‑M1/1‑150 (MQ=255)
      cTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCC                                                                                                                                             >  2:1284444‑M1/1‑150 (MQ=255)
      cTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCATCGCGGGCACGCCCGCTCC                                                                                                                                             >  2:939116‑M1/1‑150 (MQ=255)
       tCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCC                                                                                                                                            >  2:139303‑M1/1‑150 (MQ=255)
       tCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCC                                                                                                                                             >  1:696262‑M1/1‑149 (MQ=255)
       tCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCC                                                                                                                                             <  2:696262‑M1/149‑1 (MQ=255)
          aaCCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACA                                                                                                                                         <  2:906315‑M1/150‑1 (MQ=255)
            cccGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCCTTTCGCGGGCATGCCCGCTCCCACAGG                                                                                                                                       <  1:82396‑M1/150‑1 (MQ=255)
                tCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTTCCCCAGCAAAGCGCCTGCGAGAATGATCCAGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCGCAGGTACt                                                                                                                                   >  1:324194/1‑150 (MQ=255)
                           cGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTc                                                                                                                        >  2:43270/1‑150 (MQ=255)
                           cGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTAACCTAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTc                                                                                                                        <  1:866240/150‑1 (MQ=255)
                            ggCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCa                                                                                                                       <  1:140732/150‑1 (MQ=255)
                               aTCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGc                                                                                                                    >  2:490721/1‑150 (MQ=255)
                                    cTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGgatatcgcaaggcctgagggcggc                                                                                                               >  2:1197924‑M1/1‑126 (MQ=255)
                                       ggACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAg                                                                                                                      >  1:1363373/1‑140 (MQ=255)
                                       ggACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAg                                                                                                                      <  2:1363373/140‑1 (MQ=255)
                                             gggAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGgatatcgcaaggcctgagggcggcgcagcacct                                                                                                      <  2:1232875‑M1/150‑34 (MQ=255)
                                               gAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgt                                                                                                    >  1:490721‑M2/1‑145 (MQ=255)
                                                 gggCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGgatatcgcaaggcctgagggcggcgcagtacctgtgg                                                                                                  <  1:1197924‑M1/150‑38 (MQ=255)
                                                 gggCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGgatatcgcaaggcctgagggcggcgcagtacctgtgg                                                                                                  <  1:391123‑M1/150‑38 (MQ=255)
                                                       tGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgtgggagcgg                                                                                            >  2:1257618‑M2/1‑137 (MQ=255)
                                                             tGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGc                                                                                      <  2:1099405/150‑1 (MQ=255)
                                                             tGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGc                                                                                      <  2:958911/150‑1 (MQ=255)
                                                             tGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGc                                                                                      >  2:150691/1‑150 (MQ=255)
                                                               cTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGG                                                                                                                                       >  2:1245984‑M1/1‑99 (MQ=255)
                                                               cTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGG                                                                                                                                       <  1:1245984‑M1/99‑1 (MQ=255)
                                                                tGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgtg                                                                                                   <  2:1216909‑M2/134‑7 (MQ=255)
                                                                tGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgtg                                                                                                   >  1:1216909‑M2/1‑128 (MQ=255)
                                                                      ttGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAg                                                                             <  1:1110947/150‑1 (MQ=255)
                                                                       gatctcAGCAAAGCGCCTGGGAGCATGAGGTTGACATGAATGGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAgg                                                                            <  1:1143163/145‑1 (MQ=38)
                                                                         ccccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGgatatcgcaaggcctgagggcggcgcagtacctgtgg                                                                                                  >  1:130988‑M1/1‑89 (MQ=255)
                                                                         ccccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGgatatcgcaaggcctgagggcggcgcagtacctgtgg                                                                                                  <  2:130988‑M1/126‑38 (MQ=255)
                                                                         ccccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGgatatcgcaaggcctgagggcggcgcagtacctgtggga                                                                                                <  2:1335035‑M1/128‑40 (MQ=255)
                                                                         ccccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGgatatcgcaaggcctgagggcggcgcagtacctgtggga                                                                                                >  1:1335035‑M1/1‑89 (MQ=255)
                                                                          cccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGcc                                                                         >  2:1391905/1‑150 (MQ=255)
                                                                          cccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGcc                                                                         >  1:781260/1‑150 (MQ=255)
                                                                           ccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGTTATCcctgtgggagcg                                                                                             <  1:680413‑M2/129‑13 (MQ=255)
                                                                           ccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGTTATCcctgtgggagcg                                                                                             >  2:680413‑M2/1‑117 (MQ=255)
                                                                           ccAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCCTTTCGCGGGCATGCCCGCTCCCACAGGgatatcgcaaggcctgagggcggcgcagtacctgtgggagcgggcgtgcccgcgaagaggcct                                                                        <  1:1114015‑M1/150‑64 (MQ=255)
                                                                            cAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgt                                                                                                    >  2:1349430‑M2/1‑116 (MQ=255)
                                                                            cAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCg                                                                        <  1:808435/149‑1 (MQ=255)
                                                                            cAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCg                                                                        >  2:808435/1‑149 (MQ=255)
                                                                            cAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCAACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgt                                                                                                    <  1:1349430‑M2/121‑6 (MQ=255)
                                                                             aGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCCTGCCGCGTAGCCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGt                                                                      <  2:1403982/150‑1 (MQ=255)
                                                                             aGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCAAGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGTAGCGGGCATTCCCGCTAAAgggagatc                                                                      >  2:1143163/1‑145 (MQ=39)
                                                                              gCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGta                                                                     >  1:516764/1‑150 (MQ=255)
                                                                                     gcCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTgg                                                              <  2:628427/150‑1 (MQ=255)
                                                                                       cTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgtgg                                                                                                  <  1:254176‑M2/112‑8 (MQ=255)
                                                                                       cTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgtgg                                                                                                  >  2:254176‑M2/1‑105 (MQ=255)
                                                                                        tGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGccc                                                           >  2:428168/1‑150 (MQ=255)
                                                                                         gggAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGc                                                                                      >  2:502939/1‑122 (MQ=255)
                                                                                         gggAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGc                                                                                      <  1:502939/122‑1 (MQ=255)
                                                                                           gAGCATGAGGTTGCCACCGCATTGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCACGCTCCCACAGGTACTTCCCCGCCCTCAGGCCTTGCAATATCACTGTGTTAGCGGGAATGCCCGCGAAAGGGACCGTGTAGCTGACCCATa                                                        >  2:346656/1‑150 (MQ=255)
                                                                                             gCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCGCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATaaa                                                      <  1:139303/150‑1 (MQ=255)
                                                                                               aTGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCGTGCCCGCgaa                                                                                                                                             <  2:776086‑M1/61‑4 (MQ=255)
                                                                                               aTGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCGTGCCCGCgaa                                                                                                                                             >  1:776086‑M1/1‑58 (MQ=255)
                                                                                                tGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTg                                                   >  1:1353443/1‑150 (MQ=255)
                                                                                                 gAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGaaa                                                                              >  2:818213/1‑122 (MQ=255)
                                                                                                 gAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGaaa                                                                              <  1:818213/122‑1 (MQ=255)
                                                                                                  aGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCt                                                 <  1:43270/150‑1 (MQ=255)
                                                                                                    gTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGa                                               >  2:493317/1‑150 (MQ=255)
                                                                                                    gTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGa                                               <  1:493317/150‑1 (MQ=255)
                                                                                                    gTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGa                                               <  2:634995/150‑1 (MQ=255)
                                                                                                    gTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCAGGTATCGCTGGCCCATAAATTGCCGa                                               <  2:1157334/150‑1 (MQ=255)
                                                                                                      tGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATc                                             >  2:970239/1‑150 (MQ=255)
                                                                                                         cACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTAc                                          >  2:208737/1‑150 (MQ=255)
                                                                                                         cACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTAc                                          >  2:1300942/1‑150 (MQ=255)
                                                                                                           ccGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCt                                        >  2:1281238/1‑150 (MQ=255)
                                                                                                            cGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAG                                                                                                                                        <  2:690632‑M1/53‑1 (MQ=255)
                                                                                                            cGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAG                                                                                                                                        >  1:690632‑M1/1‑53 (MQ=255)
                                                                                                                         aaCCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACTGCGCCGCGCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGt                           >  1:355134/1‑150 (MQ=255)
                                                                                                                          aCCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGcc                                                                         >  1:669775/1‑103 (MQ=255)
                                                                                                                          aCCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGcc                                                                         <  2:669775/103‑1 (MQ=255)
                                                                                                                            cTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCaa                        >  1:457690/1‑150 (MQ=255)
                                                                                                                                  aGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCcctgtgggagcgggcatgcccgcgaaagggcctgcacaggttacgcaccatgcggtggcaacctcatgctcccaggcgctttgctgg                  <  1:1114015‑M2/150‑88 (MQ=255)
                                                                                                                                   ggCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCCCCGTCGAc                 >  1:922920/1‑150 (MQ=255)
                                                                                                                                   ggCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGAc                 <  1:150691/150‑1 (MQ=255)
                                                                                                                                    gCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACa                >  2:701827/1‑150 (MQ=255)
                                                                                                                                      ctctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCt                                                 <  1:1052790/115‑1 (MQ=255)
                                                                                                                                      ctctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCt                                                 >  2:1052790/1‑115 (MQ=255)
                                                                                                                                       tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATGCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc             >  1:1243524/1‑150 (MQ=255)
                                                                                                                                       tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGcc                                                                                     >  2:1259795/1‑78 (MQ=255)
                                                                                                                                       tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGcc                                                                                     <  1:1259795/78‑1 (MQ=255)
                                                                                                                                       tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc             <  1:939116/150‑1 (MQ=255)
                                                                                                                                       tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc             <  1:208737/150‑1 (MQ=255)
                                                                                                                                       tctTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGcc             >  1:955013/1‑150 (MQ=255)
                                                                                                                                           cgcgGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTca         <  2:516764/150‑1 (MQ=255)
                                                                                                                                            gcgGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGc                                                                                                               >  1:213913/1‑47 (MQ=255)
                                                                                                                                            gcgGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGc                                                                                                               <  2:213913/47‑1 (MQ=255)
                                                                                                                                            gcgGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTGGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTcac        <  1:1284444/150‑1 (MQ=255)
                                                                                                                                                  aCGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGTCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTc  >  1:249182/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCATGCTGGAAATTGCCCCAGCAAAGCGCCTGGGAGCATGAGGTTGCCACCGCATGGTGCGTAACCTGTGCAGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCCGCCCTCAGGCCTTGCGATATCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCGGTATAGCTGGCCCATAAATTGCTGATCTACCTCATCCCGGGCTGTCAACGTCGACAGCCCTCACCACTTC  >  NC_002947/1899113‑1899407

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