Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,823,335 G→C 16.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,3350GC16.0% 393.8 / 62.5 147intergenic (+50/‑180)nadE/choNAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair
Reads supporting (aligned to +/- strand):  ref base G (54/67);  new base C (13/10);  total (67/80)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 147 total aligned reads displayed.

TAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGT  >  NC_000913/1823237‑1823431
                                                                                                  |                                                                                                
tAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGa                                                                                                 <  2:2241880/100‑1 (MQ=255)
tAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGa                                                                                                 >  1:2969657/1‑100 (MQ=255)
 aaaCGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGAt                                                                                                <  2:2062446/100‑1 (MQ=255)
   aCGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgc                                                                                              <  2:1650688/100‑1 (MQ=255)
     gCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAACGAGAAATTTGTTCACGCAGGATgag                                                                                             >  1:3672044/1‑97 (MQ=255)
       cGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcggcg                                                                                          <  2:924003/100‑1 (MQ=255)
         tCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTg                                                                                        <  2:2334200/100‑1 (MQ=255)
         tCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTg                                                                                        >  1:3147688/1‑100 (MQ=255)
          ccGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGa                                                                                       >  2:1008923/1‑100 (MQ=255)
           cGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGaa                                                                                      >  2:1936642/1‑100 (MQ=255)
           cGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGaa                                                                                      >  1:2124224/1‑100 (MQ=255)
           cGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGaa                                                                                      <  1:962912/100‑1 (MQ=255)
            gCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAAc                                                                                     <  1:1028086/100‑1 (MQ=255)
            gCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAAc                                                                                     <  2:3736257/100‑1 (MQ=255)
            gCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGAAc                                                                                     <  2:3147705/100‑1 (MQ=255)
            gCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGAAc                                                                                     <  2:3600875/100‑1 (MQ=255)
             ccAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACg                                                                                    <  2:718722/100‑1 (MQ=255)
             ccAAGTACAGTTTTCGATAATTTCTGGAACAAGTAATCATTTACAAATATCGGACTGTGCGAACAAGTGTAATTTGTTCACGCCGCATCCGGCGTGAAc                                                                                     <  1:1824650/99‑1 (MQ=255)
               aaTTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGcc                                                                                  <  2:2203821/100‑1 (MQ=255)
               aaTTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGcc                                                                                  <  2:1581040/100‑1 (MQ=255)
                 ttACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCtt                                                                                <  1:2609057/100‑1 (MQ=255)
                  tACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTa                                                                               <  2:2858258/100‑1 (MQ=255)
                   aCCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTAt                                                                              >  2:1778147/1‑100 (MQ=255)
                      gTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCg                                                                           >  1:1441354/1‑100 (MQ=255)
                       ttttCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGa                                                                          <  1:2935083/100‑1 (MQ=255)
                        tttCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAc                                                                         <  1:3750753/100‑1 (MQ=255)
                        tttCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAc                                                                         <  2:2302076/100‑1 (MQ=255)
                        tttCGATGATTTATGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGACAAAGAGTAATTTGTTCACGCCGGCTGCGGCGTGAACGCCTTATTCGAc                                                                         <  1:3751589/100‑1 (MQ=255)
                         ttCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAcc                                                                        <  1:513555/100‑1 (MQ=255)
                            gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCtat                                                                     >  1:3457333/1‑100 (MQ=255)
                            gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCtat                                                                     >  2:1564654/1‑100 (MQ=255)
                            gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCtat                                                                     >  1:2654550/1‑100 (MQ=255)
                            gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCtat                                                                     >  1:2494572/1‑100 (MQ=255)
                             aagatTTCTGAAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCtata                                                                    <  1:1699061/98‑1 (MQ=255)
                                atTTCTGGAAAAAGTAATAATTTGCAGATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATaaaa                                                                 >  2:2313208/1‑100 (MQ=255)
                                 tttCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAAc                                                                >  1:1472678/1‑100 (MQ=255)
                                  ttCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACt                                                               <  1:3735215/100‑1 (MQ=255)
                                   tCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTa                                                              <  2:3669375/100‑1 (MQ=255)
                                      ggAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGc                                                           >  2:3757689/1‑100 (MQ=255)
                                        aaaaaGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCaa                                                         >  2:1272300/1‑100 (MQ=255)
                                        aaaaaGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCaa                                                         <  2:2281026/100‑1 (MQ=255)
                                         aaaaGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCaaa                                                        >  2:2115969/1‑100 (MQ=255)
                                            aGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTc                                                     <  2:2073923/100‑1 (MQ=255)
                                              taataaTTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCaa                                                   <  1:2946441/100‑1 (MQ=255)
                                                   atTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACAt                                              >  2:905229/1‑100 (MQ=255)
                                                      tGCACATATTAGATAAAGAAAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGc                                           <  1:1045321/100‑1 (MQ=25)
                                                       gCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCa                                          <  2:655688/100‑1 (MQ=255)
                                                        cacaTATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAg                                         >  2:199958/1‑100 (MQ=255)
                                                        cacaTATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAg                                         >  1:1719767/1‑100 (MQ=255)
                                                        cacaTATTGGATTGTGCGAAAAAGAGTAATTTGCTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAg                                         <  2:523982/100‑1 (MQ=255)
                                                         acaTATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAgg                                        <  1:642868/100‑1 (MQ=255)
                                                            tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGt                                     >  2:2744494/1‑100 (MQ=255)
                                                            tatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGt                                     >  2:822202/1‑100 (MQ=255)
                                                             atTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTc                                    >  2:3029747/1‑100 (MQ=255)
                                                             atTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTc                                    >  2:694125/1‑100 (MQ=255)
                                                              ttGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCg                                   <  2:1381095/100‑1 (MQ=255)
                                                               tGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGa                                  <  2:3136601/100‑1 (MQ=255)
                                                               tGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGa                                  <  1:3578669/100‑1 (MQ=255)
                                                               tGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGa                                  >  1:977887/1‑100 (MQ=255)
                                                                ggATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGaa                                 >  2:2931226/1‑100 (MQ=255)
                                                                   ttGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAg                              <  1:1027698/100‑1 (MQ=255)
                                                                   ttGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAg                              <  1:579635/100‑1 (MQ=255)
                                                                     gtgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGc                            <  1:1191856/100‑1 (MQ=255)
                                                                     gtgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGc                            >  2:1083259/1‑100 (MQ=255)
                                                                      tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAATTATGCAAATTCAATACATTGCAGGAGTCGAATAGGcc                           <  1:2038039/100‑1 (MQ=255)
                                                                      tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGcc                           <  2:822778/100‑1 (MQ=255)
                                                                      tgCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGcc                           <  2:1192333/100‑1 (MQ=255)
                                                                        cGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTg                         >  2:2144681/1‑100 (MQ=255)
                                                                        cGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTg                         <  2:1125612/100‑1 (MQ=255)
                                                                        aaaaaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGaa                                 <  1:3205394/90‑1 (MQ=255)
                                                                            aaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAg                     >  2:1051578/1‑100 (MQ=255)
                                                                               gagTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCg                  <  2:1672104/100‑1 (MQ=255)
                                                                               gagTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCg                  <  2:2356908/100‑1 (MQ=255)
                                                                               gaataaggcgTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCg                  >  2:54179/10‑100 (MQ=37)
                                                                                  tAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAg               <  1:1273056/100‑1 (MQ=255)
                                                                                   aaggcgTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGc              <  1:54177/95‑1 (MQ=37)
                                                                                   aaTTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCACATTCAATACATTACAGGAGTCGAATAGGCCGGACAGGCGTAGc              >  1:1158180/1‑100 (MQ=255)
                                                                                    aggcgTTCACGGCGCATCCGGCGTGAACGCCTTCTTCAACATCGACAAGTATGCAGATTCAAACCATTGCAGGAATCAAATCAGCCTGACAGGCGTaac              >  1:241173/5‑97 (MQ=12)
                                                                                    aTTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTGAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCa             >  2:617272/1‑100 (MQ=255)
                                                                                    aTTTGTTCACGCCGAATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGATTCGAATAGGCCTGACAGGCGTAGCa             <  1:1705123/100‑1 (MQ=255)
                                                                                      ttGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACg           <  1:3151804/100‑1 (MQ=255)
                                                                                       cgTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGt          >  2:2389613/2‑100 (MQ=255)
                                                                                       cgTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGt          >  2:2847434/2‑100 (MQ=255)
                                                                                       cgTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGt          >  2:2209120/2‑100 (MQ=255)
                                                                                       cgTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGt          >  2:1825441/2‑100 (MQ=255)
                                                                                       cgTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGt          >  2:6799/2‑100 (MQ=255)
                                                                                       cgTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGt          >  1:3164458/2‑100 (MQ=255)
                                                                                        gTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTc         >  2:2901424/1‑100 (MQ=255)
                                                                                        gTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTc         >  2:641954/1‑100 (MQ=255)
                                                                                        gTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTc         >  2:2830485/1‑100 (MQ=255)
                                                                                         ttCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCa        >  2:2783904/1‑100 (MQ=255)
                                                                                          tCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAg       >  2:3346409/1‑100 (MQ=255)
                                                                                           cACGCCGGATGTGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGa      >  1:3339104/1‑100 (MQ=255)
                                                                                           cACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGa      >  2:1416094/1‑100 (MQ=255)
                                                                                           cACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTACAGGAGTCGAATACGCCTGACAGGCGTAGCACGTCAGa      >  1:3056480/1‑100 (MQ=255)
                                                                                            aCGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGAc     <  1:1306567/100‑1 (MQ=255)
                                                                                            aCGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGAc     >  2:1528656/1‑100 (MQ=255)
                                                                                             cGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACg    >  2:3562339/1‑100 (MQ=255)
                                                                                             cGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACg    <  1:228952/100‑1 (MQ=255)
                                                                                               ccGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGgt  >  2:3290362/1‑100 (MQ=255)
                                                                                                  |                                                                                                
TAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGT  >  NC_000913/1823237‑1823431

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