Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 751,124 C→G 5.7% G349A (GGT→GCT)  ybgQ ← putative fimbrial usher protein YbgQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913751,1240CG 5.7% 295.1 / 6.7 88G349A (GGT→GCT) ybgQputative fimbrial usher protein YbgQ
Reads supporting (aligned to +/- strand):  ref base C (44/39);  new base G (3/2);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

CCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATA  >  NC_000913/750976‑751273
                                                                                                                                                    |                                                                                                                                                     
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                                                                                                                                              >  2:108042/1‑150 (MQ=255)
        ccATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCa                                                                                                                                              >  1:715319/1‑150 (MQ=255)
            tGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAAttt                                                                                                                                          <  2:85663/150‑1 (MQ=255)
              cGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTtat                                                                                                                                        <  1:369053/150‑1 (MQ=255)
              cGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTtat                                                                                                                                        >  1:754477/1‑150 (MQ=255)
              cGGCAGAATGGTAGTCATCATAAAAAATCAGCAGGCCGACGTACAACGAGGTGTTAGAGAGAATCACCCAGGAAACTTCATTGCTAAAAAAGCATTCATTTTCAGTTAGATGTGACATGGCAGGCCGTGGCTGCCCCGCGGCCAATTtat                                                                                                                                        >  2:559638/1‑150 (MQ=255)
               ggCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTtata                                                                                                                                       <  1:758297/150‑1 (MQ=255)
               ggCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTtata                                                                                                                                       <  1:350649/150‑1 (MQ=255)
               ggCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTtata                                                                                                                                       <  2:334092/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                                                                                                       >  1:316242/1‑150 (MQ=255)
                                               ggCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCg                                                                                                       >  1:81649/1‑150 (MQ=255)
                                               ggCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCg                                                                                                       >  2:862775/1‑150 (MQ=255)
                                               ggCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCg                                                                                                       >  2:876531/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:907750/150‑1 (MQ=255)
                                                           gCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCg                                                                                           <  1:1593733/150‑1 (MQ=255)
                                                           gCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCg                                                                                           >  2:619573/1‑150 (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:873362/150‑1 (MQ=255)
                                                                                                                       tttGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCt                               <  1:619573/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                   >  2:303933/1‑150 (MQ=255)
                                                                                                                                   tGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGAt                   >  1:974094/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)
                                                                                                                                               gCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCaa       <  1:303933/150‑1 (MQ=255)
                                                                                                                                                cTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGAGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTGCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCaat      <  1:1595409/150‑1 (MQ=255)
                                                                                                                                                  gACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCaataa    >  2:461292/1‑150 (MQ=255)
                                                                                                                                                   aCCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCaataat   >  2:399065/1‑150 (MQ=255)
                                                                                                                                                    cccGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAataata  <  2:81649/150‑1 (MQ=255)
                                                                                                                                                    cccGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAataata  >  2:1338681/1‑150 (MQ=255)
                                                                                                                                                    cccGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATGGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAataata  >  1:623235/1‑150 (MQ=255)
                                                                                                                                                    cccGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGACCAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGAGCCCTTAACAGACTGATTGAGGTCATCAataata  >  2:286018/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCATTTTCAGTTTGATGTGACATGGAAGGCCGTGGCTGACCCGCGGCCAATTTATAGCGAACCTGTCCCTGGCGAGTCAGGAAGGGCGTCGATGCTGCCGAAACCTGGAAATTGTTCACCCGACCATCTTCTTCCGTCACTTTGACATCCAGCGTGCCCTGAACAGACTGATTGAGGTCATCAATAATA  >  NC_000913/750976‑751273

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: ‑