Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,828,537 A→T 5.1% intergenic (+21/‑84) srlD → / → srlM sorbitol‑6‑phosphate dehydrogenase/sorbitol=responsive srl operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,828,5370AT 5.1% 544.4 / 13.5 157intergenic (+21/‑84)srlD/srlMsorbitol‑6‑phosphate dehydrogenase/sorbitol=responsive srl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (84/65);  new base T (3/5);  total (87/70)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 162 total aligned reads displayed.

AAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTG  >  NC_000913/2828389‑2828685
                                                                                                                                                    |                                                                                                                                                     
aaGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCGATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTaa                                                                                                                                                       <  2:448945/149‑1 (MQ=255)
aaGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTaa                                                                                                                                                       <  4:1262389/149‑1 (MQ=255)
aaGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTaa                                                                                                                                                       <  1:1486785/149‑1 (MQ=255)
 aGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCAGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑g                                                                                                                                                     <  4:618892/149‑1 (MQ=255)
     ccGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCCGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGAtc                                                                                                                                                 >  4:485594/1‑149 (MQ=255)
       gCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGAtctc                                                                                                                                               >  3:326407/1‑149 (MQ=255)
        cTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCTTAAT‑GGATCTcc                                                                                                                                              <  3:1268578/149‑1 (MQ=255)
         tCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCg                                                                                                                                             >  1:1397201/1‑149 (MQ=255)
          cAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTTGATCACTGTCACCGGCGGTCAGGTGATGTTCTGAACAACAGCGGAGATCCATTAA‑GGATCTCCGt                                                                                                                                            <  4:397871/149‑1 (MQ=255)
           aaaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCTTAAT‑GCATCTCCgct                                                                                                                                           <  4:1195138/149‑3 (MQ=255)
           aaaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTg                                                                                                                                           <  4:651015/149‑1 (MQ=255)
           aaaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTg                                                                                                                                           <  2:933626/149‑1 (MQ=255)
            aaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTga                                                                                                                                          <  1:1017993/149‑1 (MQ=255)
            aaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTga                                                                                                                                          <  4:987518/149‑1 (MQ=255)
            aaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTga                                                                                                                                          >  3:1364713/1‑149 (MQ=255)
              cgcgGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTgaga                                                                                                                                        <  3:215103/149‑1 (MQ=255)
                cgGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACt                                                                                                                                      <  2:364784/149‑1 (MQ=255)
                          tATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAg                                                                                                                                  >  1:828713/1‑143 (MQ=255)
                            tCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGAc                                                                                                                                       <  2:1374983/136‑1 (MQ=255)
                              aaGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGAt                                                                                                                        <  3:1179681/149‑1 (MQ=255)
                              aaGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGAt                                                                                                                        <  1:1305312/149‑1 (MQ=255)
                                  tgtgCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATgcgc                                                                                                                    <  1:1008728/149‑1 (MQ=255)
                                    tgCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTa                                                                                                                  >  3:1469387/1‑149 (MQ=255)
                                     gCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTAc                                                                                                                 >  3:1073421/1‑149 (MQ=255)
                                     gCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTAc                                                                                                                 >  3:1255297/1‑149 (MQ=255)
                                      cTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACg                                                                                                                >  1:38303/1‑149 (MQ=255)
                                       tGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGc                                                                                                               >  2:732641/1‑149 (MQ=255)
                                           tatGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCtcat                                                                                                           <  3:397999/149‑1 (MQ=255)
                                             tgctgcTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCtcatca                                                                                                         <  2:327463/149‑1 (MQ=255)
                                              gctgctGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAg                                                                                                        <  2:12054/149‑1 (MQ=255)
                                              gctgctGTTCTACGCCAGTCCTAAGGCATCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAg                                                                                                        >  3:991814/1‑149 (MQ=255)
                                               ctgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAgg                                                                                                       <  3:503720/149‑1 (MQ=255)
                                                tgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑g                                                                                                                                                     <  3:622679/102‑1 (MQ=255)
                                                tgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTga                                                                                                                                          <  1:1395869/113‑1 (MQ=255)
                                                tgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTga                                                                                                                                          >  2:1395869/1‑113 (MQ=255)
                                                tgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGc                                                                                                      >  3:932083/1‑149 (MQ=255)
                                                tgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGc                                                                                                      >  2:969872/1‑149 (MQ=255)
                                                 gctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCa                                                                                                     >  2:586864/1‑149 (MQ=255)
                                                   tgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCtcat                                                                                                           <  1:1367209/141‑1 (MQ=255)
                                                     ttCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACa                                                                                                 <  4:1342592/149‑1 (MQ=255)
                                                           gCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGa                                                                                                                         >  1:983647/1‑119 (MQ=255)
                                                             cAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTcc                                                                                         <  2:861458/149‑1 (MQ=255)
                                                                tCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGcc                                                                                      <  3:322542/149‑1 (MQ=255)
                                                                 ccTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCa                                                                                     <  4:932083/149‑1 (MQ=255)
                                                                   tAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCAGCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTGTAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACt                                                                                   >  3:1246153/1‑149 (MQ=255)
                                                                   tAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACt                                                                                   >  2:853959/1‑149 (MQ=255)
                                                                   tAAGGCGGCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACt                                                                                   >  1:1455318/1‑149 (MQ=255)
                                                                       gCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACAt                                                                               >  2:1122794/1‑149 (MQ=255)
                                                                         gtcgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTa                                                                             >  1:398385/1‑149 (MQ=255)
                                                                         gtcgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTa                                                                             >  4:195746/1‑149 (MQ=255)
                                                                         gtcgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTa                                                                             >  3:1256344/1‑149 (MQ=255)
                                                                          tcgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTaa                                                                            <  3:1124192/149‑1 (MQ=255)
                                                                          tcgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTaa                                                                            <  4:554649/149‑1 (MQ=255)
                                                                          tcgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTaa                                                                            <  2:1100863/149‑1 (MQ=255)
                                                                          tcgtACTGCAACGGACAGTCGATCAATGTCACAGGCGGTCAGGTGAGTTTCTGATCAACAGCGGAGATCGATTAA‑GGATCTCAGTGAGACTATAGATTGCCTGATCCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTaa                                                                            <  2:979263/149‑1 (MQ=255)
                                                                           cgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAg                                                                           <  4:819070/149‑1 (MQ=255)
                                                                                tGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaaaa                                                                      <  1:108195/149‑1 (MQ=255)
                                                                                tGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaaaa                                                                      <  2:310852/149‑1 (MQ=255)
                                                                                    ccGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGCCTTCCGCCACTACATTAAGGAAAAGTTa                                                                  <  4:857052/149‑1 (MQ=255)
                                                                                           gTCGAGCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGGTCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATc                                                           >  2:1350992/1‑149 (MQ=255)
                                                                                            tCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATcc                                                          >  4:792298/1‑149 (MQ=255)
                                                                                              gATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGc                                                        >  4:1166345/1‑149 (MQ=255)
                                                                                               aTCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTCTGGTATCCGCa                                                       <  3:1166345/149‑1 (MQ=255)
                                                                                                     gTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACtcatc                                                 <  3:195746/149‑1 (MQ=255)
                                                                                                       cACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCAc                                               >  3:1327790/1‑149 (MQ=255)
                                                                                                        aCCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGGTGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCAcc                                              >  1:566922/1‑149 (MQ=255)
                                                                                                        aCCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATAc                                                                                                  <  4:1097776/97‑1 (MQ=255)
                                                                                                        aCCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCAcc                                              >  1:174996/1‑149 (MQ=255)
                                                                                                          cggcggTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGt                                            <  1:1152641/149‑1 (MQ=255)
                                                                                                             cggTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGc                                         <  2:1044993/149‑1 (MQ=255)
                                                                                                               gacAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGCCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATTGTATCCGCACTCATCACCGTCGcct                                       >  4:981733/3‑148 (MQ=255)
                                                                                                                  aGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTa                                    <  4:326407/149‑1 (MQ=255)
                                                                                                                    gTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATc                                  <  2:174996/149‑1 (MQ=255)
                                                                                                                    gTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATc                                  >  1:862728/1‑149 (MQ=255)
                                                                                                                     tGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTACCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTAtag                                 >  4:1274217/1‑147 (MQ=255)
                                                                                                                       aTGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGcc                               >  1:305607/1‑149 (MQ=255)
                                                                                                                        tGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCt                              <  4:253788/149‑1 (MQ=255)
                                                                                                                              gATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGtgtg                        >  1:1205978/1‑149 (MQ=255)
                                                                                                                               aTCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCt                                                                                                                           >  4:123560/1‑49 (MQ=255)
                                                                                                                                tCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTgcg                      >  3:1407554/1‑149 (MQ=255)
                                                                                                                                   aCAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTgc                       >  2:991335/1‑145 (MQ=255)
                                                                                                                                   aCAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTgc                       <  1:991335/145‑1 (MQ=255)
                                                                                                                                   aCAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCaa                   <  3:825090/149‑1 (MQ=255)
                                                                                                                                    tcacGGAGATCCTTAAT‑GGATCTCCGCGAGACTATAGAATGCCTGATGCGCTACGCTCATCAgg                                                                                                       <  3:1068668/61‑1 (MQ=21)
                                                                                                                                    tcacGGAGATCCTTAAT‑GGATCTCCGCGAGACTATAGAATGCCTGATGCGCTACGCTCATCAgg                                                                                                       >  4:1068668/4‑64 (MQ=21)
                                                                                                                                    cAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAAc                  <  1:632760/149‑1 (MQ=255)
                                                                                                                                      acGGAGATCCTTAAT‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaaaa                                                                      <  1:700166/94‑1 (MQ=39)
                                                                                                                                      acGGAGATCCTTAAT‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaaaa                                                                      >  2:700166/2‑95 (MQ=39)
                                                                                                                                        ggAGATCCTTAAT‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATGGCCTGGTGTGCGCAACTGGc              >  1:880428/1‑149 (MQ=255)
                                                                                                                                        ggAGATCC‑TTAATGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGc              >  1:710098/1‑149 (MQ=255)
                                                                                                                                            aTCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTa          >  3:1462431/1‑149 (MQ=255)
                                                                                                                                             tCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACtc                                                    >  3:698557/1‑106 (MQ=255)
                                                                                                                                             tCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACtc                                                    <  4:698557/106‑1 (MQ=255)
                                                                                                                                             tCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGgtgt                         >  4:1433442/1‑133 (MQ=255)
                                                                                                                                                aTTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggcg      >  2:59963/1‑149 (MQ=255)
                                                                                                                                                aTTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggcg      >  2:269943/1‑149 (MQ=255)
                                                                                                                                                aTTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCCTTATCGCCTGGTGTGCGCAACTGGCCTTAggcg      >  2:404617/1‑149 (MQ=255)
                                                                                                                                                 ttAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggcgg     >  2:1162425/1‑149 (MQ=255)
                                                                                                                                                 ttAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggcgg     <  1:140637/149‑1 (MQ=255)
                                                                                                                                                    a‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTg  >  3:1356071/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTG  >  NC_000913/2828389‑2828685

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: ‑