Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 247,107 G→A 74.0% K132K (AAG→AAA yafL → putative lipoprotein and C40 family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913247,1070GA74.0% 253.2 / 115.0 174K132K (AAG→AAAyafLputative lipoprotein and C40 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (27/18);  new base A (62/66);  total (89/85)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.43e-01

Only 100 of 174 total aligned reads displayed.

TGGCAGGCGCTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCG  >  NC_000913/247009‑247205
                                                                                                  |                                                                                                  
tGGCAGGCGCTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGc                                                                                                   >  1:1895438/1‑100 (MQ=255)
 ggCAGGCGCTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAAcc                                                                                                  <  1:1046044/100‑1 (MQ=255)
      gcgcTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACg                                                                                             >  1:1462403/1‑100 (MQ=255)
       cgcTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGCGCAGCAGCTTGGCAAGCCTTACGt                                                                                            >  2:1652150/1‑100 (MQ=255)
        gcTGGCGGAAAACAACCACTTTTTACGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTc                                                                                           <  2:4098763/100‑1 (MQ=255)
        gcTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCTCCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTc                                                                                           >  2:1489888/1‑100 (MQ=255)
        gcTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTc                                                                                           <  1:3441350/100‑1 (MQ=255)
         cTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCt                                                                                          >  1:1269507/1‑100 (MQ=255)
           ggCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTgg                                                                                        >  1:2520187/1‑100 (MQ=255)
             cGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTGTgggg                                                                                      >  1:4076023/1‑100 (MQ=255)
             cGGAAAACAACCACCTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTgg                                                                                        >  2:1364580/1‑98 (MQ=255)
                aaaaCAACCACTTTTTAAGCAGCGACTATCTTCAGAACATAACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCgg                                                                                   <  2:2050903/100‑1 (MQ=255)
                  aacaacCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTa                                                                                 <  1:26807/100‑1 (MQ=255)
                  aacaacCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTa                                                                                 <  1:2459138/100‑1 (MQ=255)
                    caacCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTAcg                                                                               >  2:3971718/1‑100 (MQ=255)
                    caacCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTAcg                                                                               <  2:371038/100‑1 (MQ=255)
                    caacCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTAcg                                                                               >  2:2089042/1‑100 (MQ=255)
                      acCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTAcgcg                                                                             >  1:3073765/1‑100 (MQ=255)
                      acCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTAcgcg                                                                             >  2:870224/1‑100 (MQ=255)
                       ccACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCgg                                                                            >  1:500526/1‑100 (MQ=255)
                       ccACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCgg                                                                            >  2:428619/1‑100 (MQ=255)
                        cACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGc                                                                           <  2:2983657/100‑1 (MQ=255)
                          cTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCt                                                                         <  2:571033/100‑1 (MQ=255)
                              ttAAGCAGCGACCATCTGCATAATATTACCGAGGTGGCGATTCACCGCCTGGGGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATa                                                                     <  1:2862608/100‑1 (MQ=255)
                              ttAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATa                                                                     <  1:58265/100‑1 (MQ=255)
                              ttAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATa                                                                     <  1:1764917/100‑1 (MQ=255)
                              ttAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATa                                                                     <  1:1937333/100‑1 (MQ=255)
                               tAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATa                                                                     >  2:392304/1‑99 (MQ=255)
                               tAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATaa                                                                    >  1:1154913/1‑100 (MQ=255)
                                 agcagcGAGCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCGGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGTCGGTACGCGGCCTGATAAAg                                                                  <  2:2555548/100‑1 (MQ=255)
                                  gcagcGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAgg                                                                 <  1:3031584/100‑1 (MQ=255)
                                   cagcGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGc                                                                <  1:800519/100‑1 (MQ=255)
                                   cagcGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGc                                                                >  1:57382/1‑100 (MQ=255)
                                    agcGACCATCTACATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCt                                                               <  1:356046/100‑1 (MQ=255)
                                     gcGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCtt                                                              >  1:1344942/1‑100 (MQ=255)
                                       gACCATCTGCATAACATTACCGAAGTGGCGATTCACCACCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTg                                                            >  2:2223700/1‑100 (MQ=255)
                                        aCCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGTAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGa                                                           <  2:124689/100‑1 (MQ=255)
                                         ccATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGAc                                                          <  1:3608074/100‑1 (MQ=255)
                                         ccATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGAc                                                          <  2:3681547/100‑1 (MQ=255)
                                          cATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACt                                                         >  2:3191904/1‑100 (MQ=255)
                                             cTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTa                                                      >  2:1457142/1‑100 (MQ=255)
                                               gCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGc                                                    >  2:1911436/1‑100 (MQ=255)
                                               gCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTAATAAAGGCTTTGACTGTa                                                      >  2:58758/1‑98 (MQ=255)
                                                  tAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCggg                                                 <  1:3185729/100‑1 (MQ=255)
                                                    aCATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGggtt                                               <  1:2844602/100‑1 (MQ=255)
                                                      aTTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGggttgg                                             >  2:689261/1‑100 (MQ=255)
                                                       ttACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGggttggt                                            >  1:1842858/1‑100 (MQ=255)
                                                        tACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGggttggtt                                           >  2:724117/1‑100 (MQ=255)
                                                         aCCGAAGTGGCGATGCCCCGCCTGGAGCAGCAGCCTGGAAAAGCTTACGTCTGGGGAGGTACGCGGCCTGATAAACACTTTAACTTTAGCGggatgg                                             >  2:1612931/1‑97 (MQ=255)
                                                          ccGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGtttt                                         <  1:3995080/100‑1 (MQ=255)
                                                          ccGAAGTCGCGATTCACCGCCTGCAGCAGCAGCTCGGCAAACCCTACGTCCCGCGCGGTACGCGGCCTGATAAACGCTTTCCCTGCAGCGggctggttct                                         >  2:1172255/1‑98 (MQ=255)
                                                           cGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGttttt                                        <  1:3775860/100‑1 (MQ=255)
                                                           cGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGttttt                                        >  1:2146308/1‑100 (MQ=255)
                                                             aaGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCTTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTa                                      <  2:3888910/100‑1 (MQ=255)
                                                             aaGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTa                                      <  2:2453188/100‑1 (MQ=255)
                                                             aaGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTa                                      <  1:2621397/100‑1 (MQ=255)
                                                              aGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTAt                                     >  1:1118433/1‑100 (MQ=255)
                                                               gTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTAt                                     >  2:1736426/1‑99 (MQ=255)
                                                               gTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATg                                    <  1:3911579/100‑1 (MQ=255)
                                                                 ggCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGcc                                  <  1:3213684/100‑1 (MQ=255)
                                                                 ggCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGcc                                  >  2:3076599/1‑100 (MQ=255)
                                                                 ggCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGcc                                  >  1:2029754/1‑100 (MQ=255)
                                                                 ggCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGcc                                  <  2:1766714/100‑1 (MQ=255)
                                                                 ggCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGcc                                  >  1:1683726/1‑100 (MQ=255)
                                                                  gCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCt                                 <  2:1804801/100‑1 (MQ=255)
                                                                    gATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATACAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTac                               >  1:3354998/1‑100 (MQ=255)
                                                                    gATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTac                               >  1:2463660/1‑100 (MQ=255)
                                                                    gATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTac                               <  1:147784/100‑1 (MQ=255)
                                                                     aTTCACCGCCTGGAGCAGCAGCTTGTCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTATCGGGTTGGTTTTTTATGCCTaca                              <  2:2031832/100‑1 (MQ=255)
                                                                     aTTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTaca                              <  2:2901451/100‑1 (MQ=255)
                                                                        cACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTacaaca                           >  2:1338685/1‑100 (MQ=255)
                                                                        cACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTacaaca                           <  1:870229/100‑1 (MQ=255)
                                                                        cACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTacaaca                           >  1:13009/1‑100 (MQ=255)
                                                                        cACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTacaaca                           >  2:1930989/1‑100 (MQ=255)
                                                                         aCCGCCTGGCGCAGCAGCTTGGCACACCTTACGTCTGGGGCGGTACGCGGCCTGCTAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTacaaca                           >  2:2276956/1‑99 (MQ=255)
                                                                         aCCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTacaaca                           >  2:2769188/1‑99 (MQ=255)
                                                                         aCCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTAcaacaa                          <  1:1173596/100‑1 (MQ=255)
                                                                         aCCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTAcaacaa                          <  2:1107804/100‑1 (MQ=255)
                                                                          ccGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTacaaca                           <  2:1291714/98‑1 (MQ=255)
                                                                          ccGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAg                         <  2:2609996/100‑1 (MQ=255)
                                                                          ccGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAg                         <  1:2354888/100‑1 (MQ=255)
                                                                          ccGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAg                         >  2:3031817/1‑100 (MQ=255)
                                                                            gCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTCCAACAAGAt                       >  2:962526/1‑100 (MQ=255)
                                                                                  agcagcagcTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGa                 <  1:428622/100‑1 (MQ=255)
                                                                                   gcagcagcTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAg                <  2:3224443/100‑1 (MQ=255)
                                                                                     agcagcTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGc              >  2:3412854/1‑100 (MQ=255)
                                                                                      gcagcTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCt             >  2:2938828/1‑100 (MQ=255)
                                                                                        agcTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTGTTTATGCCTACAACAAGATCCTTGAGGCTa            >  2:3544097/1‑99 (MQ=255)
                                                                                         gcTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAg          <  2:1398452/100‑1 (MQ=255)
                                                                                         gcTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAg          <  1:2085061/100‑1 (MQ=255)
                                                                                         gcTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAg          >  1:1456048/1‑100 (MQ=255)
                                                                                         gcTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAg          <  2:2786922/100‑1 (MQ=255)
                                                                                          cTTGGCAAGCCTTACGTCTGGGGCGATATGCGGCCTGATAAAGTCTTTGACTGTAGCGGGTTGGTTTTTTATTCCTACAACAAGATCCTTGAGGCTAAGc         >  2:1975402/1‑100 (MQ=255)
                                                                                            tGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTc       >  1:264132/1‑100 (MQ=255)
                                                                                            tGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTc       >  1:263823/1‑100 (MQ=255)
                                                                                             ggCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTcc      >  2:2097680/1‑100 (MQ=255)
                                                                                               cAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCcg    >  2:443638/1‑100 (MQ=255)
                                                                                               cAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACATGATCCTTGAGGCTAAGCTCCcg    <  2:3850099/100‑1 (MQ=255)
                                                                                               cAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACATGATCCTTGAGGCTAAGCTCCcg    <  2:487551/100‑1 (MQ=255)
                                                                                                 aGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCcgcg  <  2:2191926/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
TGGCAGGCGCTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCG  >  NC_000913/247009‑247205

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