Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,807,093 A→T 8.9% intergenic (+19/‑39) proW → / → proX glycine betaine/proline ABC transporter permease/glycine betaine/proline ABC transporter periplasmic binding protein; cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,807,0930AT 8.9% 267.7 / 14.7 79intergenic (+19/‑39)proW/proXglycine betaine/proline ABC transporter permease/glycine betaine/proline ABC transporter periplasmic binding protein; cold shock protein
Reads supporting (aligned to +/- strand):  ref base A (37/35);  new base T (4/3);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

TTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCA  >  NC_000913/2806947‑2807240
                                                                                                                                                  |                                                                                                                                                   
ttGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGAc                                                                                                                                                    >  1:698112/1‑148 (MQ=255)
    gATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCg                                                                                                                                                <  2:698112/148‑1 (MQ=255)
          cGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGc                                                                                                                                             <  2:255006/145‑1 (MQ=255)
          cGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGc                                                                                                                                             >  1:255006/1‑145 (MQ=255)
          cGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCa                                                                                                                                            >  1:659458/1‑146 (MQ=255)
          cGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCa                                                                                                                                            <  2:659458/146‑1 (MQ=255)
                               gACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGa                                                                                                                     >  4:213621/1‑148 (MQ=255)
                                  gCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATa                                                                                                                  <  3:213621/148‑1 (MQ=255)
                                                     tCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCa                                                                                                                                            >  4:145432/1‑103 (MQ=255)
                                                     tCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCa                                                                                                                                            <  3:145432/103‑1 (MQ=255)
                                                      cACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACt                                                                                              >  1:745139/1‑148 (MQ=255)
                                                        cGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCa                                                                                                                                            >  4:394586/1‑100 (MQ=255)
                                                        cGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCa                                                                                                                                            <  3:394586/100‑1 (MQ=255)
                                                        cGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACttt                                                                                            >  4:260138/1‑148 (MQ=255)
                                                           aGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTg                                                                                         <  3:260138/148‑1 (MQ=255)
                                                                                tACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCg                                                                                                                                                <  1:549901/72‑1 (MQ=255)
                                                                                tACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCg                                                                                                                                                >  2:549901/1‑72 (MQ=255)
                                                                                tACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTTCGACAGCGTTTGCCACGCTTa                                                                    >  4:1075412/1‑148 (MQ=255)
                                                                                 acacCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGtt                                                                               <  2:263772/136‑1 (MQ=255)
                                                                                 acacCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGtt                                                                               >  1:263772/1‑136 (MQ=255)
                                                                                   accacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCg                                                                                       >  4:727314/1‑126 (MQ=255)
                                                                                   accacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCg                                                                                       <  3:727314/126‑1 (MQ=255)
                                                                                    ccacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtctc                                                                <  2:745139/148‑1 (MQ=255)
                                                                                      acTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCg                                                                                                          <  4:956227/104‑1 (MQ=255)
                                                                                      acTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCg                                                                                                          >  3:956227/1‑104 (MQ=255)
                                                                                      acTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGtt                                                                               >  4:974931/1‑131 (MQ=255)
                                                                                      acTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGtt                                                                               <  3:974931/131‑1 (MQ=255)
                                                                                       cTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTac                                                             <  3:1075412/148‑1 (MQ=255)
                                                                                         ggCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGcc                                                                                                                                                 >  4:36060/1‑62 (MQ=255)
                                                                                         ggCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGcc                                                                                                                                                 <  3:36060/62‑1 (MQ=255)
                                                                                            ccTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACaaa                                                        >  4:108561/1‑148 (MQ=255)
                                                                                              tGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAg                                                                                                    <  4:242080/102‑1 (MQ=255)
                                                                                              tGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAg                                                                                                    >  3:242080/1‑102 (MQ=255)
                                                                                               gTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAt                                                                   <  4:518826/134‑1 (MQ=255)
                                                                                               gTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAt                                                                   >  3:518826/1‑134 (MQ=255)
                                                                                                 tGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGttt                                                                              >  4:394046/1‑121 (MQ=255)
                                                                                                 tGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGttt                                                                              <  3:394046/121‑1 (MQ=255)
                                                                                                   gTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACa                                                                                                               <  1:390493/86‑1 (MQ=255)
                                                                                                   gTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACa                                                                                                               >  2:390493/1‑86 (MQ=255)
                                                                                                        ctgACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACtttt                                                   >  4:700164/1‑141 (MQ=255)
                                                                                                        ctgACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACtttt                                                   <  3:700164/141‑1 (MQ=255)
                                                                                                           aCCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACaaa                                                        >  1:584718/1‑133 (MQ=255)
                                                                                                           aCCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACaaa                                                        <  2:584718/133‑1 (MQ=255)
                                                                                                                      cattAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCa                                                                          >  1:1112524/1‑104 (MQ=255)
                                                                                                                      cattAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCa                                                                          <  2:1112524/104‑1 (MQ=255)
                                                                                                                       attAAGTAACTCTGCACTTGCCCGGTGACGCCTGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCaaa                              <  3:255888/147‑1 (MQ=255)
                                                                                                                       attAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACttt                                                                                            >  4:513184/1‑85 (MQ=255)
                                                                                                                       attAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACttt                                                                                            <  3:513184/85‑1 (MQ=255)
                                                                                                                       attAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCaaa                              >  4:255888/1‑147 (MQ=255)
                                                                                                                        ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATa                                                                                                                  >  2:677151/1‑62 (MQ=255)
                                                                                                                        ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATa                                                                                                                  <  1:677151/62‑1 (MQ=255)
                                                                                                                        ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtct                                                                 <  4:679215/111‑1 (MQ=255)
                                                                                                                        ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtct                                                                 >  3:679215/1‑111 (MQ=255)
                                                                                                                        ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCaa                               <  3:384139/145‑1 (MQ=255)
                                                                                                                        ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCaa                               >  4:384139/1‑145 (MQ=255)
                                                                                                                             tAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACGCTGCCAAAAAAAGGAATAACAATGCGACATAg                                                                                                    >  3:333805/1‑71 (MQ=255)
                                                                                                                             tAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCaaaaaaa                                                                                                                        <  2:1012528/51‑1 (MQ=255)
                                                                                                                             tAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCaaaaaaa                                                                                                                        >  1:1012528/1‑51 (MQ=255)
                                                                                                                             tAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAg                                                                                                    <  4:333805/71‑1 (MQ=255)
                                                                                                                             tAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGc                           <  3:1087680/144‑1 (MQ=255)
                                                                                                                             tAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGc                           >  4:1087680/1‑144 (MQ=255)
                                                                                                                              aaCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTAc                      <  3:108561/148‑1 (MQ=255)
                                                                                                                                  ctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGt                                                                                <  2:443848/86‑1 (MQ=255)
                                                                                                                                  ctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGt                                                                                >  1:443848/1‑86 (MQ=255)
                                                                                                                                  ctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGtt                  >  4:731839/1‑148 (MQ=255)
                                                                                                                                        atGCCCGGCGTCACCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAg                                                                                                    >  3:498803/2‑60 (MQ=37)
                                                                                                                                         tGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGt                                                                                                 <  2:462663/62‑1 (MQ=255)
                                                                                                                                         tGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGt                                                                                                 >  1:462663/1‑62 (MQ=255)
                                                                                                                                         tGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTg                                                                             <  2:25366/82‑1 (MQ=255)
                                                                                                                                         tGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTg                                                                             >  1:25366/1‑82 (MQ=255)
                                                                                                                                         tGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAg           >  2:980015/1‑148 (MQ=255)
                                                                                                                                             cGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATACCAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtctat                                                               <  4:761268/92‑3 (MQ=255)
                                                                                                                                             cGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACaa                                                         >  4:68601/1‑98 (MQ=255)
                                                                                                                                             cGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACaa                                                         <  3:68601/98‑1 (MQ=255)
                                                                                                                                             cGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTg                    <  4:939363/135‑1 (MQ=255)
                                                                                                                                             cGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTg                    >  3:939363/1‑135 (MQ=255)
                                                                                                                                             cGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGAGTACTTTTTGCGACAGCGTTTGCCACGCTTATctct                                                               >  3:761268/1‑92 (MQ=255)
                                                                                                                                               gTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACtttt                                                   <  3:835687/102‑1 (MQ=255)
                                                                                                                                               gTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACtttt                                                   >  4:835687/1‑102 (MQ=255)
                                                                                                                                                  aCGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCa  <  3:731839/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCA  >  NC_000913/2806947‑2807240

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