Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 247,107 G→A 49.3% 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,1070GA49.3% 9.0 / 198.9 144K132K (AAG→AAAyafLputative lipoprotein and C40 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (33/40);  new base A (37/34);  total (70/74)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 146 total aligned reads displayed.

CTGGCAGGCGCTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGC  >  NC_000913/247008‑247206
                                                                                                   |                                                                                                   
  ggCAGGCGCTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTCCCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGcc                                                                                                   <  2:1213016/100‑1 (MQ=255)
  ggCAGGCGCTGGCGGAAAACAACCACTTTCTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAAcc                                                                                                   <  1:2558931/100‑1 (MQ=255)
   gCAGGCGCTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCt                                                                                                  >  1:910826/1‑100 (MQ=255)
      ggCGCTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATTAGATTACCGCAGTGGCGATTCATCGCCTGTAGCTGCTGCTTGGCAAGCCtt                                                                                                 >  1:1448278/1‑98 (MQ=255)
      ggCGCTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTAc                                                                                               >  1:2850154/1‑100 (MQ=255)
      ggCGCTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTAc                                                                                               <  1:2777058/100‑1 (MQ=255)
       gcgcTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACg                                                                                              <  2:1429870/100‑1 (MQ=255)
        cgcTTGCGGAAAACAACCACTTTTTATGCACCGACCATCTCCATAACCTTACCTAAGTGGCGATTCACCTCCTGGAGCATCAGCTTGGCAAGCCTCACg                                                                                              <  1:3650163/99‑1 (MQ=255)
        cgcTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGt                                                                                             >  2:476638/1‑100 (MQ=255)
         gcTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTc                                                                                            <  2:702609/100‑1 (MQ=255)
         gcTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTc                                                                                            >  1:4029628/1‑100 (MQ=255)
           tGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTg                                                                                          <  1:3142871/100‑1 (MQ=255)
           tGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTg                                                                                          <  2:1328774/100‑1 (MQ=255)
            ggCGGAAAACAACCACTTTTTAAGCCGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTgg                                                                                         >  1:2005829/1‑100 (MQ=255)
             gCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTggg                                                                                        >  1:398220/1‑100 (MQ=255)
              cGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTgggg                                                                                       >  2:1016492/1‑100 (MQ=255)
               ggAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGc                                                                                      <  1:2724708/100‑1 (MQ=255)
               gatctacaacCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGc                                                                                      <  2:2185095/95‑1 (MQ=255)
                  aaaCAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGt                                                                                   >  1:3002801/1‑100 (MQ=255)
                   aacaacCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTa                                                                                  <  1:4269437/100‑1 (MQ=255)
                    acaacCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTggggagatcg                                                                                 >  1:1051965/1‑94 (MQ=255)
                    acaacCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGcagatc                                                                                 >  1:2185089/1‑95 (MQ=255)
                    acaacCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTAc                                                                                 >  1:1319236/1‑100 (MQ=255)
                     caacCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTAcg                                                                                >  2:813941/1‑100 (MQ=255)
                      aacCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTAcgc                                                                               <  2:248138/100‑1 (MQ=255)
                       acCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCGGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTAcgcg                                                                              >  2:3207765/1‑100 (MQ=255)
                         cACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGc                                                                            <  1:4204800/100‑1 (MQ=255)
                          aCTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGcc                                                                           >  1:3899622/1‑100 (MQ=255)
                           cTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCt                                                                          >  1:1757912/1‑100 (MQ=255)
                               ttAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATa                                                                      <  1:733417/100‑1 (MQ=255)
                                  agcagcGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAg                                                                   >  1:2870753/1‑100 (MQ=255)
                                   gcagcGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAgg                                                                  >  1:1421272/1‑100 (MQ=255)
                                    cagcGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGc                                                                 <  2:602181/100‑1 (MQ=255)
                                      gcGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCtt                                                               <  1:1748520/100‑1 (MQ=255)
                                      gcGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCtt                                                               <  2:2954085/100‑1 (MQ=255)
                                      gcGACCATCTGCATAACATTAACGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCATGATAAAGGCt                                                                >  1:1247779/1‑99 (MQ=255)
                                       cGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCttt                                                              <  1:1987078/100‑1 (MQ=255)
                                         aCCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGa                                                            >  1:3288603/1‑100 (MQ=255)
                                          ccATCTGCATAACATTACCGGAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGAc                                                           <  2:420493/100‑1 (MQ=255)
                                          ccATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGAc                                                           <  2:1768809/100‑1 (MQ=255)
                                          ccATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGAc                                                           <  2:588577/100‑1 (MQ=255)
                                                 cATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCg                                                    <  1:860614/100‑1 (MQ=255)
                                                   tAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCggg                                                  <  1:3303953/100‑1 (MQ=255)
                                                    aaCATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGggt                                                 <  1:2696527/100‑1 (MQ=255)
                                                     aCATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGggtt                                                >  2:2307809/1‑100 (MQ=255)
                                                      cATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGggttg                                               >  2:1513925/1‑100 (MQ=255)
                                                      cATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGggttg                                               >  1:3594705/1‑100 (MQ=255)
                                                          aCCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGttt                                           >  2:3304144/1‑100 (MQ=255)
                                                            cGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGttttt                                         <  1:4028891/100‑1 (MQ=255)
                                                            cGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGttttt                                         <  2:1681573/100‑1 (MQ=255)
                                                             gAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGtttttt                                        >  2:2626677/1‑100 (MQ=255)
                                                             gAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGtttttt                                        >  2:2242069/1‑100 (MQ=255)
                                                              aaGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTa                                       <  1:1667249/100‑1 (MQ=255)
                                                               aGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGGTTTGACGGTAGCGGGTTGGTTTTTTAt                                      >  1:998335/1‑100 (MQ=255)
                                                               aGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTa                                       <  2:2820134/99‑1 (MQ=255)
                                                                gTGGCGATTCACCGCCTGGCGCAGCAGCTTGGCAAACCTTACGTCTGGTGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTAt                                      >  1:1476276/1‑99 (MQ=255)
                                                                 tGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGc                                    <  2:2419349/100‑1 (MQ=255)
                                                                  ggCGGTTCGCGGGCTGGAGCAGCGGCTTGGCGAGCCGTAAGTGTGGGGCGGTACGCGGCCTGAGAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTa                                       <  1:1052216/96‑1 (MQ=255)
                                                                  ggCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGcc                                   <  1:286442/100‑1 (MQ=255)
                                                                  ggCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGcc                                   >  1:714510/1‑100 (MQ=255)
                                                                  ggCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGcc                                   <  2:1844088/100‑1 (MQ=255)
                                                                  ggCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGcc                                   <  1:1372257/100‑1 (MQ=255)
                                                                  ggCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGca                                   >  1:2419343/1‑99 (MQ=255)
                                                                  ggCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCAGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGcc                                   >  1:2823744/1‑100 (MQ=255)
                                                                   gCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTGTGACTGTAGCGGGTTGGTTTTTTATGCCt                                  >  1:2485988/1‑100 (MQ=255)
                                                                   gCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCt                                  >  2:2699512/1‑100 (MQ=255)
                                                                   gCGATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCt                                  >  2:2699551/1‑100 (MQ=255)
                                                                     gATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTac                                <  2:1639259/100‑1 (MQ=255)
                                                                     gATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTac                                <  1:1523386/100‑1 (MQ=255)
                                                                     gATTCACCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTac                                <  2:3465663/100‑1 (MQ=255)
                                                                         cACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCt                                  >  1:1552168/1‑94 (MQ=255)
                                                                          aCCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTAcaacaa                           <  2:1332872/100‑1 (MQ=255)
                                                                          aCCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTAcaacaa                           <  2:282871/100‑1 (MQ=255)
                                                                          aCCGCCTGGAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTAcaacaa                           >  1:3212169/1‑100 (MQ=255)
                                                                            cGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGa                         <  1:829285/100‑1 (MQ=255)
                                                                              ccTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATc                       <  1:4219748/100‑1 (MQ=255)
                                                                               cTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATcc                      <  1:589723/100‑1 (MQ=255)
                                                                                 ggAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACACCAAGATCCtt                    >  2:3971487/1‑100 (MQ=255)
                                                                                  gAGCAGCAGCTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGCCTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCtt                    >  1:564738/1‑99 (MQ=255)
                                                                                   agcagcagcTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGa                  <  2:2347983/100‑1 (MQ=255)
                                                                                    gtagcagcTCGCCAAACATTACGTGTGGGGCCGTACTCGGCGTGAGAAAGGGTTTGTCTGTCGCGGGTTGGTTTTTAATGCCGCCCACAAGATCCTTTAg                 <  1:4198812/98‑1 (MQ=255)
                                                                                    gcagcagcTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAg                 >  2:2584155/1‑100 (MQ=255)
                                                                                    gcagcagcTTGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCCACAACAAGATCCTTGAg                 <  2:3288614/100‑1 (MQ=255)
                                                                                     cagcagcTTGGCAAGCCTTACGTCTGGGGCTGTACGCGTCCTGATAAATGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAg                 <  1:808635/99‑1 (MQ=255)
                                                                                     cagcagcTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAgg                <  2:2970561/100‑1 (MQ=255)
                                                                                       gcagcTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCt              >  1:3620146/1‑100 (MQ=255)
                                                                                       gcagcTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCt              <  1:2821178/100‑1 (MQ=255)
                                                                                        cagcTTGGCAAGGCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGTTTTTTTATGCCTACAACAAGATCCTTGAGGCTa             <  1:3455429/100‑1 (MQ=255)
                                                                                        cagcTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTa             <  2:3206774/100‑1 (MQ=255)
                                                                                        cagcTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTa             >  2:2659045/1‑100 (MQ=255)
                                                                                         agcTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTaa            >  2:748506/1‑100 (MQ=255)
                                                                                         agcTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTaa            >  2:1137521/1‑100 (MQ=255)
                                                                                           cTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGc          >  2:806230/1‑100 (MQ=255)
                                                                                             tGGCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTc        >  1:2283711/1‑100 (MQ=255)
                                                                                              ggCAAACCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTcc       >  2:1284688/1‑100 (MQ=255)
                                                                                               gCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTccc      >  1:2203837/1‑100 (MQ=255)
                                                                                               gCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTccc      <  2:3934495/100‑1 (MQ=255)
                                                                                                cAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCcg     >  1:3446879/1‑100 (MQ=255)
                                                                                                cAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCcg     <  1:1117414/100‑1 (MQ=255)
                                                                                                  aGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCcgcg   <  1:528115/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CTGGCAGGCGCTGGCGGAAAACAACCACTTTTTAAGCAGCGACCATCTGCATAACATTACCGAAGTGGCGATTCACCGCCTGGAGCAGCAGCTTGGCAAGCCTTACGTCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGC  >  NC_000913/247008‑247206

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: