Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,823,338 G→C 7.1% intergenic (+53/‑177) nadE → / → cho NAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,823,3380GC 7.1% 682.6 / 33.5 197intergenic (+53/‑177)nadE/choNAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair
Reads supporting (aligned to +/- strand):  ref base G (85/98);  new base C (9/5);  total (94/103)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

Only 100 of 208 total aligned reads displayed.

AAGCGCCAACGGCCGATCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCAC  >  NC_000913/1823080‑1823612
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                  
aaGCGCCAACGGCCGAGCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCg                                                                                                                                                                                                                                                                <  1:189594/279‑1 (MQ=255)
   cgcCAACGGCCGATCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcggc                                                                                                                                                                                                                                                                                >  2:10197/1‑260 (MQ=255)
    gcCGACGGACGATCTGGAAGATGATCGCCCGTCTCTGCCAGATGAAGTGGCACTCGGCCGGACCGATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGGACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAATAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAAGTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCg                                                                                                                                                                                                                                                                <  2:609976/275‑1 (MQ=255)
       aaCGGCCGATCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGa                                                                                                                                                                                                                                                                            <  4:10197/260‑1 (MQ=255)
         cGGCCGATCTGGAAGATGATCGCCCTTCGCTGCCAGATGAAGTGGCACGCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCta                                                                                                                                                                                                                                                           <  2:198506/275‑1 (MQ=255)
                          gatCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGANGTCCGCCAATTACCGTATTCGATGATTTATGGAAAAAGTAATACTTTGCACATATTGGATTGTGCGACAAAGAGTAATTTNGACACGCCGGATGCGGCGTGAACACNTTATTCGACATCTAAAACTATGCAAATTCAATa                                                                                                                                                                                                                                      >  3:157873/1‑279 (MQ=255)
                          gatCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGAAAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAAc                                                                                                                                                                                                                                                     <  1:266700/264‑1 (MQ=255)
                                          caaaTGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACGATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAAGAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGaa                                                                                                                                                                                                                      <  1:157873/276‑1 (MQ=255)
                                                  gTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCg                                                                                                                                                                                                                                                                <  2:44994/229‑1 (MQ=255)
                                                           ggCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAg                                                                                                                                                                                                          >  2:513566/1‑274 (MQ=255)
                                                                      tGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGc                                                                                                                                                                                                                                                <  1:231874/225‑1 (MQ=255)
                                                                                 acgacTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCg                                                                                                                                                                                                                                                                >  3:255229/1‑198 (MQ=255)
                                                                                                                    gTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGt                                                                                                                                                                                                                                                                              <  8:1254378/149‑1 (MQ=255)
                                                                                                                    gTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCAAACATAAACGACGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATCTTGGATTGTTCGAAAAAGCGTAATTTGTTCACGCCAGATGCGGCGCGAACGCCTGCTTCGCCCTATAAAACTATGCAAATTCAATACATTACCGGATTAGAATAGGCGTGACAGGCGTAGACCGTCAGAAGGTATAAACTTTGTCATCGCACCCCCTCTTTTTTAATCACt                                                                                                                                                 >  4:352659/1‑274 (MQ=255)
                                                                                                                     tCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTg                                                                                                                                                                                                                                                                             <  6:1585121/149‑1 (MQ=255)
                                                                                                                     tCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTg                                                                                                                                                                                                                                                                             <  8:13630/149‑1 (MQ=255)
                                                                                                                       gCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTGTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGaa                                                                                                                                                                                                                                                                           >  8:1261936/1‑149 (MQ=255)
                                                                                                                        ccAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAAc                                                                                                                                                                                                                                                                          <  7:510621/149‑1 (MQ=255)
                                                                                                                         cAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGc                                                                                                                                           >  4:54386/1‑275 (MQ=255)
                                                                                                                          aGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGc                                                                                                                                                                                                                                                                        >  7:846400/1‑149 (MQ=255)
                                                                                                                          aGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGc                                                                                                                                                                                                                                                                        >  8:434573/1‑149 (MQ=255)
                                                                                                                            aaCAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGAAcaaat                                                                                                                                                                                                                                                                      <  5:1153444/149‑5 (MQ=255)
                                                                                                                               aaTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGTGGCGTGAACGCCTTAt                                                                                                                                                                                                                                                                   >  8:1210356/1‑149 (MQ=255)
                                                                                                                                 tAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAAAGCCTTATTCg                                                                                                                                                                                                                                                                >  5:470082/1‑149 (MQ=255)
                                                                                                                                   gggaATTGGTAACTGAAAACCAAACATAAACGCCTTCCGCAAATTACCGTTTGCGATGATTTCTGGAAAAATTAATAATTTGCACATATTGGAGTGGGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAAGAAATTGCAGAGGTCGATTAGGCCTGACAGTGGTAGAAAGTAAGAAGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTaca                                                                                                                                   <  4:507032/271‑1 (MQ=255)
                                                                                                                                   gagaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAGGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGa                                                                                                                                                                                                                                                               >  8:1339910/1‑149 (MQ=255)
                                                                                                                                   gagaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGa                                                                                                                                                                                                                                                               <  6:1111032/149‑1 (MQ=255)
                                                                                                                                    agaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAc                                                                                                                                                                                                                                                              >  6:740820/1‑149 (MQ=255)
                                                                                                                                    agaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAc                                                                                                                                                                                                                                                              <  6:856432/149‑1 (MQ=255)
                                                                                                                                      aaCTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCt                                                                                                                                                                                                                                                                      >  2:732050/1‑139 (MQ=255)
                                                                                                                                      aaCTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCt                                                                                                                                                                                                                                                            <  6:1005428/149‑1 (MQ=255)
                                                                                                                                      aaCTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCTACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGga                                                                                                                              >  4:140889/1‑275 (MQ=255)
                                                                                                                                          ggTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTAtt                                                                                                                                                                                                                                                                  <  4:732045/139‑1 (MQ=255)
                                                                                                                                           gTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATaaa                                                                                                                                                                                                                                                       >  7:1449550/1‑149 (MQ=255)
                                                                                                                                           gTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTACTCGCTTACCGCCTTGTCCACTGgatagat                                                                                                                         >  2:561768/1‑275 (MQ=255)
                                                                                                                                                tGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTAt                                                                                                                                                                                                                                                  <  6:1577334/149‑1 (MQ=255)
                                                                                                                                                tGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTAt                                                                                                                                                                                                                                                  <  5:1526377/149‑1 (MQ=255)
                                                                                                                                                 gAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATg                                                                                                                                                                                                                                                 >  8:24658/1‑149 (MQ=255)
                                                                                                                                                 gAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGTCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCATCTTTTTTAATCGCTTACCTCATGTTACACTTGCTAGAGAACCTGCAtt                                                                                                               >  1:295908/1‑279 (MQ=255)
                                                                                                                                                  aaaaCCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGc                                                                                                                                                                                                                                                <  6:521087/149‑1 (MQ=255)
                                                                                                                                                         aaCATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCa                                                                                                                                                                                                                                         <  5:1269951/149‑1 (MQ=255)
                                                                                                                                                            aTAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATa                                                                                                                                                                                                                                      <  6:70849/149‑1 (MQ=255)
                                                                                                                                                                   ccGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCa                                                                                                                                                                                                                               <  7:1489061/149‑1 (MQ=255)
                                                                                                                                                                      tCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGa                                                                                                                                                                                                                            <  7:998336/149‑1 (MQ=255)
                                                                                                                                                                         gCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTc                                                                                                                                                                                                                         <  7:1302230/149‑1 (MQ=255)
                                                                                                                                                                             aTTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAAt                                                                                                                                                                                                                     >  5:1366623/1‑149 (MQ=255)
                                                                                                                                                                                aCCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAgg                                                                                                                                                                                                                  <  5:14502/149‑1 (MQ=255)
                                                                                                                                                                                   gTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCt                                                                                                                                                                                                               <  8:79969/149‑1 (MQ=255)
                                                                                                                                                                                   gTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCt                                                                                                                                                                                                               >  6:1189945/1‑149 (MQ=255)
                                                                                                                                                                                   gTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGa                                                                                                                                                                                                                                                                            >  7:1280425/1‑88 (MQ=255)
                                                                                                                                                                                     tttCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGa                                                                                                                                                                                                             <  5:1117909/149‑1 (MQ=255)
                                                                                                                                                                                      ttCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGAc                                                                                                                                                                                                            <  8:1279346/149‑1 (MQ=255)
                                                                                                                                                                                         gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAgg                                                                                                                                                                                                         >  5:745270/1‑149 (MQ=255)
                                                                                                                                                                                          atgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGa                                                                                                                                                                                                                                                               <  8:1130724/94‑1 (MQ=255)
                                                                                                                                                                                            gatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcg                                                                                                                                                                                                                                                                                  <  4:457010/73‑1 (MQ=255)
                                                                                                                                                                                            gatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGt                                                                                                                                                                                                      <  6:1122295/149‑1 (MQ=255)
                                                                                                                                                                                                tCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATa                                                                                                                                                                                                                    <  8:889227/131‑1 (MQ=255)
                                                                                                                                                                                                    gAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTc                                                                                                                                                                                              <  8:215344/149‑1 (MQ=255)
                                                                                                                                                                                                       aaaGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                                                                                                                                    >  2:348120/1‑76 (MQ=255)
                                                                                                                                                                                                           taataaTTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                                                                                                                                                       >  5:705575/1‑149 (MQ=255)
                                                                                                                                                                                                            aataatTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgtg                                                                                                                                                                                      <  7:928718/149‑1 (MQ=255)
                                                                                                                                                                                                                 tttGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAAcc                                                                                                                                                                                 <  8:569457/149‑1 (MQ=255)
                                                                                                                                                                                                                 tttGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAAcc                                                                                                                                                                                 <  7:767417/149‑1 (MQ=255)
                                                                                                                                                                                                                 tttGCACATATTGGATTGCGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTc                                                                                                                                                                                                                         >  6:1367347/1‑109 (MQ=255)
                                                                                                                                                                                                                    gCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCttt                                                                                                                                                                              <  6:756191/149‑1 (MQ=255)
                                                                                                                                                                                                                         tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAAt                                                                                                                                                                                                                                       <  8:27875/87‑1 (MQ=255)
                                                                                                                                                                                                                         tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCAt                                                                                                                                                                         >  7:560930/1‑149 (MQ=255)
                                                                                                                                                                                                                         tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCAt                                                                                                                                                                         >  6:369375/1‑149 (MQ=255)
                                                                                                                                                                                                                          atTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATc                                                                                                                                                                        <  8:1483206/149‑1 (MQ=255)
                                                                                                                                                                                                                          atTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATCACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGaa                                           >  2:270625/1‑274 (MQ=255)
                                                                                                                                                                                                                             ggATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAc                                                                                                                                                                     >  6:960803/1‑149 (MQ=255)
                                                                                                                                                                                                                              gATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAcc                                                                                                                                                                    <  5:924065/149‑1 (MQ=255)
                                                                                                                                                                                                                                ttGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCg                                                                                                                                                                  >  7:1121011/1‑149 (MQ=255)
                                                                                                                                                                                                                                   tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAGCGCCTTa                                                                                                                                                                                                                                                                    >  4:392169/1‑48 (MQ=255)
                                                                                                                                                                                                                                   tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGt                                                                                                                                                                                                                          >  4:574461/1‑90 (MQ=255)
                                                                                                                                                                                                                                   tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACg                                                                                                                                                                                                >  2:693488/1‑116 (MQ=255)
                                                                                                                                                                                                                                   tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGt                                                                                                                                                                                               >  4:237369/1‑117 (MQ=255)
                                                                                                                                                                                                                                   tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCct                                                                                                                                                               <  5:1189945/149‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAg                                                                                                                                                                                                    <  4:293464/108‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCttt                                                                                                                                                           >  8:98548/1‑149 (MQ=255)
                                                                                                                                                                                                                                       aaaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATGCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAg                                                                                                                                                                                            <  4:693483/116‑1 (MQ=255)
                                                                                                                                                                                                                                            gagTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAccc                                                                                                                                                                   <  7:1349108/136‑1 (MQ=255)
                                                                                                                                                                                                                                              gTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCa                                                                                                                                                                          >  7:1594422/1‑127 (MQ=255)
                                                                                                                                                                                                                                                 aTTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCt                                                                                                                                                 >  8:553941/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    tGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTcc                                                                                                                                              <  6:25642/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     gTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTccc                                                                                                                                             <  5:222445/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCCACCCGCCTCTTTTTTAATCGCATACCg                                                                                                                                            >  8:1089573/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      ttCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAAACTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCg                                                                                                                                            >  5:65362/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      ttCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATCGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCg                                                                                                                                            >  5:1201752/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        cACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGcc                                                                                                                                          <  5:1583246/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        cACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGcc                                                                                                                                          >  7:374397/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         aCGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCt                                                                                                                                         <  7:437662/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         aCGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCt                                                                                                                                         >  8:844722/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          cGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                                                                                                                                                       <  8:671144/102‑1 (MQ=255)
                                                                                                                                                                                                                                                           gCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                                                                                                                                                       >  8:1688154/1‑101 (MQ=255)
                                                                                                                                                                                                                                                           gCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACAGTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt                                                                                                                                                                                       <  7:1688154/101‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gcggcgTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGaa                                                                                                                                                                                                                      <  6:693086/63‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gcggcgTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTg                                                                                                                                <  7:263487/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gcggcgTGAACGCCTTATTCGACCTATAAAACTATGCAAATTAAATACATTGCAGGAGTCGAATAGGCTTGACAGGCGTAGCACGTCAGATGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTacccac  <  2:18536/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gcggcgTGAACGCCTAATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTacccac  <  2:312120/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                  
AAGCGCCAACGGCCGATCTGGAAGATGATCGCCCTTCTCTGCCAGATGAAGTGGCACTCGGCGTGACCTATGACAATATCGACGACTATCTGGAAGGGAAAAACGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCAC  >  NC_000913/1823080‑1823612

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