Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,823,338 G→C 6.3% 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 6.3% 256.6 / 10.1 95intergenic (+53/‑177)nadE/choNAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair
Reads supporting (aligned to +/- strand):  ref base G (34/53);  new base C (2/4);  total (36/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

CAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTG  >  NC_000913/1823190‑1823486
                                                                                                                                                    |                                                                                                                                                    
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:1312364/149‑1 (MQ=255)
                     gagaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGa                                                                                                                                 <  2:277777/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                                                                                                              <  1:1507573/149‑1 (MQ=255)
                                        ccGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAAt                                                                                                              <  2:827435/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)
                                                                                                                                                  aTGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGGATAGGCCTGACAGGCGTAGCACGACAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTacac    >  1:23316/1‑149 (MQ=255)
                                                                                                                                                    gcggcgTGAACGCCTTAATAGACCTATAAAACTATGCAAATTCAAAACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCGGCCTCTTTTTTAATCGCTTCACGCCTGTTACACTg  >  2:665218/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTG  >  NC_000913/1823190‑1823486

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: ‑