Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,976,758 C→T 6.4% L54L (CTG→TTG)  galR → galactose‑inducible d‑galactose regulon transcriptional repressor; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,976,7580CT 6.4% 293.2 / 8.3 93L54L (CTG→TTG) galRgalactose‑inducible d‑galactose regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (44/43);  new base T (3/3);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

AAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTACCACAACGAACAAAA  >  NC_000913/2976611‑2976906
                                                                                                                                                   |                                                                                                                                                    
aaGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCt                                                                                                                                                     <  3:424604/149‑1 (MQ=255)
    aTGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAGCGCACGGGCGTTGgc                                                                                                                                                  <  2:87946/148‑1 (MQ=255)
     tGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAGCGCACGGGCGTTGgcg                                                                                                                                                 <  1:439708/148‑1 (MQ=255)
         gCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCgcagca                                                                                                                                            <  4:725529/149‑1 (MQ=255)
           ccGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGcagcag                                                                                                                                           <  2:883583/148‑1 (MQ=255)
           ccGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGa                                                                                                                                          >  3:309816/1‑149 (MQ=255)
                tGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACt                                                                                                                                     >  3:351551/1‑149 (MQ=255)
                        gTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAAcggt                                                                                                                             <  4:351551/149‑1 (MQ=255)
                            cAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAAcggtcgg                                                                                                                          <  2:354402/148‑1 (MQ=255)
                               tCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAAcgg                                                                                                                              <  1:157700/141‑1 (MQ=255)
                               tCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAAcgg                                                                                                                              >  2:157700/1‑141 (MQ=255)
                               tCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCt                                                                                                                       >  2:883829/1‑148 (MQ=255)
                                    aCCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCg                                                                                                                  <  2:15490/148‑1 (MQ=255)
                                    aCCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGt                                                                                                                 <  3:352613/149‑1 (MQ=255)
                                     ccGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGtt                                                                                                                <  3:166138/149‑1 (MQ=255)
                                      cGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGtt                                                                                                                >  2:819387/1‑148 (MQ=255)
                                      cGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGtt                                                                                                                >  2:204189/1‑148 (MQ=255)
                                       gTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTg                                                                                                               <  2:46775/148‑1 (MQ=255)
                                        tttCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGt                                                                                                             <  4:641330/149‑1 (MQ=255)
                                             cgcgTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGt                                                                                                        >  3:392470/1‑149 (MQ=255)
                                               cgTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGttt                                                                                                      <  3:390398/149‑1 (MQ=255)
                                                gTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTc                                                                                                     <  3:261825/149‑1 (MQ=255)
                                                   aTTAATAATTCACCGAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGGTATCACCCGAACGACAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGa                                                                                                  >  4:274007/1‑149 (MQ=255)
                                                   aTTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCg                                                                                                   <  2:717061/148‑1 (MQ=255)
                                                   aTTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGa                                                                                                  >  3:414224/1‑149 (MQ=255)
                                                    ttAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGaaa                                                                                                                                 >  3:317352/1‑117 (MQ=255)
                                                    ttAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGaaa                                                                                                                                 <  4:317352/117‑1 (MQ=255)
                                                    ttAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGAt                                                                                                 >  4:416376/1‑149 (MQ=255)
                                                    ttAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGAt                                                                                                 >  4:288451/1‑149 (MQ=255)
                                                     taataaTTCACCCAAAGCCAGCGAAGCTTCGCGGCTGGCTGTGCATTGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGAt                                                                                                 >  1:63135/1‑148 (MQ=255)
                                                        taatTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGa                                                                                                                                   >  1:37219/1‑111 (MQ=255)
                                                        taatTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGa                                                                                                                                   <  2:37219/111‑1 (MQ=255)
                                                                 cAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGt                                                                                     >  2:526127/1‑148 (MQ=255)
                                                                 cAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGt                                                                                     >  1:617288/1‑148 (MQ=255)
                                                                 cAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTg                                                                                    >  4:357250/1‑149 (MQ=255)
                                                                 cAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTg                                                                                    >  4:750768/1‑149 (MQ=255)
                                                                        aGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGt                                                                             <  3:750768/149‑1 (MQ=255)
                                                                               cTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGc                                                                       >  1:379697/1‑148 (MQ=255)
                                                                               cTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGc                                                                       <  2:379697/148‑1 (MQ=255)
                                                                                ttCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCgg                                                                     >  3:743070/1‑149 (MQ=255)
                                                                                ttCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCgg                                                                     <  3:415062/149‑1 (MQ=255)
                                                                                ttCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCgg                                                                     >  3:782607/1‑149 (MQ=255)
                                                                                 tCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTg                                                                                                               >  1:802656/1‑106 (MQ=255)
                                                                                 tCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTg                                                                                                               <  2:802656/106‑1 (MQ=255)
                                                                                          gctgTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGt                                                           <  4:309816/149‑1 (MQ=255)
                                                                                           ctgTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGt                                                           <  1:819387/148‑1 (MQ=255)
                                                                                            tgtgCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTg                                                          <  1:883829/148‑1 (MQ=255)
                                                                                             gtgCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAgg                                                            <  1:608105/145‑1 (MQ=255)
                                                                                             gtgCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAgg                                                            >  2:608105/1‑145 (MQ=255)
                                                                                                cATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCtt                                                      >  1:147494/1‑148 (MQ=255)
                                                                                                cATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCtt                                                      <  2:411059/148‑1 (MQ=255)
                                                                                                     tGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATcaca                                                >  4:117678/1‑149 (MQ=255)
                                                                                                     tGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATcaca                                                >  4:64792/1‑149 (MQ=255)
                                                                                                        aaTGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTAt                                                    <  3:200858/142‑1 (MQ=255)
                                                                                                        aaTGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTAt                                                    >  4:200858/1‑142 (MQ=255)
                                                                                                        aaTGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCg                                             >  3:466422/1‑149 (MQ=255)
                                                                                                            gAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTaa                                         <  4:889142/149‑1 (MQ=255)
                                                                                                                ctctTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGGGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAAttt                                      >  1:189877/1‑148 (MQ=255)
                                                                                                                 tctTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGt                                                                             <  3:449233/108‑1 (MQ=255)
                                                                                                                 tctTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGt                                                                             >  4:449233/1‑108 (MQ=255)
                                                                                                                 tctTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGGCTGGCCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAAttttt                                    >  4:642572/1‑149 (MQ=255)
                                                                                                                  ctTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAAtttttt                                   <  4:782607/149‑1 (MQ=255)
                                                                                                                  ctTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAAttttt                                    <  1:583049/148‑1 (MQ=255)
                                                                                                                   ttAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAAtttttt                                   <  1:526127/148‑1 (MQ=255)
                                                                                                                   ttAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTa                                  >  3:683247/1‑149 (MQ=255)
                                                                                                                    tAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTa                                  >  2:162626/1‑148 (MQ=255)
                                                                                                                    tAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTa                                  >  1:845180/1‑148 (MQ=255)
                                                                                                                    tAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTa                                  >  1:66969/1‑148 (MQ=255)
                                                                                                                      gCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGCCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTCTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTatt                                >  2:711507/1‑148 (MQ=255)
                                                                                                                        tATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTattgat                             <  3:818675/149‑1 (MQ=255)
                                                                                                                         aTCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTattgat                             <  1:711507/148‑1 (MQ=255)
                                                                                                                         aTCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTattgat                             >  1:251109/1‑148 (MQ=255)
                                                                                                                          tCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTCATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGGTTATCACACCGGTAATTTTTTattgattg                           <  4:795930/149‑1 (MQ=255)
                                                                                                                               cGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCaa                       >  1:691793/1‑148 (MQ=255)
                                                                                                                                  gcgccagcgcaCGGGCGTTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAAttt                                      <  1:65105/129‑1 (MQ=255)
                                                                                                                                  gcgccagcgcaCGGGCGTTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAAttt                                      >  2:65105/2‑130 (MQ=255)
                                                                                                                                  acgccaacgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGc                                                        >  3:666367/1‑112 (MQ=255)
                                                                                                                                  acgccaacgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGc                                                        <  4:666367/112‑1 (MQ=255)
                                                                                                                                  acgccaacgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGt                   <  3:360409/149‑1 (MQ=255)
                                                                                                                                   cgccaacgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGtt                  <  3:117678/149‑1 (MQ=255)
                                                                                                                                    gccaacgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTa                 >  4:465488/1‑149 (MQ=255)
                                                                                                                                     ccagcgcaCGGGCGTTGGCGCAGCAGACTACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTAc                >  3:776864/1‑149 (MQ=255)
                                                                                                                                     ccagcgcaCGGGCGTTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTa                 >  2:438679/5‑148 (MQ=255)
                                                                                                                                     ccaacgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTa                 <  2:66969/148‑1 (MQ=255)
                                                                                                                                     ccaacgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTAc                <  4:683247/149‑1 (MQ=255)
                                                                                                                                      caacgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTAcc               >  3:358024/1‑149 (MQ=255)
                                                                                                                                       aacgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTa                                                     >  3:246804/1‑110 (MQ=255)
                                                                                                                                       aacgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTa                                                     <  4:246804/110‑1 (MQ=255)
                                                                                                                                       aacgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTAc                <  4:667813/147‑1 (MQ=255)
                                                                                                                                       aacgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTAc                >  3:667813/1‑147 (MQ=255)
                                                                                                                                          gccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTCTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTACCACaa           >  4:63257/1‑149 (MQ=255)
                                                                                                                                           cccGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTACCACAAc          <  4:417757/149‑1 (MQ=255)
                                                                                                                                           cccGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTACCACAAc          <  3:114330/149‑1 (MQ=255)
                                                                                                                                                   cTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTACCACAACGAACaaaa  >  3:676308/1‑149 (MQ=255)
                                                                                                                                                   cTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTACCACAACGAACaaa   >  1:487914/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTACCACAACGAACAAAA  >  NC_000913/2976611‑2976906

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: