Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 340,058 T→C 5.2% intergenic (+315/‑107) yahG → / → yahI DUF1116 family protein/carbamate kinase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913340,0580TC 5.2% 293.9 / 3.9 94intergenic (+315/‑107)yahG/yahIDUF1116 family protein/carbamate kinase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (51/38);  new base C (3/2);  total (54/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

AGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCA  >  NC_000913/339913‑340204
                                                                                                                                                 |                                                                                                                                                  
aGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                                 <  4:310171/149‑1 (MQ=255)
      gcgcATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                            <  7:228231/148‑1 (MQ=255)
         cATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                         >  4:109413/1‑148 (MQ=255)
          aTCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                         >  2:139663/1‑147 (MQ=255)
             ctctaTTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTg                                                                                                                                    <  5:68946/145‑1 (MQ=255)
              ggCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                               <  2:99223/137‑1 (MQ=255)
              ggCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                               >  1:99223/1‑137 (MQ=255)
               gCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGgc                                                                                                                                  >  8:80672/1‑149 (MQ=255)
               gCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGgc                                                                                                                                  <  3:21109/149‑1 (MQ=255)
               gCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCGGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTCCGAGTGgc                                                                                                                                  >  6:360881/1‑149 (MQ=255)
               gCATTTGCCACCACTGCCGGATGCGGCGGGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGgc                                                                                                                                  <  8:111289/149‑1 (MQ=255)
                  tttGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTg                                                                                                                                    >  6:68946/1‑144 (MQ=255)
                    tGCCACCACTGCCTGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCga                                                                                                                              <  8:363482/148‑1 (MQ=255)
                        accacTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAtt                                                                                                                         >  2:377113/1‑149 (MQ=255)
                          cacTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTg                                                                                                                        >  6:181976/1‑148 (MQ=255)
                          cacTGCCGGATGCGGCGTGCACGCCTTATCCGGCCTACGAGAGGCGCGAGAATCTGTAGGCCTCAAAAGCGTAGCGCATCAGGCATTTGACACCATAGCCGGATGCGGCGTGAACGCCTTATAATGCCTACGAGTGGCGCCAGAATTg                                                                                                                        >  8:76709/1‑148 (MQ=14)
                           acTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGTTTTGCGCGAGATTCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTa                                                                                                                      <  1:192872/149‑1 (MQ=255)
                           acTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTa                                                                                                                      <  7:252881/149‑1 (MQ=255)
                               ccGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGcc                                                                                                                  >  5:57445/1‑149 (MQ=255)
                                  gATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTg                                                                                                                >  2:309382/1‑148 (MQ=255)
                                  gATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATATGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACAATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTg                                                                                                                >  8:209281/1‑148 (MQ=255)
                                            ggACGCCGTCTCCCGCCTACGAGCGGCGCGAGAATCTGTAGGCCTTATAAGCGTAGCGCATCAGGCATTTGCCACCATTGCCTGATTCGTCTTGAACGCCTCATCCGGCCTACGTGTGGCGAGATAATTGTAGGCCTGATAAGCGTAgc                                                                                                     <  2:123584/149‑1 (MQ=14)
                                             gACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAgcg                                                                                                    >  4:145623/1‑149 (MQ=255)
                                               cGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCa                                                                                                  >  5:50887/1‑149 (MQ=255)
                                                  cTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAt                                                                                                 >  6:33846/1‑147 (MQ=255)
                                                     aTCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                             <  1:309382/148‑1 (MQ=255)
                                                     aTCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGAATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                             >  8:349602/1‑148 (MQ=255)
                                                      tCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                           <  3:145623/149‑1 (MQ=255)
                                                        cGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGt                                                                                         <  8:72845/149‑1 (MQ=255)
                                                        cGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGt                                                                                         >  7:72845/1‑149 (MQ=255)
                                                        cGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCACATCAGGCAGt                                                                                         <  3:209239/149‑1 (MQ=255)
                                                         ggCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTAGTCACCATTGCCGGATGCGGCGTGACAGCCTTATCCGGCCTACGAGTGGCGCGCGAATTGTAGGCCTGATAAGCGTAGCGAATCATGCAGTc                                                                                        >  7:304008/1‑149 (MQ=14)
                                                          gCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATATGGGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCt                                                                                       <  7:76709/149‑1 (MQ=255)
                                                            cTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTgg                                                                                     >  2:370881/1‑149 (MQ=255)
                                                             tACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTa                                                                                                                      >  7:74228/1‑115 (MQ=255)
                                                             tACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCATTATCCGGCCTACGAGTGGCGCGAGAATTGTa                                                                                                                      <  8:74228/115‑1 (MQ=255)
                                                               cGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGt                                                                                  >  1:437996/1‑149 (MQ=255)
                                                                 actGGCGATAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCCCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGACCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTg                                                                                <  6:101313/147‑1 (MQ=255)
                                                                       gcgAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGc                                                                                    >  6:117121/1‑139 (MQ=255)
                                                                       gcgAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGc                                                                                    <  5:117121/139‑1 (MQ=255)
                                                                          agaATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAcc                                                                       >  5:34347/1‑149 (MQ=255)
                                                                                tGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAccc                                                                      <  6:120818/144‑1 (MQ=255)
                                                                                tGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAccc                                                                      >  5:120818/1‑144 (MQ=255)
                                                                                tGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCa                                                                 >  3:405541/1‑149 (MQ=255)
                                                                                tGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCa                                                                 >  3:333900/1‑149 (MQ=255)
                                                                                 gTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCAc                                                                >  8:68569/1‑149 (MQ=255)
                                                                                 gTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCAc                                                                >  5:109703/1‑149 (MQ=255)
                                                                                     gCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTa                                                                                                       >  2:368695/1‑106 (MQ=255)
                                                                                     gCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTa                                                                                                       <  1:368695/106‑1 (MQ=255)
                                                                                     gCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCctc                                                            >  7:211453/1‑149 (MQ=255)
                                                                                       cttaTAAGCGTAGGGCATCAGGCTTTTTTCACCATTGCGGGATGCGGCGTGAACGCCGTATCCGGCCTACGAGTGGCGCTAGATTTGTAGGCCTGATAAGCGTATCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTg                                                          <  5:229079/146‑1 (MQ=21)
                                                                                        tGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGAGCATCAGGCAGTCTGGCGTTGGTAATAACCCCATCACCCTCTg                                                          >  6:367200/1‑148 (MQ=255)
                                                                                         gATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCTAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTg                                                          >  3:308329/1‑147 (MQ=255)
                                                                                             aGCGTAGCTAATCAGTCATTTGTTACCATTGCCTGATGCGGCGTTAAGGCATTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCg                                                     <  6:10218/148‑1 (MQ=21)
                                                                                                 tAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACAt                                                <  6:57445/149‑1 (MQ=255)
                                                                                                     gcATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaac                                             >  5:192826/1‑148 (MQ=255)
                                                                                                      cATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGATAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaacaa                                           <  1:321505/149‑1 (MQ=255)
                                                                                                          aGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACAt                                        >  4:302503/1‑148 (MQ=255)
                                                                                                            gCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAccc                                                                      >  4:393658/1‑116 (MQ=255)
                                                                                                            gCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAccc                                                                      <  3:393658/116‑1 (MQ=255)
                                                                                                            gCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGa                                      >  5:253581/1‑148 (MQ=255)
                                                                                                                  gTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGc                                <  5:181976/148‑1 (MQ=255)
                                                                                                                  gTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCt                               >  8:317703/1‑149 (MQ=255)
                                                                                                                     aCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTtg                             <  7:209281/148‑1 (MQ=255)
                                                                                                                       cATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTgg                          >  4:177341/1‑149 (MQ=255)
                                                                                                                         ttGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGAt                                                                                                              <  4:332796/63‑1 (MQ=255)
                                                                                                                         ttGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGAt                                                                                                              >  3:332796/1‑63 (MQ=255)
                                                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGaaa                                    >  4:218391/3‑137 (MQ=25)
                                                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGaaa                                    <  3:218391/135‑1 (MQ=25)
                                                                                                                            ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCTAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGt                      <  4:308329/148‑1 (MQ=255)
                                                                                                                            ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                    >  2:119676/1‑118 (MQ=255)
                                                                                                                            ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                    <  1:119676/118‑1 (MQ=255)
                                                                                                                            ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGtt                     <  1:377113/149‑1 (MQ=255)
                                                                                                                            ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGtt                     >  5:34976/1‑149 (MQ=255)
                                                                                                                            ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGt                      <  3:109413/148‑1 (MQ=255)
                                                                                                                            ccGGATAAGGCGTTCACGCCGCATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAGCAGGCAGTCTGGCGTTGGTCATAACCCCCTCACCCTCTGTAGCGGACATACCACCCTGAAAGAGCTTGTGGTCGtt                     >  3:326731/1‑149 (MQ=25)
                                                                                                                             cGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAGCAGGCAGTCTGGCGGTGGTCATACCCCCCTCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTg                    >  3:18673/1‑149 (MQ=255)
                                                                                                                               gataaggcgTTCACGCCGCATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGc                   >  8:273479/6‑148 (MQ=37)
                                                                                                                                   cggcgTGAACGCCTTATCCGGCCTACGAGTGGCGCTAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATtg              >  8:172417/1‑149 (MQ=255)
                                                                                                                                   cggcgTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATtg              <  2:998/149‑1 (MQ=255)
                                                                                                                                    ggcgTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATtgg             >  5:213673/1‑149 (MQ=255)
                                                                                                                                     gcgTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaa                                              <  4:366487/115‑1 (MQ=255)
                                                                                                                                     gcgTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaa                                              >  3:366487/1‑115 (MQ=255)
                                                                                                                                     gagTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATtggt            <  2:437996/147‑1 (MQ=255)
                                                                                                                                      cgTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATtggtg           >  3:324236/1‑149 (MQ=255)
                                                                                                                                      agTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATtggtg           <  5:360881/148‑1 (MQ=255)
                                                                                                                                        tGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATtggtgg          >  3:329367/1‑148 (MQ=255)
                                                                                                                                          aaCGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCCCCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCaa       >  2:228078/1‑149 (MQ=255)
                                                                                                                                             gCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAg    >  5:83583/1‑149 (MQ=255)
                                                                                                                                             gCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAg    >  5:218498/1‑149 (MQ=255)
                                                                                                                                              ccTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGc   <  8:211453/149‑1 (MQ=255)
                                                                                                                                               cTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAaga                                  <  5:420349/117‑1 (MQ=255)
                                                                                                                                               cTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAaga                                  >  6:420349/1‑117 (MQ=255)
                                                                                                                                               cTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCa  <  3:254225/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCA  >  NC_000913/339913‑340204

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: