Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,253,248 T→G N902T (AAC→ACC)  ycgV ← putative adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,253,2480TG97.1% 235.0 / NA 68N902T (AAC→ACC) ycgVputative adhesin
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (39/27);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATATCCTCCGTTGACTTGCTTCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCCAGCAGCGATTCGTAGTGATTCAGATGTA  >  NC_000913/1253100‑1253395
                                                                                                                                                    |                                                                                                                                                   
atatCCTCCGTTGACTTGCTTCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAggtg                                                                                                                                                   <  2:643116/151‑1 (MQ=255)
    cctccGTTGACTTGCTTCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCa                                                                                                                                               <  2:348103/151‑1 (MQ=255)
        cGTTGACTTGCTTCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAAc                                                                                                                                             >  1:894071/1‑149 (MQ=255)
          ttGACTTGCTTCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAAc                                                                                                                                             >  2:228639/1‑147 (MQ=255)
           tGACTTGCTTCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGat                                                                                                                                         <  2:406704/150‑1 (MQ=255)
                      tGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATccc                                                                                                                             >  1:354978/1‑151 (MQ=255)
                        aTCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCt                                                                                                                           <  2:920914/151‑1 (MQ=255)
                         tCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTg                                                                                                                          >  2:519949/1‑151 (MQ=255)
                          cAAAGAGTTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTa                                                                                                                                 <  1:516027/143‑1 (MQ=255)
                          cAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTa                                                                                                                                 >  2:516044/1‑143 (MQ=255)
                          cAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTg                                                                                                                          >  1:644445/1‑150 (MQ=255)
                             agagGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGaaa                                                                                                                       >  1:135036/1‑150 (MQ=255)
                              cagGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGaaaa                                                                                                                      >  1:856667/2‑150 (MQ=255)
                                    aCCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTcgcg                                                                                                                <  1:766017/150‑1 (MQ=255)
                                        tGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGAt                                                                                                             >  1:195436/1‑149 (MQ=255)
                                         gCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAgc                                                                                                          <  1:351764/151‑1 (MQ=255)
                                              tAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCa                                                                                                      <  2:642402/150‑1 (MQ=255)
                                               aaTCCGCTTCGAGCTAGAATGTGTGCTGTTTGTTATACTGCGCACTGACTCCCACGCCGTTATTCCAGCCATTATCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGcc                                                                                                       >  1:611958/1‑148 (MQ=255)
                                                aTCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTc                                                                                                   >  1:99760/1‑151 (MQ=255)
                                                         gagaTAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAg                                                                                          <  2:247361/151‑1 (MQ=255)
                                                          agaTAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACAt                                                                                            >  2:648735/1‑148 (MQ=255)
                                                          agaTAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACAt                                                                                            >  2:914971/1‑148 (MQ=255)
                                                          agaTAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGa                                                                                         <  1:719433/151‑1 (MQ=255)
                                                            aTAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGAc                                                                                        >  2:578420/1‑150 (MQ=255)
                                                            aTAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGAc                                                                                        >  2:696720/1‑150 (MQ=255)
                                                            aTAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACa                                                                                       <  2:188099/151‑1 (MQ=255)
                                                             tAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACAt                                                                                      >  1:234756/1‑151 (MQ=255)
                                                               gAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTc                                                                                    >  1:695846/1‑151 (MQ=255)
                                                                 aTGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAg                                                                                  <  2:770438/151‑1 (MQ=255)
                                                                   gtgtgCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAgct                                                                                >  2:751046/1‑151 (MQ=255)
                                                                   gtgtgATGTTTGTTATACTGTGCAATGCCTCCCCCGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTc                                                                                    <  1:68088/147‑1 (MQ=255)
                                                                    tgtgCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAgct                                                                                <  1:632364/150‑1 (MQ=255)
                                                                    tgtgCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAgct                                                                                <  1:49951/150‑1 (MQ=255)
                                                                       gCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAgctggct                                                                            <  2:508584/151‑1 (MQ=255)
                                                                       gCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAgctggct                                                                            >  2:424214/1‑151 (MQ=255)
                                                                       gCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAgctgg                                                                              <  1:363035/149‑1 (MQ=255)
                                                                       gCTGTTTGTTATACTGTGCACTGACTCCAACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAgctggct                                                                            <  1:598326/151‑1 (MQ=255)
                                                                        cTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAgctggctg                                                                           >  1:205771/1‑151 (MQ=255)
                                                                         tgtttgttATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAgctg                                                                               <  1:877437/146‑1 (MQ=255)
                                                                         tgtttgttATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGt                                                                          >  2:847716/1‑151 (MQ=255)
                                                                            ttgttATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGc                                                                       >  1:703423/1‑151 (MQ=255)
                                                                               ttATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCt                                                                     >  2:514990/1‑150 (MQ=255)
                                                                                  taCTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGCGGCGTTCAACAGATATTCGGTATCCCCTGAACACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCgg                                                                 >  1:160731/1‑151 (MQ=255)
                                                                                      gtgCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGa                                                              <  2:737257/150‑1 (MQ=255)
                                                                                      gtgCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGAt                                                             >  2:617781/1‑151 (MQ=255)
                                                                                      gtgCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGAt                                                             <  2:187226/151‑1 (MQ=255)
                                                                                               cTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATAccc                                                    >  1:903374/1‑151 (MQ=255)
                                                                                                       ccGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATc                                             >  1:669511/1‑150 (MQ=255)
                                                                                                          ttATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGc                                         >  1:620054/1‑151 (MQ=255)
                                                                                                           tATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCt                                        >  2:904134/1‑151 (MQ=255)
                                                                                                                cAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCAc                                   >  2:39152/1‑151 (MQ=255)
                                                                                                                cAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCAc                                   <  2:694156/151‑1 (MQ=255)
                                                                                                                 aGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACg                                  <  1:204568/151‑1 (MQ=255)
                                                                                                                   ccATTACCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGcc                                >  1:624432/1‑151 (MQ=255)
                                                                                                                       tACCTTTGAAACTGTCCTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCca                             <  2:441140/150‑1 (MQ=255)
                                                                                                                        aCCTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCcagc                           >  1:518580/1‑151 (MQ=255)
                                                                                                                         ccTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCcagc                           >  1:828914/1‑150 (MQ=255)
                                                                                                                         ccTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCcagc                           >  1:6058/1‑150 (MQ=255)
                                                                                                                         ccTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCcag                            <  1:614964/149‑1 (MQ=255)
                                                                                                                          cTTTGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCcagc                           >  1:742240/1‑149 (MQ=255)
                                                                                                                           tttGAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCcagcag                         >  2:515722/1‑150 (MQ=255)
                                                                                                                              gAAACTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCCAGCAGCGAt                     >  1:310437/1‑151 (MQ=255)
                                                                                                                               aaaCTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCCAGCAGCGAtt                    <  2:464706/151‑1 (MQ=255)
                                                                                                                               aaaCTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCCAGCAGCGAtt                    <  1:142147/151‑1 (MQ=255)
                                                                                                                                 aCTGTACTTCTCCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCCAGCAGCGATTCg                  <  1:733571/151‑1 (MQ=255)
                                                                                                                                          ctcCCGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCCAGCAGCGATTCGTAGTGATTc         >  1:433545/1‑151 (MQ=255)
                                                                                                                                            cccGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCCAGCAGCGATTCGTAGTGAtt          <  2:378588/148‑1 (MQ=255)
                                                                                                                                            cccGGGAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCCAGCAGCGATTCGTAGTGATTc         >  1:582888/1‑149 (MQ=255)
                                                                                                                                                 gAGGTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCCAGCAGCGATTCGTAGTGATTCAGATGta  <  2:694098/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATATCCTCCGTTGACTTGCTTCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGCTGTTGCCTGCGGTGATGTCATACCCCAGTATCATGCTGGCACGCCCCAGCAGCGATTCGTAGTGATTCAGATGTA  >  NC_000913/1253100‑1253395

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