Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,828,537 A→T 5.5% 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.5% 317.0 / 7.0 92intergenic (+21/‑84)srlD/srlMsorbitol‑6‑phosphate dehydrogenase/sorbitol=responsive srl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (32/54);  new base T (3/2);  total (36/56)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

AAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTG  >  NC_000913/2828389‑2828685
                                                                                                                                                    |                                                                                                                                                    
aaGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTaa                                                                                                                                                      <  1:800775/149‑1 (MQ=255)
aaGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGAACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTaa                                                                                                                                                      <  4:718153/149‑1 (MQ=255)
 aGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAg                                                                                                                                                     <  3:243621/149‑1 (MQ=255)
 aGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAg                                                                                                                                                     <  2:340851/149‑1 (MQ=255)
    aCCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGAt                                                                                                                                                  >  1:176772/1‑149 (MQ=255)
        cTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCTTAATGGATCTcc                                                                                                                                              <  2:368343/149‑1 (MQ=255)
         tCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCg                                                                                                                                             >  1:57167/1‑149 (MQ=255)
          cAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGt                                                                                                                                            <  4:145944/149‑1 (MQ=255)
          cAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGt                                                                                                                                            <  4:118517/149‑1 (MQ=255)
          cAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGt                                                                                                                                            <  4:259720/149‑1 (MQ=255)
            aaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTga                                                                                                                                          <  3:19335/149‑1 (MQ=255)
            aaCGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTga                                                                                                                                          <  4:45309/149‑1 (MQ=255)
               gcgGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCCGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGAc                                                                                                                                       <  2:502685/149‑1 (MQ=255)
                   cTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACtata                                                                                                                                   <  4:291327/149‑1 (MQ=255)
                   cTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACtata                                                                                                                                   <  2:572976/149‑1 (MQ=255)
                   cTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACtata                                                                                                                                   <  2:852250/149‑1 (MQ=255)
                    tGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAg                                                                                                                                  <  1:296395/149‑1 (MQ=255)
                      cGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGaa                                                                                                                                <  3:444322/149‑1 (MQ=255)
                       gATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAAt                                                                                                                               <  1:875065/149‑1 (MQ=255)
                          tATCAAGATGTGCTGAATATGCTGCTGTTCTACGGCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGCTCAATGTCAACGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGcc                                                                                                                            <  2:283180/149‑1 (MQ=255)
                           aTCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCt                                                                                                                           <  1:306227/149‑1 (MQ=255)
                           aTCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCt                                                                                                                           <  3:400876/149‑1 (MQ=255)
                              aaGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGAt                                                                                                                        <  4:234241/149‑1 (MQ=255)
                                     gCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTAc                                                                                                                 >  3:802107/1‑149 (MQ=255)
                                      cTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACg                                                                                                                <  2:546736/149‑1 (MQ=255)
                                      cTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTACAGAATGCCTGATGCGCTACg                                                                                                                >  3:649711/1‑149 (MQ=255)
                                        gAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCt                                                                                                              >  1:758088/1‑149 (MQ=255)
                                              cctgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAg                                                                                                        <  1:792057/148‑1 (MQ=255)
                                                tgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGc                                                                                                      >  4:427484/1‑149 (MQ=255)
                                                tgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGc                                                                                                      >  4:199135/1‑149 (MQ=255)
                                                    gTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGTGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATg                                                                                                                       >  4:819153/1‑128 (MQ=255)
                                                    gTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGTGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATg                                                                                                                       <  3:819153/128‑1 (MQ=255)
                                                      tCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAg                                                                                                <  2:663509/149‑1 (MQ=255)
                                                      tCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAg                                                                                                <  2:20506/149‑1 (MQ=255)
                                                           gCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACtt                                                                                           >  4:172233/1‑149 (MQ=255)
                                                             cAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGa                                                                                                                         <  2:498482/117‑1 (MQ=255)
                                                             cAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGa                                                                                                                         >  1:498482/1‑117 (MQ=255)
                                                             cAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTcc                                                                                         >  2:659285/1‑149 (MQ=255)
                                                             cAGTCCTAAGGCGTCGCACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTcc                                                                                         <  3:560260/149‑1 (MQ=255)
                                                                  cTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCAc                                                                                    >  4:180532/1‑149 (MQ=255)
                                                                    aaGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTa                                                                                  <  2:741361/149‑1 (MQ=255)
                                                                     aGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGACGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTAc                                                                                 >  4:716211/1‑149 (MQ=255)
                                                                       gCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACAt                                                                               <  2:791066/149‑1 (MQ=255)
                                                                         gtcgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTa                                                                             <  3:101182/149‑1 (MQ=255)
                                                                           cgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAg                                                                           <  2:570356/149‑1 (MQ=255)
                                                                           cgtACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAg                                                                           <  2:887482/149‑1 (MQ=255)
                                                                             tACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGTGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGa                                                                         >  3:58875/1‑149 (MQ=255)
                                                                             tACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGa                                                                         >  1:804835/1‑149 (MQ=255)
                                                                                tGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTaa                                                                            >  3:86589/1‑143 (MQ=255)
                                                                                tGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTaa                                                                            <  4:86589/143‑1 (MQ=255)
                                                                                tGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaaaa                                                                      <  3:880639/149‑1 (MQ=255)
                                                                                  cACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGt                                                                    >  3:747598/1‑149 (MQ=255)
                                                                                   aCCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGtt                                                                   >  1:375408/1‑149 (MQ=255)
                                                                                    ccGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTa                                                                  <  2:758088/149‑1 (MQ=255)
                                                                                    ccGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTa                                                                  <  2:232779/149‑1 (MQ=255)
                                                                                    ccGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTa                                                                  <  3:180532/149‑1 (MQ=255)
                                                                                     cGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTAt                                                                 <  4:543797/149‑1 (MQ=255)
                                                                                     cGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTAt                                                                 <  2:57167/149‑1 (MQ=255)
                                                                                     cGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTAt                                                                 <  1:172316/149‑1 (MQ=255)
                                                                                      ggACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGTGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATg                                                                <  4:58875/149‑1 (MQ=255)
                                                                                       gACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATgg                                                               >  3:758320/1‑149 (MQ=255)
                                                                                          aGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTAt                                                            >  2:165474/1‑149 (MQ=255)
                                                                                           gTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCtcatca                                                                                                         <  2:731136/103‑1 (MQ=255)
                                                                                           gTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATc                                                           >  2:23172/1‑149 (MQ=255)
                                                                                           gTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGCATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCtcatca                                                                                                         >  1:731136/1‑103 (MQ=255)
                                                                                               aTCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTACGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCa                                                       >  2:224574/1‑149 (MQ=255)
                                                                                               aTCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCa                                                       >  3:351931/1‑149 (MQ=255)
                                                                                               aTCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCa                                                       >  3:23329/1‑149 (MQ=255)
                                                                                                tCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCAc                                                      <  4:23329/149‑1 (MQ=255)
                                                                                                  aaTGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACtc                                                    >  3:686875/1‑149 (MQ=255)
                                                                                                          cggcggTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGt                                            >  4:286191/1‑149 (MQ=255)
                                                                                                             cggTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGc                                         <  2:161921/149‑1 (MQ=255)
                                                                                                                      gATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGc                                <  2:375408/149‑1 (MQ=255)
                                                                                                                      gATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGc                                <  4:691283/149‑1 (MQ=255)
                                                                                                                        tGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCTGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCt                              <  3:699819/149‑1 (MQ=255)
                                                                                                                         gTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTg                             <  1:169268/149‑1 (MQ=255)
                                                                                                                         gTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTg                             <  2:176772/149‑1 (MQ=255)
                                                                                                                             tGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGgtgt                         <  4:686875/149‑1 (MQ=255)
                                                                                                                               aTCAACAGCGGTGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTa                                                                  <  3:762046/106‑1 (MQ=255)
                                                                                                                               aTCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCGTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTGAGGAAAAGTTATGGTATCCGCACTCAGCACCGTCGCCGTTATCGCCTGGTGTgc                       <  1:160979/149‑1 (MQ=255)
                                                                                                                               aTCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTa                                                                  >  4:762046/1‑106 (MQ=255)
                                                                                                                                  aaCAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCa                    >  3:264491/1‑149 (MQ=255)
                                                                                                                                   ctcacGGAGATCCTTAATGGATCTCCGTGAGACTATAGAATGCCTCATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCaa                   >  3:503118/5‑149 (MQ=255)
                                                                                                                                    tcacGGAGATCCTTAATGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaaaa                                                                      >  4:541946/4‑97 (MQ=38)
                                                                                                                                    tcacGGAGATCCTTAATGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaaaa                                                                      <  3:541946/94‑1 (MQ=38)
                                                                                                                                    tcacGGAGATCCTTAATGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAAc                  >  4:264612/4‑149 (MQ=255)
                                                                                                                                        ggAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGc              >  1:748335/1‑149 (MQ=255)
                                                                                                                                         gagaTCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGcc             <  1:165474/149‑1 (MQ=255)
                                                                                                                                            aTCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTa          >  2:572119/1‑149 (MQ=255)
                                                                                                                                             tCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAg         <  1:565894/149‑1 (MQ=255)
                                                                                                                                                 ttAAGGATCTCCGTGAGACTATAGAATGCCTGGTGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTaa                                                                            >  2:619822/1‑78 (MQ=255)
                                                                                                                                                 ttAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTaa                                                                            <  1:619822/78‑1 (MQ=255)
                                                                                                                                                 ttAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggcgg     >  2:305624/1‑149 (MQ=255)
                                                                                                                                                  tAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAAGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggcggc    <  4:337226/149‑1 (MQ=255)
                                                                                                                                                   aaGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCt   <  1:712287/149‑1 (MQ=255)
                                                                                                                                                    aGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTg  >  2:610195/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTG  >  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 < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: