Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,684,085 A→C 5.2% intergenic (+29/+169) yphH → / ← glyA putative DNA‑binding transcriptional regulator/serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,684,0850AC 5.2% 360.2 / 5.6 104intergenic (+29/+169)yphH/glyAputative DNA‑binding transcriptional regulator/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (44/54);  new base C (3/3);  total (47/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 108 total aligned reads displayed.

GACCGCGACGAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACA  >  NC_000913/2683937‑2684230
                                                                                                                                                    |                                                                                                                                                 
gACCGCGACGAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCga                                                                                                                                                   >  4:494597/1‑149 (MQ=255)
 aCCGCGACGAAGGACCGAGCGTGAAACCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCCCAACTGGTGCTGGGCATTGGCTATTTGTATGTTGCGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCgaa                                                                                                                                                  >  2:237163/1‑149 (MQ=255)
  ccGCGACGAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCgaag                                                                                                                                                 <  3:236993/149‑1 (MQ=255)
  ccGCGACGAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCgaag                                                                                                                                                 <  1:526443/149‑1 (MQ=255)
         gaAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCTTCGGGAAGCCt                                                                                                                                          <  4:50604/149‑1 (MQ=255)
           aGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCTTCGGGAAGCCttt                                                                                                                                        <  2:281603/149‑2 (MQ=255)
           aGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTtg                                                                                                                                        >  1:301447/1‑149 (MQ=255)
           aGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTtg                                                                                                                                        >  1:386800/1‑149 (MQ=255)
            ggACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTtga                                                                                                                                       >  3:649437/1‑149 (MQ=255)
            ggACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTtga                                                                                                                                       <  1:411722/149‑1 (MQ=255)
             gACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAAGAAGCCTtgat                                                                                                                                      <  3:571728/149‑1 (MQ=255)
               ccGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTtgatga                                                                                                                                    <  1:231521/149‑1 (MQ=255)
               ccGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTtgatga                                                                                                                                    <  1:817611/149‑1 (MQ=255)
               ccGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTtgatga                                                                                                                                    <  3:317016/149‑1 (MQ=255)
                cGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCGTgatgat                                                                                                                                   <  2:467442/149‑1 (MQ=255)
                 gAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATc                                                                                                                                  <  1:695108/149‑1 (MQ=255)
                      tGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATaaa                                                                                                                             <  3:335614/149‑1 (MQ=255)
                      tGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATaaa                                                                                                                             <  4:369671/149‑1 (MQ=255)
                       gAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAAc                                                                                                                            <  2:624391/149‑1 (MQ=255)
                           gCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGaaa                                                                                                                        <  1:483143/149‑1 (MQ=255)
                            cGACGCAAATTGGCTTTGGGCAATTAAGCGGCGCACAACAGGTGATGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCTTGGATGATCATAAACGaaaa                                                                                                                       >  2:692841/1‑149 (MQ=255)
                            cGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGaaaa                                                                                                                       >  2:77296/1‑149 (MQ=255)
                               cGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAAtt                                                                                                                    <  1:654518/149‑1 (MQ=255)
                                gCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTg                                                                                                                   >  3:200077/1‑149 (MQ=255)
                                    aTTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTg                                                                                                               >  2:403179/1‑149 (MQ=255)
                                     ttGGCTTTGGGCAATTAAGCCGCGCAGAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGa                                                                                                              <  2:603517/149‑1 (MQ=255)
                                     ttGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGa                                                                                                              <  3:806423/149‑1 (MQ=255)
                                         cTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATgcg                                                                                                          >  1:393836/1‑149 (MQ=255)
                                          tttGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATgcgc                                                                                                         >  3:298741/1‑149 (MQ=255)
                                          tttGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATgcgc                                                                                                         >  1:682117/1‑149 (MQ=255)
                                             gggCAATTAAGCCGCGCACAACGGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTAc                                                                                                      <  2:333234/149‑1 (MQ=255)
                                             gggCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTAc                                                                                                      <  1:166821/149‑1 (MQ=255)
                                               gCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGc                                                                                                    <  3:351328/149‑1 (MQ=255)
                                                 aaTTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCtt                                                                                                  >  3:346148/1‑149 (MQ=255)
                                                  aTTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTa                                                                                                 >  3:578703/1‑149 (MQ=255)
                                                   ttAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTAt                                                                                                <  1:207487/149‑1 (MQ=255)
                                                      aGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAg                                                                                             <  3:204692/149‑1 (MQ=255)
                                                       gCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAgg                                                                                            <  4:624227/149‑1 (MQ=255)
                                                        ccGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGc                                                                                           <  3:800181/149‑1 (MQ=255)
                                                              acaacaGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTacac                                                                                     >  4:816085/1‑149 (MQ=255)
                                                                 acaGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                                  <  4:744112/149‑1 (MQ=255)
                                                                 acaGGGGCTGGGAATTGGTTATTTGGATGTTGGGGCGCAGTTACGACAGATTTGATGGCGCGATAATGTAGAAAGGCTTCCCGCAGGAAGCCTTGATGATCAGAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                                  <  4:344529/149‑1 (MQ=38)
                                                                  caGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGgag                                                                                 <  3:775417/149‑1 (MQ=255)
                                                                     gTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGAtt                                                                              <  4:200077/149‑1 (MQ=255)
                                                                     gTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGAtt                                                                              <  3:462550/149‑1 (MQ=255)
                                                                     gTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGAtt                                                                              <  2:823520/149‑1 (MQ=255)
                                                                      tGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTg                                                                             >  2:403824/1‑149 (MQ=255)
                                                                         tGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCaa                                                                          <  1:612043/149‑1 (MQ=255)
                                                                          gggAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGc                                                                                                                  >  3:85565/1‑108 (MQ=255)
                                                                          gggAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAat                                                                         >  1:728596/1‑149 (MQ=255)
                                                                            gAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAatat                                                                       <  3:816085/149‑1 (MQ=255)
                                                                                  gCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGaa                                                                 >  1:228695/1‑149 (MQ=255)
                                                                                  gCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGaa                                                                 <  1:403824/149‑1 (MQ=255)
                                                                                   cTATTTGGATGTTGAGGCGCCGTTGCGACAGTTTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAAt                                                                >  4:728103/1‑149 (MQ=255)
                                                                                    tATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAAtt                                                               >  2:36098/1‑149 (MQ=255)
                                                                                    tATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAAtt                                                               >  2:690839/1‑149 (MQ=255)
                                                                                      tttGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTg                                                             >  1:94129/1‑149 (MQ=255)
                                                                                       ttGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATCACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGTTCATAAACGAGAAATTGCCTGATGCGCTACGCTTATCAGGCCTATACGGTGATTGCAATATATTGAATTTgg                                                            >  2:643729/1‑148 (MQ=38)
                                                                                          tATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCaaa                                                         <  2:69631/149‑1 (MQ=255)
                                                                                            tGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGt                                                       >  1:154345/1‑149 (MQ=255)
                                                                                            tGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGt                                                       >  4:531675/1‑149 (MQ=255)
                                                                                             gTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGtt                                                      >  1:764216/1‑149 (MQ=255)
                                                                                               tGAGGCGCAGTTACTACAGATTGGATGGCTCGATAACGTAGAAAGGCTTCCAGAAGGAAACCTTGATGATCATTAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGtttt                                                    <  1:643729/149‑1 (MQ=38)
                                                                                               tGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGTTGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGtttt                                                    <  4:498747/149‑1 (MQ=255)
                                                                                                 aGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGt                                                  <  4:668178/149‑1 (MQ=255)
                                                                                                     gcAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCTATATATTGAATTTGCAAAGTTTTGTAGGc                                              >  4:155542/1‑149 (MQ=255)
                                                                                                           aCGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa                                        <  1:584583/149‑1 (MQ=255)
                                                                                                            cGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATaa                                       <  4:367994/149‑1 (MQ=255)
                                                                                                            cGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATaa                                       >  4:487728/1‑149 (MQ=255)
                                                                                                            cGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATaa                                       <  3:283351/149‑1 (MQ=255)
                                                                                                            cGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATaa                                       <  2:506576/149‑1 (MQ=255)
                                                                                                             gACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGCTTGCAATATATTGAATTTGCAAAGGTTTGTATGCCGGATAAg                                      >  2:432797/1‑149 (MQ=255)
                                                                                                             gACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                                      <  3:713901/149‑1 (MQ=255)
                                                                                                              aCAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAgg                                     >  4:278188/1‑149 (MQ=255)
                                                                                                               cAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGc                                    <  3:660282/149‑1 (MQ=255)
                                                                                                                aGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGtttt                                                    >  2:423347/1‑132 (MQ=255)
                                                                                                                aGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGtttt                                                    <  1:423347/132‑1 (MQ=255)
                                                                                                                   tttGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTc                                >  3:571274/1‑149 (MQ=255)
                                                                                                                                 aaCGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                  <  4:346148/149‑1 (MQ=37)
                                                                                                                                  aCGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                 <  3:586971/149‑1 (MQ=255)
                                                                                                                                  aCGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                 >  3:748433/1‑149 (MQ=255)
                                                                                                                                    gTATAAAGGCTTCCCGAATGAAGGCAAGATGATCATAAAAGAAAAATTGTCTGATGCGCTACGCTTATCAGGCGTACTCGGCGATTGCAAGATATTGAATTCGCAAAGTTTTGTAGGCCGGATAAGGCGTTTGGGCCGCATACGGTATg               <  2:407715/149‑1 (MQ=18)
                                                                                                                                    gTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                      >  4:770367/1‑142 (MQ=255)
                                                                                                                                    gTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                      <  3:770367/142‑1 (MQ=255)
                                                                                                                                     tAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa              >  4:352152/1‑149 (MQ=255)
                                                                                                                                      atcaaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa             >  2:140433/4‑149 (MQ=21)
                                                                                                                                      atcaaGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa             >  3:766080/4‑149 (MQ=255)
                                                                                                                                      aGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa             >  3:766954/1‑149 (MQ=255)
                                                                                                                                       tcaaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac            >  4:16214/3‑149 (MQ=255)
                                                                                                                                        caaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGtttt                                                    >  2:93306/2‑108 (MQ=25)
                                                                                                                                          aGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac         <  4:591509/149‑1 (MQ=255)
                                                                                                                                          aGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac         <  4:298741/149‑1 (MQ=255)
                                                                                                                                          aGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac         <  4:649437/149‑1 (MQ=255)
                                                                                                                                             cTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCGACACGGAGATTGCAATATATTGAATTAGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAg      <  2:514873/149‑1 (MQ=255)
                                                                                                                                              ttCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGc     >  1:217381/1‑149 (MQ=255)
                                                                                                                                              ttCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTAACAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGc     >  4:398128/1‑149 (MQ=255)
                                                                                                                                                cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcac   >  2:631056/1‑149 (MQ=255)
                                                                                                                                                 ccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcaca  >  2:355050/1‑149 (MQ=255)
                                                                                                                                                  cGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATc                                                                                               <  3:578364/55‑1 (MQ=255)
                                                                                                                                                  cGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATc                                                                                               >  4:578364/1‑55 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GACCGCGACGAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACA  >  NC_000913/2683937‑2684230

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