Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,477,879 G→T 5.1% L225I (CTA→ATA)  nadD ← nicotinate‑nucleotide adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,477,8790GT 5.1% 806.2 / 22.9 237L225I (CTA→ATA) nadDnicotinate‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (133/92);  new base T (4/8);  total (137/100)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 240 total aligned reads displayed.

ACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGG  >  NC_002947/5477730‑5478028
                                                                                                                                                     |                                                                                                                                                     
  cTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGtt                                                                                                                                                     >  2:807026/1‑150 (MQ=255)
  cTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGtt                                                                                                                                                     >  1:548893/1‑150 (MQ=255)
         tCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGTTCTAGTTCGTgcgc                                                                                                                                              >  2:641761/1‑150 (MQ=255)
             gggCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGAACTATATCGTGCGcccaa                                                                                                                                          <  1:3935902/150‑4 (MQ=255)
             gggCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGAACTATATCGTGCGcccaa                                                                                                                                          <  2:1824982/150‑4 (MQ=255)
              ggCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGa                                                                                                                                         >  1:1527535/1‑150 (MQ=255)
                 cTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGatat                                                                                                                                      <  1:2327387/150‑1 (MQ=255)
                  ttGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGAtata                                                                                                                                     <  2:4618062/150‑1 (MQ=255)
                   tGACGTCTTCCAGCGCTGCCTTCGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAg                                                                                                                                    <  1:1845524/150‑1 (MQ=255)
                     aCGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGc                                                                                                                                  <  1:3015141/150‑1 (MQ=255)
                       gTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCa                                                                                                                                <  2:3966065/150‑1 (MQ=255)
                       gTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCa                                                                                                                                <  2:3319655/150‑1 (MQ=255)
                       gTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCa                                                                                                                                <  1:1991165/150‑1 (MQ=255)
                          ttCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTa                                                                                                                             >  2:5076197/1‑150 (MQ=255)
                             cAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCg                                                                                                                                          >  2:5205233/1‑134 (MQ=255)
                             cAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCg                                                                                                                                          <  1:5205233/134‑1 (MQ=255)
                              aGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGc                                                                                                                         >  2:1781983/1‑150 (MQ=255)
                               gcgcTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCg                                                                                                                        <  2:648432/150‑1 (MQ=255)
                                 gcTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTTAATTTTCTTCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCACACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTc                                                                                                                      <  2:5397098/150‑1 (MQ=255)
                                  cTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTcc                                                                                                                     <  2:1423322/150‑1 (MQ=255)
                                   tGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCg                                                                                                                    >  1:3883123/1‑150 (MQ=255)
                                    gCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCgg                                                                                                                   >  1:5547441/1‑150 (MQ=255)
                                     ccTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGc                                                                                                                  >  1:1748360/1‑150 (MQ=255)
                                       ttGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCAc                                                                                                                >  1:1365118/1‑150 (MQ=255)
                                       ttGGTCACTGCGACCAGTTCGTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCAc                                                                                                                <  2:4265622/150‑1 (MQ=255)
                                         ggTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCa                                                                                                              <  2:474969/150‑1 (MQ=255)
                                          gTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAg                                                                                                             <  2:2900358/150‑1 (MQ=255)
                                           tCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATCTAGGCCAGTACGGCGTCCGGCACCAgg                                                                                                            >  2:2998330/1‑150 (MQ=255)
                                               tGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAAcc                                                                                                        >  1:2720172/1‑150 (MQ=255)
                                                gCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCt                                                                                                       <  1:28974/150‑1 (MQ=255)
                                                gCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCt                                                                                                       <  1:1720145/150‑1 (MQ=255)
                                                 cGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCATTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGaa                                                                                                          <  2:4143126/146‑1 (MQ=255)
                                                 cGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCATTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGaa                                                                                                          >  1:4143126/1‑146 (MQ=255)
                                                    ccAGTTCTTCGCCGTACATGTTCGGCTGGGTCATATAACACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCAcc                                                                                                   <  2:5436286/150‑1 (MQ=255)
                                                    ccAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCAcc                                                                                                   >  1:3241177/1‑150 (MQ=255)
                                                     cAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCg                                                                                                  >  1:4582102/1‑150 (MQ=255)
                                                           ttcGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGc                                                                                            >  2:5620389/1‑150 (MQ=255)
                                                             cGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCg                                                                                          <  1:1358324/150‑1 (MQ=255)
                                                               ccGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCt                                                                                        <  2:4466922/150‑1 (MQ=255)
                                                                cGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTg                                                                                       >  1:253493/1‑150 (MQ=255)
                                                                  tAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGc                                                                                     >  2:759789/1‑150 (MQ=255)
                                                                  tAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGc                                                                                     >  1:3292019/1‑150 (MQ=255)
                                                                  tAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGc                                                                                     <  1:2072084/150‑1 (MQ=255)
                                                                  tAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGc                                                                                     <  1:3174574/150‑1 (MQ=255)
                                                                    aaTTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGTATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGACAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGACGCTGGCca                                                                                   >  2:4369098/1‑150 (MQ=255)
                                                                     aTTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCcag                                                                                  <  1:4293094/150‑1 (MQ=255)
                                                                            gCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGt                                                                           <  1:5566867/150‑1 (MQ=255)
                                                                            gCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGt                                                                           <  1:2741142/150‑1 (MQ=255)
                                                                                   cATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCCGAGGCACCGTCAGTTGGGCGCACGATCTAGTGCGTGCGCCTCGATATCGGCCAGTACGGCGTCCGGCACCAGGCACCGCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCt                                                                    >  2:3562432/1‑150 (MQ=255)
                                                                                     tataAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCtgt                                                                  <  1:1118181/150‑1 (MQ=255)
                                                                                     tataAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCtgt                                                                  <  2:1101016/150‑1 (MQ=255)
                                                                                     tataAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCtgt                                                                  >  1:5349742/1‑150 (MQ=255)
                                                                                      ataAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCtgtg                                                                 >  2:1037147/1‑150 (MQ=255)
                                                                                          aaCTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGc                                                             >  2:4945086/1‑150 (MQ=255)
                                                                                           aCTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGAtt                                                                                               <  2:3446494/115‑1 (MQ=255)
                                                                                             tCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGa                                                          <  2:1641464/150‑1 (MQ=255)
                                                                                               gTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGaca                                                        >  2:4867028/1‑150 (MQ=255)
                                                                                               gTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGaca                                                        >  1:4726891/1‑150 (MQ=255)
                                                                                               gTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGaca                                                        >  2:1505920/1‑150 (MQ=255)
                                                                                               gTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGaca                                                        >  2:1639892/1‑150 (MQ=255)
                                                                                               gTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGaca                                                        >  1:5131635/1‑150 (MQ=255)
                                                                                                   aaCTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGc                                                    >  1:4658438/1‑150 (MQ=255)
                                                                                                    aCTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCa                                                   >  1:1080182/1‑150 (MQ=255)
                                                                                                     cTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCaa                                                  <  1:4809531/150‑1 (MQ=255)
                                                                                                     cTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCaa                                                  >  1:504117/1‑150 (MQ=255)
                                                                                                      tCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAg                                                 >  2:1936631/1‑150 (MQ=255)
                                                                                                        gTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAgcg                                               >  1:2682180/1‑150 (MQ=255)
                                                                                                        gTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAgcg                                               >  1:3465496/1‑150 (MQ=255)
                                                                                                         ttGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAgcgg                                              >  2:2034448/1‑150 (MQ=255)
                                                                                                          tGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAgcggc                                             <  2:2887749/150‑1 (MQ=255)
                                                                                                          tGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAgcggc                                             <  1:1187247/150‑1 (MQ=255)
                                                                                                            gATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGt                                           <  1:3796438/150‑1 (MQ=255)
                                                                                                              tGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCt                                         <  1:4020339/150‑1 (MQ=255)
                                                                                                               gAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTg                                        >  1:3698189/1‑150 (MQ=255)
                                                                                                               gAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTg                                        >  2:3393915/1‑150 (MQ=255)
                                                                                                                aGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGc                                       <  1:4569236/150‑1 (MQ=255)
                                                                                                                 ggCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGcc                                      >  2:2083425/1‑150 (MQ=255)
                                                                                                                 ggCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGcc                                      >  2:5018777/1‑150 (MQ=255)
                                                                                                                  gcgcCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCa                                     >  1:3627337/1‑150 (MQ=255)
                                                                                                                  gcgcCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCa                                     >  1:4489800/1‑150 (MQ=255)
                                                                                                                    gcCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGa                                   >  1:334710/1‑150 (MQ=255)
                                                                                                                    gcCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGa                                   >  1:4904165/1‑150 (MQ=255)
                                                                                                                         agGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGacgaac                              <  2:4104200/150‑1 (MQ=255)
                                                                                                                            cACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGAcgaacgaa                           <  1:3222486/150‑1 (MQ=255)
                                                                                                                             aCCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGaaa                          <  1:2476245/150‑1 (MQ=255)
                                                                                                                             aCCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGaaa                          >  2:2696706/1‑150 (MQ=255)
                                                                                                                              ccGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAat                         >  1:3600629/1‑150 (MQ=255)
                                                                                                                                  cAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATAttt                     <  2:283112/150‑1 (MQ=255)
                                                                                                                                     cgaggCGCACGAACTATATCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGcc                                                                                           <  2:240717/74‑1 (MQ=37)
                                                                                                                                      tGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCccg                 <  1:3176839/150‑1 (MQ=255)
                                                                                                                                      tGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCccg                 <  1:4878272/150‑1 (MQ=255)
                                                                                                                                       gggCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCccgc                <  2:2499226/150‑1 (MQ=255)
                                                                                                                                       aggCGCACGAACTATATCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCccgc                >  2:4439750/2‑150 (MQ=255)
                                                                                                                                        ggCGCACGAACTATATCGTGCGCCTCGATATAGGCAAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCccgcc               >  1:640749/1‑150 (MQ=255)
                                                                                                                                         gcgcACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCccgccg              >  2:5309379/1‑150 (MQ=255)
                                                                                                                                             aCGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCgccgggcc          >  2:1324711/1‑150 (MQ=255)
                                                                                                                                               gATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCgccgggccgg        >  1:4369345/1‑150 (MQ=255)
                                                                                                                                               gATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCgccgggccgg        >  1:2195011/1‑150 (MQ=255)
                                                                                                                                                atatAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATGTCCCGCCGGGCCGGa       <  2:4821335/150‑1 (MQ=255)
                                                                                                                                                    aGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATg   >  2:277383/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGG  >  NC_002947/5477730‑5478028

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: