Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,823,335 G→C 7.0% intergenic (+50/‑180) 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,3350GC 7.0% 251.5 / 10.1 96intergenic (+50/‑180)nadE/choNAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair
Reads supporting (aligned to +/- strand):  ref base G (33/54);  new base C (3/4);  total (36/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

CTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACA  >  NC_000913/1823188‑1823483
                                                                                                                                                   |                                                                                                                                                    
cTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTAACGCCGGa                                                                                                                                                     >  1:926693/1‑149 (MQ=255)
  cAACATTTCGCCAGAACAATAGAGAACTGTTATCTGAAAACCGAACATAAACGCCCTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATg                                                                                                                                                   <  2:1559512/149‑1 (MQ=255)
     cAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACCTACACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcgg                                                                                                                                                <  1:1416749/149‑1 (MQ=255)
        gTCGCCAGAACAATAGAGAACTGGTATCTGAACACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGAAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCAGGATGcgtcgt                                                                                                                                             >  2:1042940/1‑149 (MQ=255)
        gTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGAATGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGt                                                                                                                                             <  1:366185/149‑1 (MQ=255)
         gcGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTg                                                                                                                                            <  1:1421112/148‑1 (MQ=255)
          cGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAACAGTAATTTGTTCACGCCGCATCCGGCGTGa                                                                                                                                           <  2:356074/149‑1 (MQ=255)
           gCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGaa                                                                                                                                          >  1:593329/1‑149 (MQ=255)
           gCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGaa                                                                                                                                          <  2:570651/149‑1 (MQ=255)
            ccAGAACAATAGAGAACTGGTATCTGAAAACAGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAAc                                                                                                                                         <  2:372359/149‑1 (MQ=255)
             cAGAACAATAGAGAACTGGTTTCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACg                                                                                                                                        >  2:1098389/1‑149 (MQ=255)
             cAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTGGATGATATCTGGATAAAGTAATAATTTGCACATTTTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTTAAcc                                                                                                                                        >  1:1300677/1‑148 (MQ=255)
              aGAACACTAGAGAACTGGTATCTGAAACCCGAACATAAACCCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGAAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGc                                                                                                                                       <  2:764374/149‑1 (MQ=255)
              aGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAAAAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGc                                                                                                                                       <  2:810753/149‑1 (MQ=255)
                aaCAATTCAGAACTGGTATCTGAAACCCGAACATAACCGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCt                                                                                                                                     <  2:1341386/149‑1 (MQ=255)
                 aCAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGAAcaaatt                                                                                                                                    <  2:595428/149‑6 (MQ=255)
                   aaTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                  >  2:274532/1‑149 (MQ=255)
                     tAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACGGTTTTCGAGGAGTTCTGGAAAAAGTAATAATTTGCCCATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGACCGCATTATTc                                                                                                                                >  2:143649/1‑149 (MQ=255)
                       gagaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGa                                                                                                                              <  2:277777/149‑1 (MQ=255)
                       gagaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGa                                                                                                                              <  2:1312364/149‑1 (MQ=255)
                        agaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAc                                                                                                                             <  1:303154/149‑1 (MQ=255)
                        agaACTGGTATCTGAAAACCGAACATAAACGCCCTCCGCCAAATACCGTTTTCGATGAATTCTGGAAAAAGTACTAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAc                                                                                                                             <  1:955988/149‑1 (MQ=255)
                        agaACTGGTATCTGAAAACCGAACATAAACCCCGTCCGCCAAATACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAc                                                                                                                             <  1:381232/149‑1 (MQ=255)
                        agaACGGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGAGGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAc                                                                                                                             <  2:1290209/149‑1 (MQ=255)
                         gaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAcc                                                                                                                            <  1:255329/149‑1 (MQ=255)
                         caaCTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATATATTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAcc                                                                                                                            <  2:1239852/148‑1 (MQ=255)
                          aaCTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCt                                                                                                                           >  1:1003464/1‑149 (MQ=255)
                             tGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCtata                                                                                                                        >  2:284483/1‑149 (MQ=255)
                               gTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTTGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATaaa                                                                                                                      >  2:55026/1‑149 (MQ=255)
                                tATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATaaaa                                                                                                                     <  2:815327/149‑1 (MQ=255)
                                  tCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCCGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACt                                                                                                                   <  1:924362/149‑1 (MQ=255)
                                  tCTGAAAACCGAACATAAACGACGTCCGCCAATTACCGTTTTCGAAGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCTAAAAAGAGTAATTTGTTAACGCCGCATGCGGCGTGAACTCATTATTCGACCTATAAAACt                                                                                                                   >  2:152822/1‑149 (MQ=255)
                                   cTGAAAACCGAACATTAACGCCGTACGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTCGATTGTGCGAAAAAGAGTAATTTGTTCACGCAGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTa                                                                                                                  <  1:866179/149‑1 (MQ=255)
                                   cTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTa                                                                                                                  <  1:1505952/149‑1 (MQ=255)
                                   cTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTa                                                                                                                  <  1:57757/149‑1 (MQ=255)
                                    tGAAAACCGAACATAAACGCCGTCCGCCAATGACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTCTAAAACTAt                                                                                                                 >  1:1156168/1‑149 (MQ=255)
                                     gAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCTACCTATAAAACTATg                                                                                                                >  2:685496/1‑149 (MQ=255)
                                     gAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATg                                                                                                                <  1:485469/149‑1 (MQ=255)
                                        aaCCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCaa                                                                                                             >  1:655572/1‑149 (MQ=255)
                                          ccGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAAt                                                                                                           <  2:827435/149‑1 (MQ=255)
                                          ccGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAAt                                                                                                           <  1:1507573/149‑1 (MQ=255)
                                          ccGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGAAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAAt                                                                                                           <  2:197786/149‑1 (MQ=255)
                                           cGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAAtt                                                                                                          >  2:110047/1‑149 (MQ=255)
                                           cGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTATGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGAAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAAtt                                                                                                          <  1:891306/149‑1 (MQ=255)
                                            gAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTc                                                                                                         <  1:512683/149‑1 (MQ=255)
                                              aCATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCaa                                                                                                       <  1:97414/149‑1 (MQ=255)
                                               cATAACCGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGACGGATGCGGCGTGAACGCATTATTCGACCTATAAAACTATGCAAATTCAAt                                                                                                      >  2:290414/1‑149 (MQ=255)
                                               cATAAACGCCTTCCGCCAATTAGCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACTTATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAAt                                                                                                      <  1:1220858/149‑1 (MQ=255)
                                                  aaaCGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGAGTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTAAATACa                                                                                                   >  1:114735/1‑149 (MQ=255)
                                                        cGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAg                                                                                             >  2:1415145/1‑149 (MQ=255)
                                                             gCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTc                                                                                        <  1:1317695/149‑1 (MQ=255)
                                                             gCCAATTACCGTTTTCGAAGATTTCTGGAAAAAGTAATAATTTGCACAAATTGGATTGTCCGAAAAAGAGTAATTTGATCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTc                                                                                        <  1:873146/149‑1 (MQ=255)
                                                              ccAATTACCGTTTTCGATGATTTCTGGAAAAACTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGTGGCGTGAACGCCTTATTCGAGCTATAAAACTATGCAAATTAAAAACATTGCAGGAGTCg                                                                                       >  1:566370/1‑149 (MQ=255)
                                                                    aCCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAgg                                                                                 <  1:1355406/149‑1 (MQ=255)
                                                                    aCCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATATGTTCACGCCGGATGCGGCGTGAACGCCTTATTCTACCTATAAAACTATGCAAATTCAAGACATTGCAGGAGTCGAATAgg                                                                                 >  2:1526781/1‑149 (MQ=255)
                                                                     ccGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCTAATAGGc                                                                                >  2:1111301/1‑149 (MQ=255)
                                                                         tttCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGa                                                                            <  1:1312049/149‑1 (MQ=255)
                                                                             gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAgg                                                                        >  2:764984/1‑149 (MQ=255)
                                                                              atgatTTCTGGAAAAAGTAATAATTTGCACATATTTGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAACAGGCCTGACAGGc                                                                       <  1:1166609/149‑1 (MQ=255)
                                                                              atgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGc                                                                       >  2:1279485/1‑149 (MQ=255)
                                                                                tatTTCTGGATAAAGTAATAATTTGAACATATAGGATTGTTCGAAAAAGAGTAATTTGTTCACCCCGGATGCGGCGTGAACGCCTTAATCGACCTATAAAACTATGCAAATTCAAAACATTGCAGGAGTCGAATAGGCCTGACAGGCGt                                                                     <  1:866399/148‑1 (MQ=38)
                                                                                  tttCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATTCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAg                                                                   <  1:139261/149‑1 (MQ=255)
                                                                                   ttCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGc                                                                  >  2:965792/1‑149 (MQ=255)
                                                                                     atGGAAAAAGTAATAATTTGCAAATATTGGATTGTGAGAAAAAGAGTAATTTGTTCACGCAGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGACGAATAGGCCTGACAGGCGTAGCAc                                                                <  1:1076510/148‑1 (MQ=255)
                                                                                            aaGAAATAATTTGCACATATAGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGAc                                                         <  1:935316/149‑1 (MQ=255)
                                                                                             aGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACg                                                        <  1:802143/149‑1 (MQ=255)
                                                                                              gTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACgg                                                       >  1:719876/1‑149 (MQ=255)
                                                                                                      ttGCACATATTCGATTGTGCGAAAATGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCt                                               <  1:612338/149‑1 (MQ=255)
                                                                                                        gaacaTAATGGATTGTGCGAAAAAGAGTAATATCTTCACGCCGTATGCGGCGTGAACGCCTTATACGACCTAAAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCttt                                             <  1:870866/147‑1 (MQ=38)
                                                                                                          acaTATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGCACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGCCGGTGTAACCTTAGt                                           >  2:1232102/1‑149 (MQ=255)
                                                                                                           caTAATGGATACTGCGAACAAGAGTACTTTGTTCACGCCGGTCGCCGCGTGAACGCCTTATTCGACCTATTAAACTATGCAAATTTATTCCATTGCATGAGTCGAATAGGCCTGTCAGGCGTTGCACGTCAGACGGTGTAACCTTTggc                                          <  2:1083200/149‑3 (MQ=255)
                                                                                                           aatatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTc                                          <  1:385881/148‑1 (MQ=255)
                                                                                                             tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGGAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCAt                                        <  1:152822/149‑1 (MQ=255)
                                                                                                             tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCAt                                        <  1:284483/149‑1 (MQ=255)
                                                                                                              atTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATc                                       >  1:157362/1‑149 (MQ=255)
                                                                                                              atTGGATTGTGAGTAATTGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCATATTCGTCCTATAAAAGTTTGCAATTTCAATTCATTGCAGGTGTAGAATAGGCCTGACAGGCGTAGCTTGTCAGACGGTGTAACCTTTGTCATc                                       <  2:689864/149‑1 (MQ=25)
                                                                                                               ttGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCg                                      >  2:238273/1‑149 (MQ=255)
                                                                                                                tGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGa                                     <  2:506217/149‑1 (MQ=255)
                                                                                                                 ggATTGTGAGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAc                                    <  2:1074147/149‑1 (MQ=255)
                                                                                                                     tgtgCTAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGc                                <  1:1185110/149‑1 (MQ=255)
                                                                                                                       tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGTCCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCct                              >  1:877019/1‑149 (MQ=255)
                                                                                                                       tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGccg                              >  1:488246/1‑148 (MQ=255)
                                                                                                                       tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCct                              >  1:400190/1‑149 (MQ=255)
                                                                                                                            aaaaGAGTAATTTGTTCACGCCGTATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCtttt                         <  1:274532/149‑1 (MQ=255)
                                                                                                                            aaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGGGTAACCTTTGTCATCGACCCGCCTCtttt                         <  2:1129600/149‑1 (MQ=255)
                                                                                                                             aaaGAGTAATTTGTTCACGCCGGATTCGGCGTGAACGCCTTATTCGACCTATATAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCttttt                        <  1:371985/149‑1 (MQ=255)
                                                                                                                             aaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCttttt                        <  2:1237029/149‑1 (MQ=255)
                                                                                                                             aaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCttttt                        <  2:1124570/149‑1 (MQ=255)
                                                                                                                                     aTTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCt                >  1:1516164/1‑149 (MQ=255)
                                                                                                                                      tttGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCtt               <  2:1003464/149‑1 (MQ=255)
                                                                                                                                         gTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATAGATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGACAGACGGTGTAACCTTTGTCATCGAACCGCGGCTTTTTTAATCGCTTccc            >  1:134983/1‑149 (MQ=38)
                                                                                                                                          ttCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAACACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCg           >  1:838288/1‑149 (MQ=255)
                                                                                                                                          ttCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCg           >  1:1277776/1‑149 (MQ=255)
                                                                                                                                            cACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAACCTATGCAAATTCAATACAGTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATAGACCCGCCTCTTTTTTAATCGCGTCACGcc         >  1:857208/1‑149 (MQ=38)
                                                                                                                                             aCGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCt        <  1:1098389/149‑1 (MQ=255)
                                                                                                                                                ccGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGtt     <  1:1250139/149‑1 (MQ=255)
                                                                                                                                                   gATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTaca  <  2:157362/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACA  >  NC_000913/1823188‑1823483

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