Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,828,537 A→T 5.8% 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.8% 359.8 / 7.2 118intergenic (+21/‑84)srlD/srlMsorbitol‑6‑phosphate dehydrogenase/sorbitol=responsive srl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (65/45);  new base T (5/2);  total (70/48)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 119 total aligned reads displayed.

AAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTG  >  NC_000913/2828389‑2828685
                                                                                                                                                    |                                                                                                                                                     
aaGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑g                                                                                                                                                     <  1:417817/150‑1 (MQ=255)
    accgagcAAACGCGGCTGAGATTATCAAGATGTGCTGAATATGCTTCTGGTCTACGCCAGTCCTAAGGCGTCGTACTGCTCGGGACAGTCGAGCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGAtc                                                                                                                                                 <  2:1590081/144‑1 (MQ=255)
    aCCGCTCAACCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGTTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGAtc                                                                                                                                                 <  2:896781/150‑1 (MQ=255)
      cGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGAtctc                                                                                                                                               >  2:1056517/1‑150 (MQ=255)
        cTCAAACGCGGCTGCGATTATCAAGATGTTCTTAATATGCTGCTGTTCTACGCCAGTCCTAATGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCg                                                                                                                                             <  1:1495997/150‑1 (MQ=255)
         tCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCTTAAT‑GGATCTCCgc                                                                                                                                            <  1:253866/150‑2 (MQ=255)
          cAAACGCGGCTGCGATTATCTAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCTTAAT‑GGATCTCCgct                                                                                                                                           <  2:1387196/150‑3 (MQ=255)
            aaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTgag                                                                                                                                         <  1:1129706/150‑1 (MQ=255)
            aaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTgag                                                                                                                                         <  2:506754/150‑1 (MQ=255)
                cgGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACta                                                                                                                                     >  1:614763/1‑150 (MQ=255)
                     gCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGaa                                                                                                                                >  2:1195807/1‑150 (MQ=255)
                         ttATCAAGATGTGCTGAATTTGCGGCTGTTCTACGCCAGGCCTAAGGCGGCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGcc                                                                                                                            <  2:1227875/150‑1 (MQ=255)
                         ttATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTTCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGcc                                                                                                                            <  2:1115687/150‑1 (MQ=255)
                           aTCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTgag                                                                                                                                         >  1:390120/1‑135 (MQ=255)
                           aTCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTgag                                                                                                                                         <  2:390120/135‑1 (MQ=255)
                            tCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGa                                                                                                                         >  1:1760561/1‑150 (MQ=255)
                              aaGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATg                                                                                                                       >  1:481425/1‑150 (MQ=255)
                               aGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATgc                                                                                                                      >  2:1628710/1‑150 (MQ=255)
                               aGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATgc                                                                                                                      <  1:228805/150‑1 (MQ=255)
                                gATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATgcg                                                                                                                     >  1:978356/1‑150 (MQ=255)
                                gATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGCCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATgcg                                                                                                                     >  2:664543/1‑150 (MQ=255)
                                 aTGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATgcgc                                                                                                                    <  1:601229/150‑1 (MQ=255)
                                   gtgCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTa                                                                                                                  >  2:1367874/1‑150 (MQ=255)
                                    tgCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTAc                                                                                                                 >  2:854820/1‑150 (MQ=255)
                                     gCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACg                                                                                                                >  1:1695333/1‑150 (MQ=255)
                                     gCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACg                                                                                                                >  1:1133773/1‑150 (MQ=255)
                                      cTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGc                                                                                                               >  1:1603987/1‑150 (MQ=255)
                                        gAATATGCTGCTGTTCAACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGAGCAATGTCACCGGCGGTCAGGTGAAGTACTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCAGATGCGCTACgata                                                                                                             >  1:1002757/1‑147 (MQ=255)
                                             tgctgcTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAg                                                                                                        <  1:1400017/150‑1 (MQ=255)
                                               ctgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGAGCGTAGCGCATCAGGc                                                                                                      <  1:1012656/150‑1 (MQ=255)
                                                tgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGAGCTACGCTCATCAGGCa                                                                                                     >  2:869005/1‑150 (MQ=255)
                                                 gctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCAt                                                                                                    >  2:375776/1‑150 (MQ=255)
                                                    gTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTga                                                                                                                                          >  2:1570200/1‑109 (MQ=255)
                                                    gTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTga                                                                                                                                          <  1:1570200/109‑1 (MQ=255)
                                                        tACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGAc                                                                                             >  2:1016067/1‑150 (MQ=255)
                                                        tACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGAc                                                                                             >  1:90922/1‑150 (MQ=255)
                                                          cGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTTCGCTCATCAGGCATACAGGACtt                                                                                           >  2:357089/1‑150 (MQ=255)
                                                            ccAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTcc                                                                                         >  2:899789/1‑150 (MQ=255)
                                                                 ccTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCAc                                                                                    >  1:324005/1‑150 (MQ=255)
                                                                  cTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACt                                                                                   <  1:1703234/150‑1 (MQ=255)
                                                                    aaGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTAc                                                                                 <  2:481425/150‑1 (MQ=255)
                                                                        cgtcgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTa                                                                             <  1:1641272/150‑1 (MQ=255)
                                                                        cgtcgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTa                                                                             <  2:1760561/150‑1 (MQ=255)
                                                                           cgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAgg                                                                          <  2:1296995/150‑1 (MQ=255)
                                                                           cgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAgg                                                                          >  2:1438116/1‑150 (MQ=255)
                                                                           cgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAgg                                                                          <  2:1424554/150‑1 (MQ=255)
                                                                              aCTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaaa                                                                       >  2:290499/1‑150 (MQ=255)
                                                                              aCTGCACCGGACAGTCGATCAATGTCACCGGCGGGCAGGTGATGTGCTGATCAACAGCGGAGATCCATTAC‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaaa                                                                       <  1:38578/150‑1 (MQ=255)
                                                                               cTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaaaa                                                                      >  1:1064309/1‑150 (MQ=255)
                                                                                 gCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGt                                                                    <  1:762001/150‑1 (MQ=255)
                                                                                           gTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCTGAGATCCATTAA‑GGATCTCCGTGAGACTATATAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATcc                                                          <  1:871309/150‑1 (MQ=255)
                                                                                           gTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATcc                                                          >  1:811649/1‑150 (MQ=255)
                                                                                           gTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATcc                                                          >  2:1066562/1‑150 (MQ=255)
                                                                                              gATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCa                                                       >  2:691817/1‑150 (MQ=255)
                                                                                              gATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCa                                                       >  1:1692920/1‑150 (MQ=255)
                                                                                              gATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCa                                                       >  1:58071/1‑150 (MQ=255)
                                                                                                tCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACt                                                     <  2:1010764/150‑1 (MQ=255)
                                                                                                     gTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGGCACTACATTAAGGAAAAGTTATGGTATCCGCACtcatca                                                >  2:927337/1‑150 (MQ=255)
                                                                                                     gTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACtcatca                                                <  2:389016/150‑1 (MQ=255)
                                                                                                     gTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACtcatca                                                >  2:1149293/1‑150 (MQ=255)
                                                                                                        aCCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCg                                             >  2:1555295/1‑150 (MQ=255)
                                                                                                         ccGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGt                                            >  1:1008941/1‑150 (MQ=255)
                                                                                                         ccGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGt                                            >  1:458188/1‑150 (MQ=255)
                                                                                                             cggTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGcc                                        <  1:622834/150‑1 (MQ=255)
                                                                                                              ggTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCg                                       <  2:222369/150‑1 (MQ=255)
                                                                                                               gTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGt                                      <  2:902326/150‑1 (MQ=255)
                                                                                                               gTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGt                                      <  1:425858/150‑1 (MQ=255)
                                                                                                                 cAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATTCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTa                                    <  2:1127676/150‑1 (MQ=255)
                                                                                                                 cAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTa                                    >  2:846814/1‑150 (MQ=255)
                                                                                                                   ggTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATc                                  <  2:1599613/150‑1 (MQ=255)
                                                                                                                         gTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTgg                            >  1:1385948/1‑150 (MQ=255)
                                                                                                                             tGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGtgtg                        >  1:628697/1‑150 (MQ=255)
                                                                                                                             tGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGtgtg                        >  2:377317/1‑150 (MQ=255)
                                                                                                                             tGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGtgtg                        <  1:1195807/150‑1 (MQ=255)
                                                                                                                              gATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTgc                       >  2:851135/1‑150 (MQ=255)
                                                                                                                              gATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTgc                       >  2:1095025/1‑150 (MQ=255)
                                                                                                                                 cAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCa                    >  1:137451/1‑150 (MQ=255)
                                                                                                                                 cAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCa                    >  1:434013/1‑150 (MQ=255)
                                                                                                                                  tctcacGGAGAT‑CATTAATGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCAGCAGGCATACAGGACTTCCGCAACTACATTAAGGAAAAGTTATGGTAACCGCACTCATCACCGTCGCCCTTAACGCCTGGAGTGCGCaa                   >  2:175408/6‑150 (MQ=255)
                                                                                                                                   aCAGCGGAGATCCATTAA‑TGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAAc                  <  2:1182649/150‑1 (MQ=255)
                                                                                                                                    cAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACt                 >  2:1515048/1‑150 (MQ=255)
                                                                                                                                     cacGGAGATCCTTAAT‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTg                >  2:591390/3‑150 (MQ=255)
                                                                                                                                     cacGGAGATCCTTAAT‑GGATCTCCGCGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTg                >  1:827772/3‑150 (MQ=255)
                                                                                                                                     aGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTg                >  1:470142/1‑150 (MQ=255)
                                                                                                                                      gCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTgg               <  1:1056517/150‑1 (MQ=255)
                                                                                                                                      gCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACAGTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCTCCTGGTGTGCGCAACTgg               >  1:450359/1‑150 (MQ=255)
                                                                                                                                      acGGAGATCCTTAAT‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTgg               >  1:557085/2‑150 (MQ=255)
                                                                                                                                        ggAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGcc             <  2:324005/150‑1 (MQ=255)
                                                                                                                                        ggAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGcc             <  1:290499/150‑1 (MQ=255)
                                                                                                                                           gaTCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTa          >  2:702261/1‑150 (MQ=255)
                                                                                                                                           caTCCAATAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAACGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTa          <  1:282312/149‑1 (MQ=255)
                                                                                                                                             tCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAgg        <  1:1613649/150‑1 (MQ=255)
                                                                                                                                               cATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggcg      >  2:713912/1‑150 (MQ=255)
                                                                                                                                               cATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggcg      >  1:21275/1‑150 (MQ=255)
                                                                                                                                               cATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggcg      >  2:1460598/1‑150 (MQ=255)
                                                                                                                                               cATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggcg      >  2:1198002/1‑150 (MQ=255)
                                                                                                                                               cATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGATATGTTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggcg      >  1:1612949/1‑150 (MQ=255)
                                                                                                                                                  tAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCt   <  1:1742857/150‑1 (MQ=255)
                                                                                                                                                  tAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCt   >  1:993502/1‑150 (MQ=255)
                                                                                                                                                   aa‑GGTTCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATAGAGGACTTGCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTg  <  2:1002757/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAA‑GGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTG  >  NC_000913/2828389‑2828685

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