Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 428,188 G→T 100% R181L (CGC→CTC)  secD → SecYEG protein translocase auxillary subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913428,1880GT98.3% 565.0 / ‑1.5 171R181L (CGC→CTC) secDSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (72/96);  total (74/97)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 173 total aligned reads displayed.

CAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGCGGT  >  NC_000913/428091‑428284
                                                                                                 |                                                                                                
cAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCTCTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCctcg                                                                                                <  1:2381339/100‑1 (MQ=255)
   aaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATg                                                                                             <  1:3276816/100‑1 (MQ=255)
    atctCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGc                                                                                            <  1:3690456/97‑1 (MQ=255)
      atCGACAGCCTACGCAGTAACCTGCCCCAAAAGGCCATCCCATATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCt                                                                                           <  2:3436895/99‑1 (MQ=255)
        cGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTaaa                                                                                        >  2:1835968/1‑100 (MQ=255)
        cGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTaaa                                                                                        <  2:404905/100‑1 (MQ=255)
         gATAGCCGACGCAGTGACCTGCGAGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTaaa                                                                                        <  2:3816710/99‑1 (MQ=255)
          aTAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGc                                                                                      >  1:2881238/1‑100 (MQ=255)
          aTAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGc                                                                                      >  1:2504129/1‑100 (MQ=255)
           tAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCt                                                                                     >  1:1154979/1‑100 (MQ=255)
            aGCCTACGCAGCGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAACGAAAACAACTACGGCCTGAGCATCAATTTCCTCGATACTAAAGcac                                                                                    >  2:663031/1‑98 (MQ=255)
             gCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCg                                                                                   >  2:2389021/1‑100 (MQ=255)
             gCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCg                                                                                   <  1:2354148/100‑1 (MQ=255)
             gCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCg                                                                                   <  2:1688935/100‑1 (MQ=255)
              ccTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGt                                                                                  >  1:666882/1‑100 (MQ=255)
               cTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtg                                                                                 <  1:1931128/100‑1 (MQ=255)
                 aCGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtgat                                                                               >  1:3239628/1‑100 (MQ=255)
                 aCGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtgat                                                                               <  2:510226/100‑1 (MQ=255)
                    cAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGaa                                                                            >  2:1271476/1‑100 (MQ=255)
                     aGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAg                                                                           >  2:1325867/1‑100 (MQ=255)
                      gTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGc                                                                          <  1:1049073/100‑1 (MQ=255)
                      gTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGc                                                                          <  2:448412/100‑1 (MQ=255)
                         aCCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCAt                                                                       <  2:553474/100‑1 (MQ=255)
                          ccTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCAtt                                                                      >  1:232431/1‑100 (MQ=255)
                          ccTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATACTACAGCTCGGAATGAAGCCAtt                                                                      >  2:2613901/1‑100 (MQ=255)
                           cTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTg                                                                     <  1:1344388/100‑1 (MQ=255)
                           cTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTg                                                                     <  2:2969734/100‑1 (MQ=255)
                           cTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTg                                                                     <  2:3099201/100‑1 (MQ=255)
                             gcgcGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCAtt                                                                      <  2:3293775/97‑1 (MQ=255)
                             gcgcGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCg                                                                   >  2:2281039/1‑100 (MQ=255)
                               gcgAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTa                                                                 <  2:4138098/100‑1 (MQ=255)
                                cgAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTGCCTCGTTGCTAAAGATCGTGATGAAGCCATTGCGTAt                                                                >  2:2793822/1‑100 (MQ=255)
                                ccaaaaGGGCATCCCTTATACCACTGTTCCTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGt                                                                  <  2:747726/96‑1 (MQ=255)
                                 gAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATc                                                               <  2:3026018/100‑1 (MQ=255)
                                   aaaGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTg                                                             <  1:2617919/100‑1 (MQ=255)
                                   aaaGGGCATCCCGTATACCACTATTCGTAAAGAAAACAACTACAGCATGAGCATCACTTTCCTCGATGCTAATGCTCGTGAGGAAGCCATTGCGTATCTg                                                             <  2:1633478/100‑1 (MQ=255)
                                     aGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAg                                                           <  2:1557753/100‑1 (MQ=255)
                                       ggCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCa                                                         <  1:2125010/100‑1 (MQ=255)
                                       ggCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCa                                                         <  2:2847544/100‑1 (MQ=255)
                                       ggCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCa                                                         >  1:4038732/1‑100 (MQ=255)
                                       ggCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCa                                                         <  1:375418/100‑1 (MQ=255)
                                        gCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCaa                                                        >  1:2657275/1‑100 (MQ=255)
                                        gCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCaa                                                        >  1:3055458/1‑100 (MQ=255)
                                         cATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCAATTGCGTATCTGAGCAAg                                                       >  1:213610/1‑100 (MQ=255)
                                           tCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgcg                                                     <  1:417666/100‑1 (MQ=255)
                                           tCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgcg                                                     >  2:224443/1‑100 (MQ=255)
                                             ccGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGcc                                                   >  1:1827459/1‑100 (MQ=255)
                                             ccGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGcc                                                   >  1:721/1‑100 (MQ=255)
                                                 ataCCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCCCTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATcc                                               >  2:2257989/1‑100 (MQ=255)
                                                 ataCCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATcc                                               >  1:3965994/1‑100 (MQ=255)
                                                  taCCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCg                                              >  1:284151/1‑100 (MQ=255)
                                                    ccACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGa                                            >  1:98521/1‑100 (MQ=255)
                                                    ccACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGa                                            >  2:3157514/1‑100 (MQ=255)
                                                     cACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAc                                           >  1:3678383/1‑100 (MQ=255)
                                                     cACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAc                                           <  2:3657754/100‑1 (MQ=255)
                                                     cACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAc                                           <  1:1953414/100‑1 (MQ=255)
                                                     cACTGTTCGTAAAGAAAACAACGACGGCCTGAGCATCACTGTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAc                                           <  2:3941930/100‑1 (MQ=255)
                                                      aCTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAcc                                          <  1:3554421/100‑1 (MQ=255)
                                                         gTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTgg                                       <  1:1501066/100‑1 (MQ=255)
                                                          ttCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGt                                      <  2:885206/100‑1 (MQ=255)
                                                          ttCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGt                                      <  1:3929135/100‑1 (MQ=255)
                                                          ttCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGt                                      <  1:198443/100‑1 (MQ=255)
                                                           tCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTg                                     <  2:1773992/100‑1 (MQ=255)
                                                            cGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGa                                    >  2:3594870/1‑100 (MQ=255)
                                                             gTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAt                                   <  2:2542445/100‑1 (MQ=255)
                                                             gTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAt                                   <  1:2221081/100‑1 (MQ=255)
                                                             gTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAt                                   >  2:3014234/1‑100 (MQ=255)
                                                              tAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAtt                                  <  1:3270322/100‑1 (MQ=255)
                                                               aaaGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTa                                 <  2:1180146/100‑1 (MQ=255)
                                                               aaaGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTa                                 <  1:143619/100‑1 (MQ=255)
                                                                aaGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTag                                <  2:586474/100‑1 (MQ=255)
                                                                aaGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTag                                <  1:598072/100‑1 (MQ=255)
                                                                 aGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagc                               <  2:704961/100‑1 (MQ=255)
                                                                 aGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagc                               <  2:833166/100‑1 (MQ=255)
                                                                  gAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagca                              <  1:2571487/100‑1 (MQ=255)
                                                                  gAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagca                              <  2:1375169/100‑1 (MQ=255)
                                                                    aaaCAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagcagc                            >  1:441510/1‑100 (MQ=255)
                                                                    aaaCAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagcagc                            >  2:3795458/1‑100 (MQ=255)
                                                                    aaaCAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagcagc                            >  1:3419386/1‑100 (MQ=255)
                                                                    aaaCAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGACCAAGCGCCATCCGGACCTGGTGATTagcatc                            <  1:2782645/100‑3 (MQ=255)
                                                                     aacaacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGcc                           >  1:741476/1‑100 (MQ=255)
                                                                      acaacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCa                          <  2:105909/100‑1 (MQ=255)
                                                                       caacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAg                         <  1:1210572/100‑1 (MQ=255)
                                                                       caacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAg                         <  2:2863400/100‑1 (MQ=255)
                                                                       caacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGACTAGCAGCCAg                         >  2:2704796/1‑100 (MQ=255)
                                                                         acTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAggg                       >  1:1568297/1‑100 (MQ=255)
                                                                          cTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgc                      <  2:3055439/100‑1 (MQ=255)
                                                                           tACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgca                     <  1:1758616/100‑1 (MQ=255)
                                                                            aCGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcag                    <  1:1293112/100‑1 (MQ=255)
                                                                             cGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcag                    <  2:3274221/99‑1 (MQ=255)
                                                                             cGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCAGTGCGTATCTGAGCAAGCGCCACCCGGACCTGGTGATTAGCAGCCAGGgcag                    <  2:4078890/99‑1 (MQ=255)
                                                                              ggCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcagca                  <  1:557047/100‑1 (MQ=255)
                                                                               gCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcagcan                 <  2:1869707/100‑2 (MQ=255)
                                                                               gCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCaa                 >  1:1597847/1‑100 (MQ=255)
                                                                               gCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCaa                 <  1:3879039/100‑1 (MQ=255)
                                                                                ccTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAAc                <  2:1643712/100‑1 (MQ=255)
                                                                                 cTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAAcc               <  2:3467467/100‑1 (MQ=255)
                                                                                        tCACTTTCCTCGATGCTAAAGCTCGTGATGACGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcg        >  2:3487758/1‑100 (MQ=255)
                                                                                         cACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCAGCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcgt       <  1:505252/100‑1 (MQ=255)
                                                                                             ttCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGCgg   <  2:1997523/100‑1 (MQ=255)
                                                                                                 |                                                                                                
CAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGCGGT  >  NC_000913/428091‑428284

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