Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,253,830 C→G 8.7% S708T (AGT→ACT)  ycgV ← putative adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,253,8300CG 8.7% 764.0 / 52.3 220S708T (AGT→ACT) ycgVputative adhesin
Reads supporting (aligned to +/- strand):  ref base C (107/93);  new base G (11/8);  total (118/102)
Fisher's exact test for biased strand distribution p-value = 8.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.28e-02

Only 100 of 220 total aligned reads displayed.

CGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAG  >  NC_000913/1253684‑1253977
                                                                                                                                                  |                                                                                                                                                   
     ccGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCaaa                                                                                                                                              <  4:951211/149‑1 (MQ=255)
     ccGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCaaa                                                                                                                                              <  4:1268651/149‑1 (MQ=255)
     ccGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCaaa                                                                                                                                              <  2:58390/149‑1 (MQ=255)
             tCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCa                                                                                                                                      >  1:917811/1‑149 (MQ=255)
             tCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCa                                                                                                                                      >  4:1378752/1‑149 (MQ=255)
                 cTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGa                                                                                                                                           <  2:1170687/140‑1 (MQ=255)
                        ccGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGGGTCCAGGCTTCCCCCa                                                                                                                           >  3:241393/1‑149 (MQ=255)
                        ccGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCa                                                                                                                           >  1:185643/1‑149 (MQ=255)
                           gCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATaa                                                                                                                        <  3:141264/149‑1 (MQ=255)
                               tGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                    <  2:1651338/149‑1 (MQ=255)
                               tGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                    <  4:1007047/149‑1 (MQ=255)
                                  aTGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCa                                                                                                                 >  3:789338/1‑149 (MQ=255)
                                       ttGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTACAGGCTTCCCCCATAACTgc                                                                                                                    >  3:405757/1‑141 (MQ=255)
                                          aGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGata                                                                                                         <  3:953775/149‑1 (MQ=255)
                                           gCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGatat                                                                                                        >  1:700748/1‑149 (MQ=255)
                                            ccAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATAtt                                                                                                       >  3:477458/1‑149 (MQ=255)
                                            ccAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATAtt                                                                                                       >  1:274625/1‑149 (MQ=255)
                                             cAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTa                                                                                                      >  4:412790/1‑149 (MQ=255)
                                               atatACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCAGTGGCaaactg                                                                                                                                           >  3:986503/1‑107 (MQ=255)
                                               atatACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGa                                                                                                           >  3:836329/1‑142 (MQ=255)
                                                  taCAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCa                                                                                                                                      >  1:729997/1‑112 (MQ=255)
                                                  taCAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCa                                                                                                                                      <  2:729997/112‑1 (MQ=255)
                                                  taCAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTc                                                                                                 <  3:1204727/149‑1 (MQ=255)
                                                  taCAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCACTGTAACCCATGTCAAAGCCGCTCAGTTTGCAACTGGCAAAGCAGTCCAGGCTTCCCCCATAACAGCGCAACCAGATATTACCGTc                                                                                                 >  4:1589259/1‑149 (MQ=255)
                                                   aCAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCt                                                                                                <  1:670269/149‑1 (MQ=255)
                                                     aGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCttt                                                                                              >  2:1607709/1‑149 (MQ=255)
                                                     aGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCttt                                                                                              >  2:31887/1‑149 (MQ=255)
                                                      gACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTa                                                                                             <  2:131692/149‑1 (MQ=255)
                                                           aCATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                    >  4:1350786/1‑121 (MQ=255)
                                                             aTACAACGGCATTACATCAGAGAGACGTTTATCTCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTGTTTACTCTGAt                                                                                      <  3:957492/149‑1 (MQ=255)
                                                             aTACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGAt                                                                                      >  2:986340/1‑149 (MQ=255)
                                                               aCAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGAttt                                                                                    >  1:102147/1‑149 (MQ=255)
                                                               aCAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGAttt                                                                                    >  3:328064/1‑149 (MQ=255)
                                                                 aaCGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATAtt                                                                                                       >  3:1052065/1‑128 (MQ=255)
                                                                 aaCGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATAtt                                                                                                       <  4:1052065/128‑1 (MQ=255)
                                                                  aCGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAATGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGc                                                                                 >  1:1227474/1‑149 (MQ=255)
                                                                   cGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCa                                                                                <  2:170613/149‑1 (MQ=255)
                                                                   cGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCa                                                                                <  3:1269642/149‑1 (MQ=255)
                                                                    ggCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCaaa                                                                                                                                              <  4:1342008/86‑1 (MQ=255)
                                                                    ggCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCaaa                                                                                                                                              >  3:1342008/1‑86 (MQ=255)
                                                                      cATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGt                                                                             >  3:21907/1‑149 (MQ=255)
                                                                        ttACATCAGAGAGGCGTTTATCCCCTCAAAACTGGATACCGCTGTAACCCATGTCAAAGCCGATCAGTCTGCCAATGGTAAAGTAGTCTATGCTTCCCCCACAACTGCGCAAGCAGATAGTCCCGCCCCTCCTCTTATTTCGCAGGTCa                                                                           <  2:1371642/149‑1 (MQ=255)
                                                                           cATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCAccc                                                                        <  4:328064/149‑1 (MQ=255)
                                                                            aTCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCa                                                                       >  3:134961/1‑149 (MQ=255)
                                                                            aTCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCa                                                                       >  3:1156472/1‑149 (MQ=255)
                                                                             tCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCATTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCAATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCAt                                                                      >  4:1607536/1‑149 (MQ=255)
                                                                             tCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATa                                                                                                                         <  2:811833/98‑1 (MQ=255)
                                                                             tCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATa                                                                                                                         >  1:811833/1‑98 (MQ=255)
                                                                             tCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCAt                                                                      >  3:893446/1‑149 (MQ=255)
                                                                                 agagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGt                                                                  >  2:659578/1‑149 (MQ=255)
                                                                                  gagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTAcc                                                                                                    >  1:1455200/1‑114 (MQ=255)
                                                                                  gagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACGCATCCGtt                                                                 <  3:1553348/149‑1 (MQ=255)
                                                                                  gagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGtt                                                                 <  4:1294822/149‑1 (MQ=255)
                                                                                  gagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGtt                                                                 >  1:1551047/1‑149 (MQ=255)
                                                                                   agaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                    <  3:969053/97‑1 (MQ=255)
                                                                                   agaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                    >  4:969053/1‑97 (MQ=255)
                                                                                     aCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAATGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCa                                                              >  1:1209349/1‑149 (MQ=255)
                                                                                     aCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCa                                                              >  4:605866/1‑149 (MQ=255)
                                                                                      cGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCTCTGCCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCAt                                                             <  2:1227474/149‑1 (MQ=255)
                                                                                        tttATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCa                                                           <  4:866682/149‑1 (MQ=255)
                                                                                               ccaccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACg                                                    <  1:1673552/149‑1 (MQ=255)
                                                                                               ccaccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACg                                                    >  2:461169/1‑149 (MQ=255)
                                                                                                caccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACgg                                                   >  3:1602650/1‑149 (MQ=255)
                                                                                                caccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACgg                                                   >  1:528450/1‑149 (MQ=255)
                                                                                                 accaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGt                                                  >  3:1532777/1‑149 (MQ=255)
                                                                                                     aaCTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTc                                              >  4:99384/1‑149 (MQ=255)
                                                                                                        tGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAAc                                           <  4:499120/149‑1 (MQ=255)
                                                                                                        tGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAAc                                           >  2:1491402/1‑149 (MQ=255)
                                                                                                          gATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACAt                                         >  1:74440/1‑149 (MQ=255)
                                                                                                           aTACCGCTGTAGCCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATa                                        >  3:1353091/1‑149 (MQ=255)
                                                                                                           aTACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATa                                        <  1:1147204/149‑1 (MQ=255)
                                                                                                           aTACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATa                                        <  3:1378752/149‑1 (MQ=255)
                                                                                                              ccGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCAc                                                                          >  2:1421155/1‑112 (MQ=255)
                                                                                                               cGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTc                                    <  3:761351/149‑1 (MQ=255)
                                                                                                                  tGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAAt                                 <  1:461169/149‑1 (MQ=255)
                                                                                                                   gTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATa                                        >  4:1181298/1‑141 (MQ=255)
                                                                                                                     aaCCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAg                              <  2:883600/149‑1 (MQ=255)
                                                                                                                     aaCCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAg                              <  3:1668114/149‑1 (MQ=255)
                                                                                                                        ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTaa                           <  1:986340/149‑1 (MQ=255)
                                                                                                                         cATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAg                                       >  4:767694/1‑136 (MQ=255)
                                                                                                                          aTGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAAcc                         <  4:898815/149‑1 (MQ=255)
                                                                                                                          aTGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAAcc                         <  4:169972/149‑1 (MQ=255)
                                                                                                                          aTGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAAcc                         >  3:1309086/1‑149 (MQ=255)
                                                                                                                             tCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGAc                      >  4:1493493/1‑149 (MQ=255)
                                                                                                                               aaaGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGAt                                                                                      >  1:1417267/1‑83 (MQ=255)
                                                                                                                               aaaGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAAc                          >  4:1545273/1‑142 (MQ=255)
                                                                                                                                  gCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTga                 <  3:589190/149‑1 (MQ=255)
                                                                                                                                   ccGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTgag                >  2:276723/1‑149 (MQ=255)
                                                                                                                                       ctcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTga                 <  4:1083569/140‑1 (MQ=255)
                                                                                                                                       ctcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTga                 >  3:1083569/5‑144 (MQ=255)
                                                                                                                                        tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGc                                                               <  1:495238/94‑1 (MQ=255)
                                                                                                                                        tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAg                              <  1:408492/127‑1 (MQ=255)
                                                                                                                                        tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAg                              >  2:408492/4‑130 (MQ=255)
                                                                                                                                        tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAt           >  3:1520405/4‑149 (MQ=255)
                                                                                                                                        tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAt           >  1:1640741/4‑149 (MQ=255)
                                                                                                                                        tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAt           >  2:941816/4‑149 (MQ=255)
                                                                                                                                         aGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGt                                                                                                  >  1:1303045/1‑61 (MQ=255)
                                                                                                                                         aGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAtt          <  3:289879/149‑1 (MQ=255)
                                                                                                                                                cACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACTGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCa   >  4:1363780/1‑149 (MQ=255)
                                                                                                                                                 aCTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAg  >  2:152975/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAG  >  NC_000913/1253684‑1253977

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