Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,976,758 C→T 5.3% 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 5.3% 369.5 / 7.2 112L54L (CTG→TTG) galRgalactose‑inducible d‑galactose regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (49/57);  new base T (3/3);  total (52/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

Only 100 of 115 total aligned reads displayed.

AAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTACCACAACGAACAA  >  NC_000913/2976610‑2976904
                                                                                                                                                    |                                                                                                                                                  
aaaGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCGCCCGAACGCCAACGCCCGTgcgc                                                                                                                                                    >  4:732112/1‑149 (MQ=255)
aaaGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTgcgc                                                                                                                                                    >  3:930988/1‑149 (MQ=255)
 aaGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCt                                                                                                                                                   >  2:1059412/1‑149 (MQ=255)
 aaGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCt                                                                                                                                                   >  4:336621/1‑149 (MQ=255)
  aGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTg                                                                                                                                                  >  2:946454/1‑149 (MQ=255)
    gATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGgc                                                                                                                                                >  3:118707/1‑149 (MQ=255)
        tAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAGCGCACGGGCGTTGgcgttc                                                                                                                                            <  3:336252/149‑4 (MQ=255)
         aGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAGCGCACGGGCGTTGgcgttcg                                                                                                                                           <  2:43343/149‑5 (MQ=255)
                  ggCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTg                                                                                                                                  <  4:324024/149‑1 (MQ=255)
                   gCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGa                                                                                                                                 >  1:182376/1‑149 (MQ=255)
                     aGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGaaa                                                                                                                               <  4:867053/149‑1 (MQ=255)
                         gTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAAcggt                                                                                                                           <  1:833086/149‑1 (MQ=255)
                              aGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACggtcggtc                                                                                                                      <  1:650106/149‑1 (MQ=255)
                               gTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCt                                                                                                                     >  2:803418/1‑149 (MQ=255)
                               gTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACACTGAAACGGTCGGTCt                                                                                                                     <  3:831955/148‑1 (MQ=255)
                               gTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACACTGAAACGGTCGGTCt                                                                                                                     >  4:831955/1‑148 (MQ=255)
                                 cGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTgg                                                                                                                   <  3:220432/149‑1 (MQ=255)
                                 cGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTgg                                                                                                                   <  1:507909/149‑1 (MQ=255)
                                  gCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAGCGCACGGGCGTTGgc                                                                                                                                                <  4:103440/119‑1 (MQ=255)
                                   ccACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCATACCCCTGAAACGGTCGGTCTGGTc                                                                                                                 <  4:503349/149‑1 (MQ=255)
                                   ccACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTc                                                                                                                 <  2:195614/149‑1 (MQ=255)
                                    cACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCg                                                                                                                >  1:661294/1‑149 (MQ=255)
                                     aCCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGt                                                                                                               >  1:353834/1‑149 (MQ=255)
                                       cGTTTCCCGGGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTg                                                                                                             <  1:452145/149‑1 (MQ=255)
                                          ttCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTg                                                                                                          >  3:949192/1‑149 (MQ=255)
                                           tCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGa                                                                                                         <  3:358610/149‑1 (MQ=255)
                                           tCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGa                                                                                                         <  4:1007764/149‑1 (MQ=255)
                                            cccGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGAt                                                                                                        <  1:803418/149‑1 (MQ=255)
                                                  tCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTcc                                                                                                  >  4:529353/1‑149 (MQ=255)
                                                  tCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTc                                                                                                   <  2:443687/148‑1 (MQ=255)
                                                     ttAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGaaa                                                                                                                               <  1:165913/117‑1 (MQ=255)
                                                     ttAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGAt                                                                                               >  3:278892/1‑149 (MQ=255)
                                                              cACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTc                                                                                      <  2:74294/149‑1 (MQ=255)
                                                              cACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTc                                                                                      >  4:838501/1‑149 (MQ=255)
                                                              cACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTc                                                                                      <  4:540261/149‑1 (MQ=255)
                                                               aCCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGt                                                                                           >  1:2102/1‑143 (MQ=255)
                                                               aCCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGt                                                                                           <  2:2102/143‑1 (MQ=255)
                                                                 ccAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCg                                                                                     <  4:330023/147‑1 (MQ=255)
                                                                 ccAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGt                                                                                   <  1:454457/149‑1 (MQ=255)
                                                                 ccAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGt                                                                                   >  4:362424/1‑149 (MQ=255)
                                                                   aaaGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGc                                                                                 >  4:614955/1‑149 (MQ=255)
                                                                   aaaGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGc                                                                                 >  1:688949/1‑149 (MQ=255)
                                                                   aaaGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGc                                                                                 >  2:35747/1‑149 (MQ=255)
                                                                   aaaGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGc                                                                                 <  4:817151/149‑1 (MQ=255)
                                                                       ccAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATg                                                                             <  4:118707/149‑1 (MQ=255)
                                                                        cAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATgg                                                                            >  1:144412/1‑149 (MQ=255)
                                                                         aGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATg                                                                             >  4:1003230/1‑147 (MQ=255)
                                                                         aGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATg                                                                             <  3:1003230/147‑1 (MQ=255)
                                                                         aGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCACAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGt                                                                           <  3:732112/149‑1 (MQ=255)
                                                                          gCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTg                                                                          >  1:588474/1‑149 (MQ=255)
                                                                          gCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTg                                                                          >  1:118559/1‑149 (MQ=255)
                                                                             aaGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGaaa                                                                       <  3:419698/149‑1 (MQ=255)
                                                                                cTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCg                                                                    <  3:336621/149‑1 (MQ=255)
                                                                                 ttCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCgg                                                                   >  2:774791/1‑149 (MQ=255)
                                                                                 ttCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCgg                                                                   <  4:325298/149‑1 (MQ=255)
                                                                                   cccGGCTTGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTc                                                                 >  2:569416/1‑149 (MQ=255)
                                                                                   cccGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTc                                                                 <  4:278892/149‑1 (MQ=255)
                                                                                    ccGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCg                                                                <  1:1059412/149‑1 (MQ=255)
                                                                                      ggctggctGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGaa                                                              <  3:772261/149‑1 (MQ=255)
                                                                                           gctgTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGt                                                         >  3:543769/1‑149 (MQ=255)
                                                                                           gctgTGCATAGTGCAATGGAGGCTCTTAGCTATCACCCGACCGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGt                                                         <  1:278388/149‑1 (MQ=255)
                                                                                                gCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCtt                                                    >  1:349539/1‑149 (MQ=255)
                                                                                                 cATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTa                                                   <  2:951334/149‑1 (MQ=255)
                                                                                                  aTAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCAc                                                                                                                                    >  1:752201/1‑67 (MQ=255)
                                                                                                  aTAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCAc                                                                                                                                    <  2:752201/67‑1 (MQ=255)
                                                                                                     gTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGt                                                                  <  1:415856/130‑1 (MQ=255)
                                                                                                     gTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGt                                                                  >  2:415856/1‑130 (MQ=255)
                                                                                                       gCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCacac                                             <  1:738517/149‑1 (MQ=255)
                                                                                                          aTGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCgg                                          >  2:636765/1‑149 (MQ=255)
                                                                                                          aTGGAGTCTCTTAGCTATCACCCGAACGACAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCgg                                          <  3:122679/149‑1 (MQ=255)
                                                                                                              aGTCTCTTAGCGATCACCCGATCGCCAACGCCCGTGCGCTTGCGCAGCAGTCCACTGATACCGTCTGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCACTGGTGAAAGCGGTCGAACAGGTGGCTCATCACACCGGTaaa                                      <  2:182376/149‑2 (MQ=255)
                                                                                                                  tctTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAAttttt                                  >  2:194493/1‑149 (MQ=255)
                                                                                                                  tctTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAAttttt                                  >  3:245726/1‑149 (MQ=255)
                                                                                                                    ttAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTa                                <  1:636765/149‑1 (MQ=255)
                                                                                                                      aGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTatt                              >  1:267574/1‑149 (MQ=255)
                                                                                                                      aGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTatt                              <  3:734058/149‑1 (MQ=255)
                                                                                                                        cTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTattga                            <  3:529353/149‑1 (MQ=255)
                                                                                                                        cTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTattga                            <  2:144412/149‑1 (MQ=255)
                                                                                                                         tATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGaa                                                              >  2:159781/1‑114 (MQ=255)
                                                                                                                         tATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGaa                                                              <  1:159781/114‑1 (MQ=255)
                                                                                                                         tATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTattgat                           <  1:520331/149‑1 (MQ=255)
                                                                                                                          aTCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTattgatt                          <  4:17728/149‑1 (MQ=255)
                                                                                                                          aTCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTattgatt                          <  2:180105/149‑1 (MQ=255)
                                                                                                                           tCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTattgattg                         <  4:543769/149‑1 (MQ=255)
                                                                                                                            cACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTgg                        <  2:136183/149‑1 (MQ=255)
                                                                                                                             aCCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTcc                                                                                                  <  4:106699/74‑1 (MQ=255)
                                                                                                                             aCCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTcc                                                                                                  >  3:106699/1‑74 (MQ=255)
                                                                                                                             aCCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGc                       >  1:853450/1‑149 (MQ=255)
                                                                                                                             aCCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGc                       >  4:399313/1‑149 (MQ=255)
                                                                                                                             aCCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGc                       >  3:890904/1‑149 (MQ=255)
                                                                                                                              cccGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCa                      >  1:552059/1‑149 (MQ=255)
                                                                                                                               tgctgcgccagcgcaCGGGCGTTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCaa                     >  4:59848/6‑149 (MQ=255)
                                                                                                                               ccGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCaa                     >  4:884351/1‑149 (MQ=255)
                                                                                                                                 gAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACg                   >  4:824314/1‑149 (MQ=255)
                                                                                                                                        aacgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTACca            >  1:342917/1‑149 (MQ=255)
                                                                                                                                         acgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCg                                           <  1:103921/117‑1 (MQ=255)
                                                                                                                                         acgccCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCg                                           >  2:103921/1‑117 (MQ=255)
                                                                                                                                             ccGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAACGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTACCACAACg       >  3:364054/1‑149 (MQ=255)
                                                                                                                                               gTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGGAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTACCACAACGaa     >  2:422627/1‑149 (MQ=255)
                                                                                                                                                  cgcTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTACCACAACGAACaa  <  2:661294/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAAGGATGTAGCCCGACTGGCAGGCGTTTCAGTCGCCACCGTTTCCCGCGTCATTAATAATTCACCCAAAGCCAGCGAAGCTTCCCGGCTGGCTGTGCATAGTGCAATGGAGTCTCTTAGCTATCACCCGAACGCCAACGCCCGTGCGCTGGCGCAGCAGACCACTGAAACGGTCGGTCTGGTCGTTGGTGATGTTTCCGATCCGTTTTTCGGTGCAATGGTGAAAGCGGTCGAACAGGTGGCTTATCACACCGGTAATTTTTTATTGATTGGCAACGGTTACCACAACGAACAA  >  NC_000913/2976610‑2976904

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