Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,477,880 T→A 5.4% E224D (GAA→GAT nadD ← nicotinate‑nucleotide adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,477,8800TA 5.4% 794.9 / 22.4 234E224D (GAA→GATnadDnicotinate‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (128/92);  new base A (5/8);  total (133/101)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

CTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGG  >  NC_002947/5477732‑5478028
                                                                                                                                                    |                                                                                                                                                    
cTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGtt                                                                                                                                                     >  2:807026/1‑150 (MQ=255)
cTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGtt                                                                                                                                                     >  1:548893/1‑150 (MQ=255)
       tCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGTTCTAGTTCGTgcgc                                                                                                                                              >  2:641761/1‑150 (MQ=255)
         cTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGAACTATATCGTGCGccc                                                                                                                                            <  1:1869516/150‑2 (MQ=255)
           gggCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGAACTATATCGTGCGcccaa                                                                                                                                          <  1:3935902/150‑4 (MQ=255)
           gggCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGAACTATATCGTGCGcccaa                                                                                                                                          <  2:1824982/150‑4 (MQ=255)
               cTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGatat                                                                                                                                      <  1:2327387/150‑1 (MQ=255)
                ttGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGAtata                                                                                                                                     <  2:4618062/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)
                        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)
                                     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:1720263/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                                                                                     <  1:3174574/150‑1 (MQ=255)
                                                                tAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGc                                                                                     >  2:759789/1‑150 (MQ=255)
                                                                tAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGc                                                                                     >  1:3292019/1‑150 (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)
                                                                           cTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTc                                                                          >  1:1614936/1‑150 (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                                                        >  1:4726891/1‑150 (MQ=255)
                                                                                             gTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGaca                                                        >  1:5131635/1‑150 (MQ=255)
                                                                                             gTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGaca                                                        >  2:1505920/1‑150 (MQ=255)
                                                                                             gTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGaca                                                        >  2:1639892/1‑150 (MQ=255)
                                                                                             gTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGaca                                                        >  2:4867028/1‑150 (MQ=255)
                                                                                              ttCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGacac                                                       >  1:1822431/1‑150 (MQ=255)
                                                                                                 aaCTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGc                                                    >  1:4658438/1‑150 (MQ=255)
                                                                                                  aCTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCa                                                   >  1:1080182/1‑150 (MQ=255)
                                                                                                   cTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCaa                                                  >  1:504117/1‑150 (MQ=255)
                                                                                                   cTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCaa                                                  <  1:4809531/150‑1 (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:5018777/1‑150 (MQ=255)
                                                                                                               ggCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGcc                                      >  2:2083425/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)
                                                                                                                                    tGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCccg                 <  1:2046512/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:2195011/1‑150 (MQ=255)
                                                                                                                                             gATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCgccgggccgg        >  1:4369345/1‑150 (MQ=255)
                                                                                                                                              atatAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATGTCCCGCCGGGCCGGa       <  2:4821335/150‑1 (MQ=255)
                                                                                                                                                  aGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATg   >  2:277383/1‑150 (MQ=255)
                                                                                                                                                   tttCGTTGGCCTCGTTTTAGGGCTGTTCGGCGTCCGGCACCATGAAACTCAACGTTTTGCCGCTGGCCAGCAGCTGTCGGATATGTGTAGCCGACACCGCAAGCGGCTTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATgg  <  1:2190397/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGACATGG  >  NC_002947/5477732‑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: