Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC41.2% 145.1 / 25.5 80T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (42/5);  new base C (0/33);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 1.57e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376356‑1376655
                                                                                                                                                         |                                                                                                                                                      
gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAacac                                                                                                                                                       >  1:480652/1‑151 (MQ=255)
     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCg                                                                                                                                                  <  2:426780/151‑1 (MQ=255)
        ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCg                                                                                                                                               <  2:130926/151‑1 (MQ=255)
        ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCg                                                                                                                                               <  2:1283858/151‑1 (MQ=255)
          gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACCCCCTCGCCGAt                                                                                                                                             <  2:449662/151‑1 (MQ=255)
             gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATGT‑GATCAACACCCTCGCCGATGcc                                                                                                                                          >  1:1268509/1‑151 (MQ=255)
                      gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                                 >  1:24953/1‑151 (MQ=255)
                       cgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                                 >  1:863797/1‑150 (MQ=255)
                        gtggCCTGGCCGGCGCCCGAAAAGTGGGGAGTGATCCCCTGCGCGTACGGAAAGGGGTCGCCGCGCGTACCGAACCGGTGG‑TTTATCTCGAC‑CCGTTAAC‑CCGCTATCGGGATCATAT‑GATAAACCCCCTCGCCGATGCCCGCCGTTACAt                                                                                                                               <  2:517709/148‑1 (MQ=255)
                        gTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACAt                                                                                                                               >  2:68990/1‑151 (MQ=255)
                             cTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                          >  2:360933/1‑151 (MQ=255)
                              tGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                          >  2:1382964/1‑150 (MQ=255)
                              tGGACGGCGACCGAAAAAGGGTGAGTGATTCCCTGCGCGTACTGAAACAGGTCCCCGCGCGTACCGGACCCGTGG‑TGTATCTCGAA‑CCGTTAAC‑CCGCTATCAGTATCATAT‑GATCAACCCCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                           <  1:982218/149‑1 (MQ=255)
                               ggcgggCGCCCGAAAAAGGGTGAGTGTTCCCCTGCGCGTACGGGAACGGGTCCCCGCGCGTCCCGGACCCGTGG‑TGTATCTCAAC‑CCGTTAAC‑CCGCTATCAGGATCATAT‑GACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGaa                                                                                                                        <  1:1142076/147‑1 (MQ=255)
                                gACGGCGACCGCAAAATGGTGTGTGATTCCCTGCGCGTACTGGAACAGGACGCCGCGCGTACCGGAACCGGGG‑TGTATCTCGAA‑CCGGTAAA‑TCGCTATCAGTATCATAT‑GATAAACACCCTCGCCGATGCCCGACGTTACATCGTCGa                                                                                                                         >  2:308795/1‑149 (MQ=255)
                                gACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGaa                                                                                                                        >  2:494915/1‑150 (MQ=255)
                                gACGGCGACCGCAAAATGGAGAGTGATTACCTGCGCGTACTGGCACAGGTCGCCGCGCGAACCGGAACCGGGG‑TGTATCTCGAA‑CCGTTAAA‑ACGCTATCAGGATCATAT‑GATAAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                         >  2:112031/1‑149 (MQ=255)
                                   ggCGCCCGAAAATTGGTGAGTTATCCCCTGCGTTTCCGGGACCAGGCCCCCGCGCGTCCGGGACCGGTGG‑TGTCCCCCAAC‑CCTTTAAC‑CCGCTACAAGGATCATAT‑AACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                     <  2:828329/150‑1 (MQ=255)
                                   ggCGCCCGAAAATGGGTGGGTGATTCCCGGCGCGTACGGAACCGGGTCGCCGCGCTGCCCGGACCCGTTG‑TTTTTCTCGAA‑CCTTAAAC‑CCTCTATAAGGAAAAAAT‑GATCAACCCCCGCGCCGATCCCCGCCGTTACATCGTCGaaaa                                                                                                                      <  2:187439/149‑1 (MQ=255)
                                   ggCGACCGAAAAATGGGGAGTTTTCCCCGGCGTTTCCGGGAACAGGTGGCCGCGCG‑CCCGGGACCCTGG‑TTTTTCCCCAACCCGTTAAC‑CCGCTATCAGGATAATAG‑GACAACCCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                    <  2:1423955/151‑1 (MQ=255)
                                       aCCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                                >  2:1331085/1‑151 (MQ=255)
                                        ccGAAAAATGGGGAGTGATTCCCTGCGCGTACGGAAACGGGTCGCCGCGCGTCCCGGACCGGTGG‑TTTATCTCGAC‑CCTTTAAC‑CCGCTATCAGGATCATAT‑GACAAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                                 <  1:1287078/149‑1 (MQ=255)
                                          gCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                              >  1:145781/1‑150 (MQ=255)
                                            aaaaTGGTGGTTGATTCCCTGCTCGTACTGGAACAGCCCGCCGCGGGTCCCGGAACGGTGT‑TGTTTCCCGAA‑CCTTTAAC‑CCGCT‑TCCGGGTTAAATGGATCACCCCCCTGGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                                 <  2:337433/145‑1 (MQ=255)
                                            aaaaTGGGGATTGTTTCCCTGCGCTTACGGAAACGGTTCGCCGCGCGTACC‑GAACCGGGGTTTTATCTC‑AA‑CCCTTAAACCCGCTATCAGTATCATAT‑AATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                              <  1:290522/148‑1 (MQ=255)
                                              aaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGACCCGTGG‑TGTATCTCGAA‑CCTTTAAC‑CCGCTATCAGAATCATAT‑GATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                         <  1:940381/151‑1 (MQ=255)
                                              aaTGGTGAGTGATTCCCTGCGCGTACTGAAACAGGTCGCCGCGCGTACCGGACCCGTGG‑TTTATCCCGAC‑CCTTTAAC‑CCGCTATCAGGATCATAT‑GACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                         <  1:174337/151‑1 (MQ=255)
                                                tGGTGAGTGATCCCCTGCGCGTACGGAAACAGGCCGCGGCGCGTCCCGGAACGGTGG‑TGTACCTCGAA‑CCGTAAAA‑CCGTTATCAGGATCATAT‑GATCAACCCCCCCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                        <  1:615940/150‑1 (MQ=255)
                                                  gtgggtgaTCCCCTGCGGGTACGGGAAGGGGTCGCCGCGCGCCCGGGACCCGTGT‑TGTATCTCAAC‑CCGTAAAC‑CCGATATAAGGATCATAT‑GAACACCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                         <  2:237574/147‑1 (MQ=255)
                                                   tgagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                     >  1:1263650/1‑150 (MQ=255)
                                                    gagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                     >  1:850453/1‑149 (MQ=255)
                                                      gttttccccTGCGCGTAGGGAAACGGGCCGCCGGGGGTCCGGGAACGGTGG‑TGTATCCCGAC‑CCGTAAAC‑CCGTTATCAGGATAATAT‑GACCACCCCCCTCGCCAATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGa                                                                                                  <  2:802644/144‑1 (MQ=255)
                                                                     tACTGGAACGGGTCGCCGCGCGTCCCGGACCCGTGT‑TGTATCTCGAC‑CCGTAAAC‑CCGTTATAAGGATCATAT‑GATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                   <  1:1160031/150‑1 (MQ=255)
                                                                          gAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                              >  1:935052/1‑150 (MQ=255)
                                                                          gAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                              >  1:1222456/1‑150 (MQ=255)
                                                                           accAGGTCGCCGCGCGTACGGGACCCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                            <  1:1351860/149‑1 (MQ=255)
                                                                           aaCAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                            >  1:1254645/1‑151 (MQ=255)
                                                                            acggtccgccgcGCGTCCCGGACCGGGGT‑TGTCTCTCGAA‑CCGTTACC‑CCGCTAGCGGAATCATAT‑AACAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                               <  1:300885/141‑1 (MQ=255)
                                                                             cggatcgcgcggcgTCCCGGAACTGGGT‑TTTTTCTCGAA‑CCGTAAAC‑CCGTTACTAGGAACATAT‑TTTAACCCCCCCCCCCTAAGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                           <  1:1403259/139‑1 (MQ=255)
                                                                              aGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                          >  1:478346/1‑150 (MQ=255)
                                                                                 tCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                                      >  1:304787/1‑151 (MQ=255)
                                                                                   gccgcGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                    >  2:1339728/1‑151 (MQ=255)
                                                                                   gccgcGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                    >  1:1203561/1‑151 (MQ=255)
                                                                                   gccgcGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                    >  1:937837/1‑151 (MQ=255)
                                                                                     cgcgcgTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                   >  1:892336/1‑150 (MQ=255)
                                                                                      gcgcgTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                   >  1:776641/1‑149 (MQ=255)
                                                                                      gcgcgTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                   >  1:327558/1‑149 (MQ=255)
                                                                                        gtttcccggcgcggtgg‑tgTACCTGAAA‑CTGTTAAC‑CCGTTACAAGGATAATAT‑GACCAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                          <  2:846952/127‑1 (MQ=255)
                                                                                        gcgTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  1:903665/1‑150 (MQ=255)
                                                                                        gcgTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  1:738842/1‑150 (MQ=255)
                                                                                        gcgTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  1:415894/1‑150 (MQ=255)
                                                                                         cgTCCCGAAACTTTGC‑TGTATCTCGCA‑CTTTTAAA‑CCGCTATAAGAATTTTAT‑GATAACCACCCCCCCCGCGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCCCATGAatat                                                                      <  1:535276/143‑1 (MQ=255)
                                                                                             cgggaccggggg‑ttTATCTC‑AA‑CCCTTAAACCCGCTATCAGGATCATAT‑GACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            <  1:68991/139‑1 (MQ=255)
                                                                                             ccGGACCGGTGG‑TTTATCTCGAA‑CCTTTAAC‑CCGCTATAAGGATCATAT‑GATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                <  1:557092/145‑1 (MQ=255)
                                                                                                   ccGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                    >  1:234227/1‑151 (MQ=255)
                                                                                                     gtgg‑tgTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                                  >  1:491387/1‑151 (MQ=255)
                                                                                                     gggg‑tgTATCTCGAC‑CCTTTAAC‑CCGTTATCAGGATAATAT‑GATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                       <  2:456495/144‑1 (MQ=255)
                                                                                                      tgg‑tgTATCTCGAC‑CCGTTAAC‑CCGCAATCAGGATCATAT‑GATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                     <  1:992838/147‑1 (MQ=255)
                                                                                                      gtgttttATCTCAAA‑CCTTTAAC‑CCGATATAAGGATCATAT‑GATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                    <  2:676658/143‑1 (MQ=255)
                                                                                                      ggtgttttcctcGAC‑CCTTTAAC‑CCGCTATCAGGATCATAT‑AATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                    <  2:1080214/140‑1 (MQ=255)
                                                                                                        gttgtACCTCGAC‑CCGTAAAC‑CCGCTACCAGAATCATAG‑AAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                                <  2:1396628/149‑1 (MQ=255)
                                                                                                          tgtATCCCGAC‑CCTTAAAC‑CCGCTATCAGGATCATAT‑GATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                               <  1:1540449/150‑1 (MQ=255)
                                                                                                              tcccGAC‑CCGTTAAC‑CCGTTATCAGGATCATAT‑GACCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                           <  1:196132/147‑1 (MQ=255)
                                                                                                                tcGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                                         >  1:1012869/1‑150 (MQ=255)
                                                                                                                  gAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                      >  1:104181/1‑151 (MQ=255)
                                                                                                                   aa‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                      >  1:248920/1‑150 (MQ=255)
                                                                                                                    acccGTAAAC‑CCGCTACAAGAATCAAAT‑GACAACCCCCCCCACCGTTCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCa                                        <  2:1180829/146‑1 (MQ=255)
                                                                                                                      ccGTTAAA‑CCGCTATCGGGTTCATAT‑GATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                                   <  2:525514/151‑1 (MQ=255)
                                                                                                                            aa‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                             >  1:225895/1‑151 (MQ=255)
                                                                                                                             acccGCTATCAGGATCATAT‑GACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                             <  1:436874/149‑1 (MQ=255)
                                                                                                                             a‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                             >  1:206721/1‑150 (MQ=255)
                                                                                                                                 gCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                         >  1:144564/1‑151 (MQ=255)
                                                                                                                                   tATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                       >  1:663762/1‑151 (MQ=255)
                                                                                                                                      cAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                    >  2:700196/1‑151 (MQ=255)
                                                                                                                                       aGGCTCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg                   >  2:626492/1‑151 (MQ=255)
                                                                                                                                           tCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt               >  1:260145/1‑151 (MQ=255)
                                                                                                                                              tat‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg            >  2:293291/1‑151 (MQ=255)
                                                                                                                                                  gATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat         >  2:183079/1‑151 (MQ=255)
                                                                                                                                                   aaaaaCCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat         <  1:367999/147‑1 (MQ=255)
                                                                                                                                                   aTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat         >  2:760138/1‑150 (MQ=255)
                                                                                                                                                      aaCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCata        <  1:699846/148‑1 (MQ=255)
                                                                                                                                                        cacCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg     <  2:983350/149‑1 (MQ=255)
                                                                                                                                                         acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  1:691057/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                      
GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAA‑CCGTTAAA‑CCGCTATCAGGATCATAT‑GATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376356‑1376655

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