Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,343,454 T→G Y239S (TAT→TCT)  yfaL ← adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,343,4540TG98.5% 248.3 / ‑3.5 66Y239S (TAT→TCT) yfaLadhesin
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (26/39);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

ACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTAT  >  NC_000913/2343306‑2343601
                                                                                                                                                    |                                                                                                                                                   
ccccGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGaa                                                                                                                                                   <  1:541935/150‑1 (MQ=255)
    gccgcAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCtt                                                                                                                                               >  1:610057/1‑151 (MQ=255)
      cgcAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCg                                                                                                                                             >  2:1115373/1‑151 (MQ=255)
         agagGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCa                                                                                                                                          <  1:454382/151‑1 (MQ=255)
            ggaaggaCCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGTTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAgt                                                                                                                                        <  2:375342/150‑1 (MQ=255)
            ggaaggaCCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAgtg                                                                                                                                       >  1:1034581/1‑151 (MQ=255)
                         ccATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAgg                                                                                                                          >  1:1042683/1‑151 (MQ=255)
                          cATAGCCTGCTGCGTTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGc                                                                                                                         <  1:845847/151‑1 (MQ=255)
                          cATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGc                                                                                                                         <  2:595675/151‑1 (MQ=255)
                           aTAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCa                                                                                                                        <  1:835909/151‑1 (MQ=255)
                               ccTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAgtgt                                                                                                                    >  2:1022665/1‑151 (MQ=255)
                                  gctgcGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGtt                                                                                                                   <  2:154179/149‑1 (MQ=255)
                                     gctcttttGTTCTCAACGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAGGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTAttt                                                                                                               <  1:299262/144‑1 (MQ=255)
                                           ttGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATc                                                                                                         <  2:767121/150‑1 (MQ=255)
                                            tGTTCTCCTCCACAAACACGGCCCCCGTCCGGCTTTGGCCGTCATAAACAAAAAAAACAATAAGATAGGGGGCCGGCACGCTATTGGTAATTTACCCCCCGGCAGAACCTTCGGCCGTGGTATTTGTAAAGGCAGTGTTATTTACTATCg                                                                                                        <  2:972619/150‑1 (MQ=255)
                                            tGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCg                                                                                                        <  1:1126582/150‑1 (MQ=255)
                                                 tcATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAAc                                                                                                  >  2:739731/1‑151 (MQ=255)
                                                  cATCGACTAACACGCCTCCGTTCTGGCTTTAGCTGTCAACAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATaa                                                                                                   <  1:60663/149‑1 (MQ=255)
                                                  cATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGtata                                                                                                    <  1:95920/148‑1 (MQ=255)
                                                   aTCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCt                                                                                                >  1:1075031/1‑151 (MQ=255)
                                                   aTCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCt                                                                                                <  2:727686/151‑1 (MQ=255)
                                                     cGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGa                                                                                              <  2:100767/151‑1 (MQ=255)
                                                      gACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGa                                                                                              >  1:773921/1‑150 (MQ=255)
                                                        cTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCCCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAgg                                                                                           >  1:119781/1‑151 (MQ=255)
                                                          aaCCCGCCTCCGTTCTGCCTGTAGCTGTCATCAACACAAATATCAATAAAATAGGGAGCCGTCACGCTTTTGGTATATTTCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGAt                                                                                         <  2:31980/151‑1 (MQ=255)
                                                          aaCACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAgg                                                                                           <  2:754564/149‑1 (MQ=255)
                                                           acacGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGAtg                                                                                        >  1:270663/1‑151 (MQ=255)
                                                             acGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGAtgct                                                                                      >  1:996794/1‑151 (MQ=255)
                                                              cGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGAtgctt                                                                                     >  2:1119244/1‑151 (MQ=255)
                                                                   ccGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTg                                                                                 <  2:1101986/150‑1 (MQ=255)
                                                                       tCTGGCTGTAGCTGTAATCAACAGAAATATCAATAAGATAGGGGGCCGCCCCGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGc                                                                             <  1:166383/150‑1 (MQ=255)
                                                                          ggCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCta                                                                           <  2:496813/149‑1 (MQ=255)
                                                                          ggCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCta                                                                           <  2:432834/149‑1 (MQ=255)
                                                                           gCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCtatt                                                                         >  1:296481/1‑150 (MQ=255)
                                                                            cTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCtatta                                                                        <  1:763918/150‑1 (MQ=255)
                                                                               tAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCtattattat                                                                    >  2:360247/1‑151 (MQ=255)
                                                                                aGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATc                                                                   >  1:444265/1‑151 (MQ=255)
                                                                                   tGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGt                                                                >  1:458486/1‑151 (MQ=255)
                                                                                         caACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACAt                                                           <  1:495084/150‑1 (MQ=255)
                                                                                             aGAAATATCAATAAGATAGGGATCCGTCACGCTATTGATATAGATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGAt                                                       <  2:587244/150‑1 (MQ=255)
                                                                                                aaTATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGcc                                                   >  1:447987/1‑151 (MQ=255)
                                                                                                   atCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGcccc                                                 <  1:718318/150‑1 (MQ=255)
                                                                                                   atCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGcccc                                                 <  2:75232/150‑1 (MQ=255)
                                                                                                    tCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCcg                                                <  1:702823/150‑1 (MQ=255)
                                                                                                            aTAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCcgccatcgcc                                        >  1:430890/1‑150 (MQ=255)
                                                                                                             tAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCcgccatcgccat                                      >  1:79299/1‑151 (MQ=255)
                                                                                                             tAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCcgccatcgccat                                      >  1:267876/1‑151 (MQ=255)
                                                                                                               gggAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCgccatcgccatc                                     >  2:1115918/1‑150 (MQ=255)
                                                                                                                ggAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCAc                                   <  1:598908/151‑1 (MQ=255)
                                                                                                                 gAGCCGTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACt                                  >  2:752384/1‑151 (MQ=255)
                                                                                                                      gTCACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAAc                             >  1:1124004/1‑151 (MQ=255)
                                                                                                                        cACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACtt                           <  1:325070/151‑1 (MQ=255)
                                                                                                                        cACGCTATTGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACtt                           >  1:683104/1‑151 (MQ=255)
                                                                                                                               ttGGTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATg                    >  1:12137/1‑151 (MQ=255)
                                                                                                                                  gTATATATCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGtat                 <  1:786033/151‑1 (MQ=255)
                                                                                                                                   tatataTCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGtat                 <  1:68574/150‑1 (MQ=255)
                                                                                                                                    atatatCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTatat               <  1:670395/151‑1 (MQ=255)
                                                                                                                                     tatatCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATg              >  2:488914/1‑151 (MQ=255)
                                                                                                                                        atCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGcc            <  2:45868/150‑1 (MQ=255)
                                                                                                                                         tCGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCt           <  1:396781/150‑1 (MQ=255)
                                                                                                                                          cGCCCCGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTg          <  2:148571/150‑1 (MQ=255)
                                                                                                                                              ccGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGgttgt     <  1:190509/151‑1 (MQ=255)
                                                                                                                                              ccGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGgttgt     <  2:136417/151‑1 (MQ=255)
                                                                                                                                               cGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGgttgtt    >  2:213373/1‑151 (MQ=255)
                                                                                                                                               cGCCAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGgttgtt    <  2:132219/151‑1 (MQ=255)
                                                                                                                                                 ccAGAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTAt  <  2:775384/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATCAACAGAAATATCAATAAGATAGGGAGCCGTCACGCTATTGGTATATATCGCCCCGCCATAACCTTCGGCAGTGTTATTTGTAAAGGCAGTGTTATTTACTATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTAT  >  NC_000913/2343306‑2343601

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