Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,477,985 C→T 5.4% T189T (ACG→ACA)  nadD ← nicotinate‑nucleotide adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,477,9850CT 5.4% 904.6 / 25.9 279T189T (ACG→ACA) nadDnicotinate‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (149/115);  new base T (13/2);  total (162/117)
Fisher's exact test for biased strand distribution p-value = 2.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 290 total aligned reads displayed.

GGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGATGTGACAGTGTTGCAG  >  NC_002947/5477837‑5478134
                                                                                                                                                    |                                                                                                                                                     
 gATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGt                                                                                                                                                     >  1:4749254/1‑150 (MQ=255)
 gATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGt                                                                                                                                                     >  1:5534686/1‑150 (MQ=255)
 gATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGt                                                                                                                                                     <  2:6255603/150‑1 (MQ=255)
   tGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCt                                                                                                                                                   >  1:4226836/1‑150 (MQ=255)
   tGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCt                                                                                                                                                   >  1:568335/1‑150 (MQ=255)
    gAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTg                                                                                                                                                  >  1:2630250/1‑150 (MQ=255)
     aGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGc                                                                                                                                                 <  2:6303332/150‑1 (MQ=255)
      ggCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGcc                                                                                                                                                >  1:3696291/1‑150 (MQ=255)
      ggCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGcc                                                                                                                                                >  2:1592278/1‑150 (MQ=255)
       gcgcCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCa                                                                                                                                               <  1:5430298/150‑1 (MQ=255)
       gcgcCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCa                                                                                                                                               <  1:2508206/150‑1 (MQ=255)
       gcgcCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCa                                                                                                                                               >  1:17800/1‑150 (MQ=255)
          ccAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGac                                                                                                                                            <  2:7105562/150‑1 (MQ=255)
                    cGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAata                                                                                                                                  <  1:2259727/150‑1 (MQ=255)
                    cGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAata                                                                                                                                  <  2:4197435/150‑1 (MQ=255)
                      tCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGACAGCCAGACGAACGAAATAtt                                                                                                                                >  1:6018139/1‑150 (MQ=255)
                       cAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATAttt                                                                                                                               <  2:305975/150‑1 (MQ=255)
                       cAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCCGCCAGACGAACGAAATAttt                                                                                                                               >  1:2877383/1‑150 (MQ=255)
                            aggCGCACGAACTATATCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCccgc                                                                                                                          >  2:3309897/2‑150 (MQ=255)
                                gcACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCgccggg                                                                                                                      >  1:7401610/1‑150 (MQ=255)
                                  aCGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCgccgggcc                                                                                                                    <  1:443570/150‑1 (MQ=255)
                                    gATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATAttt                                                                                                                               <  1:5030271/137‑1 (MQ=255)
                                    gATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATAttt                                                                                                                               >  2:5030271/1‑137 (MQ=255)
                                     atatAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGa                                                                                                                 <  1:2929620/150‑1 (MQ=255)
                                     atatAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGa                                                                                                                 <  1:5977187/150‑1 (MQ=255)
                                         aGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATg                                                                                                             >  1:7086334/1‑150 (MQ=255)
                                                gcgcCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTgg                                                                                                      <  2:2208498/150‑1 (MQ=255)
                                                 cgcCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTggg                                                                                                     <  2:5887991/150‑1 (MQ=255)
                                                 cgcCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTggg                                                                                                     <  2:6001128/150‑1 (MQ=255)
                                                  gcCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGa                                                                                                    >  1:3052417/1‑150 (MQ=255)
                                                   ccTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGAt                                                                                                   >  2:1465184/1‑150 (MQ=255)
                                                   ccTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGAt                                                                                                   <  1:5259791/150‑1 (MQ=255)
                                                    cacGAGATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATc                                                                                                  >  2:3445879/3‑150 (MQ=255)
                                                    cTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATc                                                                                                  >  1:7163622/1‑150 (MQ=255)
                                                     tCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCg                                                                                                 >  1:6433938/1‑150 (MQ=255)
                                                     tCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCg                                                                                                 >  1:3009116/1‑150 (MQ=255)
                                                      cGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGc                                                                                                <  2:4527589/150‑1 (MQ=255)
                                                          ataGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCtctc                                                                                            <  1:5120907/150‑1 (MQ=255)
                                                          ataGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCtctc                                                                                            <  2:5199579/150‑1 (MQ=255)
                                                          ataGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCtctc                                                                                            <  2:1343544/150‑1 (MQ=255)
                                                             ggCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGa                                                                                         <  2:4447892/150‑1 (MQ=255)
                                                              gCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAg                                                                                        >  2:2081989/1‑150 (MQ=255)
                                                              gCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAg                                                                                        >  2:3662868/1‑150 (MQ=255)
                                                              gCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAg                                                                                        >  1:4883865/1‑150 (MQ=255)
                                                                cAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAgcg                                                                                      >  2:3123281/1‑150 (MQ=255)
                                                                    aCGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCa                                                                                  <  2:7326862/150‑1 (MQ=255)
                                                                    aCGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCa                                                                                  >  1:1047798/1‑150 (MQ=255)
                                                                      ggCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGc                                                                                >  2:3815902/1‑150 (MQ=255)
                                                                       gCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGcc                                                                               <  1:120777/150‑1 (MQ=255)
                                                                         gTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCaa                                                                             <  2:899559/150‑1 (MQ=255)
                                                                         gTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCaa                                                                             <  1:6907439/150‑1 (MQ=255)
                                                                         gTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCaa                                                                             <  1:4312428/150‑1 (MQ=255)
                                                                            cGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAg                                                                          >  2:7440383/1‑150 (MQ=255)
                                                                            cGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAg                                                                          >  2:2922474/1‑150 (MQ=255)
                                                                              gCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAgcg                                                                                      >  2:1301720/1‑136 (MQ=255)
                                                                                       aCCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGGTGCGCCACt                                                               >  1:7210658/1‑150 (MQ=255)
                                                                                         cTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCg                                                             >  2:534451/1‑150 (MQ=255)
                                                                                          tCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGt                                                            >  2:5253938/1‑150 (MQ=255)
                                                                                          tCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGt                                                            >  1:6167137/1‑150 (MQ=255)
                                                                                          tCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGt                                                            <  2:6642800/150‑1 (MQ=255)
                                                                                          tCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGt                                                            >  2:5614362/1‑150 (MQ=255)
                                                                                          tCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGt                                                            >  2:6062917/1‑150 (MQ=255)
                                                                                          tCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGGTGCGCAACTCGt                                                            >  1:3588846/1‑150 (MQ=255)
                                                                                           cACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTc                                                           >  2:4780490/1‑150 (MQ=255)
                                                                                           cACCGATTTGCCGCTGGCCACCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGAACGCTCTCCGAGCGCGCACCCAACAGGTTGCGCAACTCgac                                                           <  2:1138666/150‑3 (MQ=255)
                                                                                            aCCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATc                                                                                                  <  2:4880054/110‑1 (MQ=255)
                                                                                             ccGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAg                                                         >  2:5710258/1‑150 (MQ=255)
                                                                                              cGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAgg                                                        <  1:1704430/150‑1 (MQ=255)
                                                                                              cGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAgg                                                        <  1:1985726/150‑1 (MQ=255)
                                                                                              cGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAgg                                                        <  2:6713997/150‑1 (MQ=255)
                                                                                              cGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAgg                                                        <  1:1251144/150‑1 (MQ=255)
                                                                                              cGATTTGCCGCTGGCCAGCAGCTGGCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATAGGTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAgg                                                        <  2:6362602/150‑1 (MQ=255)
                                                                                                          ggCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTc                                            <  1:1381176/150‑1 (MQ=255)
                                                                                                          ggCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTc                                            <  2:1840540/150‑1 (MQ=255)
                                                                                                           gCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCg                                           <  1:1011414/150‑1 (MQ=255)
                                                                                                           gCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCg                                           <  1:6765135/150‑1 (MQ=255)
                                                                                                                    tGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGa                                  <  2:167934/150‑1 (MQ=255)
                                                                                                                      tCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGAATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACg                                >  2:6688011/1‑150 (MQ=255)
                                                                                                                        ggATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGtt                              >  1:3254066/1‑150 (MQ=255)
                                                                                                                        ggATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGtt                              <  1:1428254/150‑1 (MQ=255)
                                                                                                                               tgtAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTAc                       >  1:3506643/1‑150 (MQ=255)
                                                                                                                               tctggcagACGCCGCTTGCGGTGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTAc                       >  1:3943523/8‑150 (MQ=255)
                                                                                                                                   gcagACGCCGCTTGCGGTGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAgg                   >  2:5786377/4‑150 (MQ=255)
                                                                                                                                   gCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAgg                   <  1:4840402/150‑1 (MQ=255)
                                                                                                                                   gCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGAAATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAgg                   >  2:3909718/1‑150 (MQ=255)
                                                                                                                                    ccGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGa                  >  2:2770629/1‑150 (MQ=255)
                                                                                                                                     cGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGAt                 >  1:5475851/1‑150 (MQ=255)
                                                                                                                                     agACGCCGCTTGCGGTGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTAc                                               >  2:4305825/2‑120 (MQ=255)
                                                                                                                                     agACGCCGCTTGCGGTGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGAt                 >  1:6103674/2‑150 (MQ=255)
                                                                                                                                      gACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGGTCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGAtg                >  2:1784052/1‑150 (MQ=255)
                                                                                                                                      gACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGAtg                <  1:2340197/150‑1 (MQ=255)
                                                                                                                                      gACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGAtg                >  2:5175498/1‑150 (MQ=255)
                                                                                                                                      gACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGAtg                >  1:3213246/1‑150 (MQ=255)
                                                                                                                                      gACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGAtg                >  1:4987121/1‑150 (MQ=255)
                                                                                                                                       cgagcGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGAtgt               >  2:1810572/5‑150 (MQ=255)
                                                                                                                                           cGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGATGTGACa           >  2:2474944/1‑150 (MQ=255)
                                                                                                                                           cGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGATGTGACa           <  2:7209709/150‑1 (MQ=255)
                                                                                                                                            gCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGATGTGACAg          <  1:2081989/150‑1 (MQ=255)
                                                                                                                                                 cggcgTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGGTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGATGTGACAGTGTTg     <  2:2461476/150‑1 (MQ=255)
                                                                                                                                                   gcgTCCGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGATGTGACAGTGTTgca   <  2:2877383/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGGCGGTGGGATCGCTCTCCGAGCGCGCAGCCAACAGGTTGCGCAACTCGTCAGGGGGTTCTACGTCGGCATCCGGACGTTGCAGTACCAGGATGTGACAGTGTTGCAG  >  NC_002947/5477837‑5478134

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