Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,652,372 C→T 5.2% intergenic (‑85/‑122) yfhM ← / → sseA bacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein/3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,652,3720CT 5.2% 257.0 / 2.6 76intergenic (‑85/‑122)yfhM/sseAbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein/3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base T (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

GTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGC  >  NC_000913/2652228‑2652519
                                                                                                                                                |                                                                                                                                                   
gttgtCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCATCGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCAt                                                                                                                                                 <  2:663694/149‑1 (MQ=255)
gttgtCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCAt                                                                                                                                                 <  4:235882/149‑1 (MQ=255)
      gCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATc                                                                                                                                           <  1:111822/149‑1 (MQ=255)
      gCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATc                                                                                                                                           >  2:406031/1‑149 (MQ=255)
      gCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAACTTACGTTATGTAATTCAGTATCACGCCCATACAATc                                                                                                                                           <  1:299877/149‑1 (MQ=255)
       cACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATcc                                                                                                                                          >  3:516380/1‑149 (MQ=255)
        aCCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGTCCATACAATCCg                                                                                                                                         >  1:509043/1‑149 (MQ=255)
        aCCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGTCCATACAATCCg                                                                                                                                         >  1:514498/1‑149 (MQ=255)
         cccTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGa                                                                                                                                        >  3:334818/1‑149 (MQ=255)
         cccTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGa                                                                                                                                        >  3:570943/1‑149 (MQ=255)
          ccTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGaa                                                                                                                                       <  1:406031/149‑1 (MQ=255)
           cTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAAt                                                                                                                                      <  4:783776/149‑1 (MQ=255)
             gccagcgccagcATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAAttt                                                                                                                                    <  4:334818/149‑1 (MQ=255)
               cagcgccagcATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGTCCATACAATCCGAATTTTa                                                                                                                                  <  2:509043/149‑1 (MQ=255)
               cagcgccagcATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGTCCATACAATCCGAATTTTa                                                                                                                                  <  2:514498/149‑1 (MQ=255)
               cagcgccagcATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTa                                                                                                                                  >  2:725408/1‑149 (MQ=255)
               cagcgccagcATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTa                                                                                                                                  <  4:648473/149‑1 (MQ=255)
                agcgccagcATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAg                                                                                                                                 <  4:516380/149‑1 (MQ=255)
                 gcgccagcATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGc                                                                                                                                                     >  4:222266/1‑128 (MQ=255)
                 gcgccagcATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGc                                                                                                                                                     <  3:222266/128‑1 (MQ=255)
                          ttAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATtgt                                                                                                                       <  2:518633/149‑1 (MQ=255)
                           tAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATtgtg                                                                                                                      <  2:481948/149‑1 (MQ=255)
                             gCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGc                                                                                                                    >  3:551894/1‑149 (MQ=255)
                                  cAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGa                                                                                                               >  4:126071/1‑149 (MQ=255)
                                       ggcTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGcc                                                                                                          <  4:464156/149‑1 (MQ=255)
                                          tACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCaaa                                                                                                       >  4:507329/1‑149 (MQ=255)
                                              cgTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAg                                                                                                   <  2:279228/149‑1 (MQ=255)
                                                 aaCTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAg                                                                                                            >  3:718661/1‑137 (MQ=255)
                                                 aaCTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAg                                                                                                            <  4:718661/137‑1 (MQ=255)
                                                        tcatttttcatCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGc                                                                                         <  2:807158/149‑1 (MQ=255)
                                                          atttttcatCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCtt                                                                                       <  3:464183/149‑1 (MQ=255)
                                                          atttttcatCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCtt                                                                                       >  1:447601/1‑149 (MQ=255)
                                                            ttttcatCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGt                                                                                     <  2:447601/149‑1 (MQ=255)
                                                                   cccTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCtttt                                                                              <  3:380889/149‑1 (MQ=255)
                                                                   cccTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCtttt                                                                              <  1:152202/149‑1 (MQ=255)
                                                                    ccTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTc                                                                             >  1:460562/1‑149 (MQ=255)
                                                                    ccTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTc                                                                             <  1:407216/149‑1 (MQ=255)
                                                                      tGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCAc                                                                           >  2:1840/1‑149 (MQ=255)
                                                                          acaATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGcc                                                                                                          >  4:100682/1‑114 (MQ=255)
                                                                          acaATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTaaa                                                                       >  4:713158/1‑149 (MQ=255)
                                                                          acaATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAAAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGcc                                                                                                          <  3:100682/114‑1 (MQ=255)
                                                                             aTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACAt                                                                    >  4:711290/1‑149 (MQ=255)
                                                                                gCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCg                                                                 >  1:293526/1‑149 (MQ=255)
                                                                                 cccGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGa                                                                >  3:309207/1‑149 (MQ=255)
                                                                                  ccGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCttt                                                                               >  1:29557/1‑133 (MQ=255)
                                                                                  ccGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCttt                                                                               <  2:29557/133‑1 (MQ=255)
                                                                                    gTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACaa                                                             <  3:16739/149‑1 (MQ=255)
                                                                                     ttAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAAt                                                            <  3:126071/149‑1 (MQ=255)
                                                                                             gTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGcc                                                                                                 >  1:594897/1‑104 (MQ=255)
                                                                                             gTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGcc                                                                                                 <  2:594897/104‑1 (MQ=255)
                                                                                                   taatTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCg                                              >  2:531429/1‑149 (MQ=255)
                                                                                                     attattATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGt                                                                                     >  2:522701/1‑108 (MQ=255)
                                                                                                     attattATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGt                                                                                     <  1:522701/108‑1 (MQ=255)
                                                                                                            tGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATg                                     <  3:711290/149‑1 (MQ=255)
                                                                                                              cGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGga                                   <  4:570943/149‑1 (MQ=255)
                                                                                                              cGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGga                                   <  1:1840/149‑1 (MQ=255)
                                                                                                                tGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATtgtg                                                                                                                      <  2:745542/64‑1 (MQ=255)
                                                                                                                tGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATtgtg                                                                                                                      >  1:745542/1‑64 (MQ=255)
                                                                                                                tGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGgaga                                 >  2:303256/1‑149 (MQ=255)
                                                                                                                 gagaATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGAt                                >  1:763228/1‑149 (MQ=255)
                                                                                                                 gagaATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGAt                                >  2:385184/1‑149 (MQ=255)
                                                                                                                 gagaATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTAGCTTTTCACTAAACATCCGACAATTTAGCCTAACTCCGGCAAAAATGGAGAt                                >  4:286441/1‑149 (MQ=255)
                                                                                                                       ttACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTAt                          <  3:713158/149‑1 (MQ=255)
                                                                                                                         aCGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGt                        <  2:460562/149‑1 (MQ=255)
                                                                                                                            ttATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCa                     >  2:490601/1‑149 (MQ=255)
                                                                                                                                 tAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACa                >  4:553941/1‑149 (MQ=255)
                                                                                                                                   aTTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATg              >  4:427816/1‑149 (MQ=255)
                                                                                                                                   aTTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATg              >  3:26483/1‑149 (MQ=255)
                                                                                                                                    ttCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATgg             <  4:309207/149‑1 (MQ=255)
                                                                                                                                       aGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGttt          >  3:487535/1‑149 (MQ=255)
                                                                                                                                       aGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGttt          >  1:372397/1‑149 (MQ=255)
                                                                                                                                         tATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCttt                                                                                            <  4:294393/65‑1 (MQ=255)
                                                                                                                                         tATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCttt                                                                                            >  3:294393/1‑65 (MQ=255)
                                                                                                                                           tCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATg              <  3:60174/141‑1 (MQ=255)
                                                                                                                                           tCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATg              >  4:60174/1‑141 (MQ=255)
                                                                                                                                           tCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAg      >  3:239942/1‑149 (MQ=255)
                                                                                                                                             aCGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGa    <  2:763228/149‑1 (MQ=255)
                                                                                                                                               gCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGc  >  1:461768/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGC  >  NC_000913/2652228‑2652519

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