Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 751,117 T→C 5.5% P351P (CCA→CCG ybgQ ← putative fimbrial usher protein YbgQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913751,1170TC 5.5% 274.2 / 4.9 90P351P (CCA→CCGybgQputative fimbrial usher protein YbgQ
Reads supporting (aligned to +/- strand):  ref base T (47/38);  new base C (3/2);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

TATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTC  >  NC_000913/750969‑751263
                                                                                                                                                    |                                                                                                                                                  
tatTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTg                                                                                                                                                   >  2:907750/1‑150 (MQ=255)
tatTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTCGATTTTCAGTTTGATGTGACATGGAAGGCCGTg                                                                                                                                                   >  1:402995/1‑150 (MQ=255)
  ttttGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGc                                                                                                                                                 <  1:822034/150‑1 (MQ=255)
  ttttGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTCACATGGAAGGCCGTGGc                                                                                                                                                 >  2:558612/1‑150 (MQ=255)
       accGATACCAATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGAcc                                                                                                                                            >  2:1595409/2‑150 (MQ=255)
        ccGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGAccc                                                                                                                                           >  2:1469221/1‑150 (MQ=255)
            tACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCgg                                                                                                                                       >  1:370058/1‑150 (MQ=255)
             aCCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGCGGGTCAGCCACGGc                                                                                                                                      <  1:1316158/150‑1 (MQ=255)
             aCCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGCGGGTCAGCCACGGc                                                                                                                                      <  1:1345573/150‑1 (MQ=255)
              cccATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGcc                                                                                                                                     <  2:923342/150‑1 (MQ=255)
               ccATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCa                                                                                                                                    >  1:715319/1‑150 (MQ=255)
               ccATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCa                                                                                                                                    >  2:108042/1‑150 (MQ=255)
                   tGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAAttt                                                                                                                                <  2:85663/150‑1 (MQ=255)
                     cGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTtat                                                                                                                              >  1:754477/1‑150 (MQ=255)
                     cGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTtat                                                                                                                              <  1:369053/150‑1 (MQ=255)
                     cGGCAGAATGGTAGTCATCATAAAAAATCAGCAGGCCGACGTACAACGAGGTGTTAGAGAGAATCACCCAGGAAACTTCATTGCTAAAAAAGCATTCATTTTCAGTTAGATGTGACATGGCAGGCCGTGGCTGCCCCGCGGCCAATTtat                                                                                                                              >  2:559638/1‑150 (MQ=255)
                      ggCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTtata                                                                                                                             <  1:350649/150‑1 (MQ=255)
                      ggCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTtata                                                                                                                             <  2:334092/150‑1 (MQ=255)
                      ggCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTtata                                                                                                                             <  1:758297/150‑1 (MQ=255)
                       gCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGGGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAg                                                                                                                            <  1:1035160/150‑1 (MQ=255)
                        cAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGc                                                                                                                           >  1:594342/1‑150 (MQ=255)
                           aaTGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGaa                                                                                                                        <  2:948600/150‑1 (MQ=255)
                               gTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCCGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGGGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTg                                                                                                                    <  2:578772/150‑1 (MQ=255)
                                   tcatcatcaGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTccc                                                                                                                >  2:1384596/1‑150 (MQ=255)
                                    catcatcaGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTcccg                                                                                                               >  2:794818/1‑149 (MQ=255)
                                    catcatcaGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCt                                                                                                               >  2:866174/1‑150 (MQ=255)
                                    catcatcaGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGATGCCCCGCGGCCAATTTATAGCGAACGTGTcccg                                                                                                               >  1:907450/1‑149 (MQ=255)
                                       catcaGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGc                                                                                                            >  1:1277320/1‑150 (MQ=255)
                                        atcaGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCg                                                                                                           >  1:728785/1‑150 (MQ=255)
                                        atcaGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCACCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCg                                                                                                           >  2:147390/1‑150 (MQ=255)
                                          caGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAg                                                                                                         >  1:1137393/1‑150 (MQ=255)
                                          caGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAg                                                                                                         >  2:651737/1‑150 (MQ=255)
                                           aGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGt                                                                                                        >  2:1363401/1‑150 (MQ=255)
                                            gAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTc                                                                                                       >  1:802528/1‑150 (MQ=255)
                                                      ggCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCg                                                                                             >  2:862775/1‑150 (MQ=255)
                                                      ggCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCg                                                                                             >  2:876531/1‑150 (MQ=255)
                                                      ggCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCg                                                                                             >  1:316242/1‑150 (MQ=255)
                                                      ggCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCg                                                                                             >  1:81649/1‑150 (MQ=255)
                                                         cgccgTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCg                                                                                          <  2:742906/150‑1 (MQ=255)
                                                           ccgTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTCGATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGAt                                                                                        <  2:402995/150‑1 (MQ=255)
                                                            cgTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGAtg                                                                                       <  2:181597/150‑1 (MQ=255)
                                                             gTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGAtgc                                                                                      <  2:540405/150‑1 (MQ=255)
                                                              tACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGAtgct                                                                                     >  2:866337/1‑150 (MQ=255)
                                                              tACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGAtgcg                                                                                     >  2:1476346/1‑149 (MQ=255)
                                                                cAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGAtgctgc                                                                                   >  1:1420997/1‑150 (MQ=255)
                                                                 aGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGcc                                                                                  <  2:728785/150‑1 (MQ=255)
                                                                  gCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCg                                                                                 <  1:1593733/150‑1 (MQ=255)
                                                                  gCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCg                                                                                 >  2:619573/1‑150 (MQ=255)
                                                                  gCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCg                                                                                 <  1:907750/150‑1 (MQ=255)
                                                                        tgtTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCt                                                                           <  2:1420997/150‑1 (MQ=255)
                                                                         gtTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACACGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTg                                                                          >  1:500541/1‑150 (MQ=255)
                                                                           ttGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGa                                                                        >  2:873362/1‑150 (MQ=255)
                                                                                  cATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGtt                                                                 >  2:760730/1‑150 (MQ=255)
                                                                                    tCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCa                                                               <  1:866174/150‑1 (MQ=255)
                                                                                      ccccAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCAcc                                                             >  2:992227/1‑150 (MQ=255)
                                                                                             aaaCTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGAtgct                                                                                     >  2:956667/1‑119 (MQ=255)
                                                                                             aaaCTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGAtgct                                                                                     <  1:956667/119‑1 (MQ=255)
                                                                                               aCTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCAtct                                                    <  1:108042/150‑1 (MQ=255)
                                                                                                cTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCAtctt                                                   <  1:1431917/150‑1 (MQ=255)
                                                                                                   taTTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGAAATGGAAGGACGTGGATGACACGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATcttctt                                                >  2:1332581/2‑150 (MQ=255)
                                                                                                       gCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGt                                            <  1:904563/150‑1 (MQ=255)
                                                                                                          aaaaaaGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCAc                                         >  1:1415667/1‑150 (MQ=255)
                                                                                                              aaGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTg                                     <  1:760730/150‑1 (MQ=255)
                                                                                                               aGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGa                                    <  2:802528/150‑1 (MQ=255)
                                                                                                                  tttCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACAt                                 >  2:1300244/1‑150 (MQ=255)
                                                                                                                   ttCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGTGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATc                                >  2:1104872/1‑150 (MQ=255)
                                                                                                                      aTTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAg                             <  1:1363401/150‑1 (MQ=255)
                                                                                                                        tttCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCg                           <  2:370058/150‑1 (MQ=255)
                                                                                                                             gTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGccc                      <  2:316242/150‑1 (MQ=255)
                                                                                                                              tttGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCt                     <  1:619573/150‑1 (MQ=255)
                                                                                                                              tttGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCt                     <  1:873362/150‑1 (MQ=255)
                                                                                                                                  aTGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAAc                 >  2:1489702/1‑150 (MQ=255)
                                                                                                                                  aTGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAAc                 <  2:1415667/150‑1 (MQ=255)
                                                                                                                                   tgtgACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTAACTTTGACATCCAGCGTGCCCTGAACa                >  2:1159623/1‑150 (MQ=255)
                                                                                                                                     tgACATGGAAGGCCGTGGCTGACGCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGa              <  1:558612/150‑1 (MQ=255)
                                                                                                                                     tgACACGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGa              <  2:500541/150‑1 (MQ=255)
                                                                                                                                      gACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGAc             >  1:79318/1‑150 (MQ=255)
                                                                                                                                      gACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGAc             <  2:1295984/150‑1 (MQ=255)
                                                                                                                                       aCATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACt            >  1:552082/1‑150 (MQ=255)
                                                                                                                                        cATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGGCCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTg           >  1:1045734/1‑150 (MQ=255)
                                                                                                                                          tGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGAt         >  1:974094/1‑150 (MQ=255)
                                                                                                                                          tGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGAt         >  2:303933/1‑150 (MQ=255)
                                                                                                                                           ggAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGAtt        >  2:479400/1‑150 (MQ=255)
                                                                                                                                            gAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTg       <  2:1137393/150‑1 (MQ=255)
                                                                                                                                             aaGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGa      <  1:1013405/150‑1 (MQ=255)
                                                                                                                                              tggCCGCGGGTCAGCCACGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAg     >  1:1052822/2‑150 (MQ=255)
                                                                                                                                               ggCCGCGGGTCAGCCACGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAgg    >  1:1255956/1‑150 (MQ=255)
                                                                                                                                               ggCCGCGGGTCAGCCACGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAgg    >  1:1062550/1‑150 (MQ=255)
                                                                                                                                               ggCAGTGGCTGACCCGCGTCCAATTTATAGCGAACCTGTGCCTTTAGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTTTAATTTGTTCACACTACCATCTTCTTCCGTCAATTTTACATTCAGCGTTCCCTCAACAGACTGATTTAgg    <  2:320905/150‑1 (MQ=255)
                                                                                                                                                 ccGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGtc  <  2:594342/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTC  >  NC_000913/750969‑751263

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