Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,702,835 G→A 15.5% L12L (CTG→TTG)  dppF ← dipeptide/heme ABC transporter ATPas

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,702,8350GA15.5% 472.5 / 64.5 193L12L (CTG→TTG) dppFdipeptide/heme ABC transporter ATPas
Reads supporting (aligned to +/- strand):  ref base G (75/88);  new base A (17/13);  total (92/101)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

Only 100 of 199 total aligned reads displayed.

GGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATA  >  NC_000913/3702736‑3702934
                                                                                                   |                                                                                                   
ggTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAg                                                                                                     >  1:501357/1‑100 (MQ=255)
 gTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGCTAGCCTGCAACAGc                                                                                                    >  1:3755599/1‑100 (MQ=255)
 gTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGc                                                                                                    >  1:2991342/1‑100 (MQ=255)
  ttAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCg                                                                                                   >  1:4193488/1‑100 (MQ=255)
   tAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCATGTCTATATCCTGCAACAGCgg                                                                                                  >  1:1658039/1‑100 (MQ=255)
     aacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgtt                                                                                                >  2:112768/1‑100 (MQ=255)
     aacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgtt                                                                                                >  1:4067674/1‑100 (MQ=255)
     aacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgtt                                                                                                >  2:3199148/1‑100 (MQ=255)
     aacgaaacgCCATCCAGCGCTTTAACAAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGTTCGATAGCCTGCAACAGAgg                                                                                                  >  1:2937888/1‑98 (MQ=255)
       cgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttgt                                                                                              >  2:2371968/1‑100 (MQ=255)
       caaaacgCCATCCAGCGCTTTAACCAGACGTTTCTGCTCGCACATGCCTTACTTCACCGCATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgtt                                                                                                <  1:930875/96‑1 (MQ=255)
         aaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCCAACATGCCTTTCTTCACCGGATAATGTTTTTTCCGGTCGATAGCCTGCAACAGCGGttgttg                                                                                            >  1:1559234/1‑100 (MQ=255)
          aacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGc                                                                                           <  2:184905/100‑1 (MQ=255)
          aacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGc                                                                                           <  1:4167541/100‑1 (MQ=255)
           acgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTACAACAACCGCTGTTGCa                                                                                          <  1:1180714/100‑1 (MQ=255)
            cgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCTAACATGCCTTTCTCCACCGGATAATGTTTTCTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAg                                                                                         <  1:2519487/100‑1 (MQ=255)
            cgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAg                                                                                         <  2:1245182/100‑1 (MQ=255)
             gCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATCATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAgg                                                                                        <  1:2602363/100‑1 (MQ=255)
             gCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAACCGCTGTTGCAgg                                                                                        <  1:554209/100‑1 (MQ=255)
              ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAggg                                                                                       >  1:457359/1‑100 (MQ=255)
              ccATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAggg                                                                                       <  1:1374204/100‑1 (MQ=255)
               cATCCAGCGCTTTAACCAGATGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGt                                                                                      <  2:3929605/100‑1 (MQ=255)
                aTCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCATTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTg                                                                                     >  2:4291670/1‑100 (MQ=255)
                 tCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTgg                                                                                    >  1:3389829/1‑100 (MQ=255)
                 tCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATACCCTGCAACAACCGCTGTTGCAggctat                                                                                    <  2:1189876/100‑5 (MQ=37)
                   cAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGcc                                                                                  >  2:3038430/1‑100 (MQ=255)
                   cAGCGCGTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGcc                                                                                  <  2:1779523/100‑1 (MQ=255)
                    aGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGcc                                                                                  <  1:3676950/99‑1 (MQ=255)
                      cgcTTTAACCAGACGTTCCGGCGCGGACATGCCTTTCTTCACCGGATAATGTTTTTGCAGGTCGATAGCCTGCAAGAGCGGTTGTTGCAGGGTGGCctct                                                                               <  1:1950806/100‑1 (MQ=255)
                        cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                             <  1:2075232/100‑1 (MQ=255)
                        cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                             <  2:1871713/100‑1 (MQ=255)
                        cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                             >  2:1946037/1‑100 (MQ=255)
                        cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                             >  2:2879139/1‑100 (MQ=255)
                        cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                             <  1:1725245/100‑1 (MQ=255)
                        cTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                             <  2:4003083/100‑1 (MQ=255)
                          ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTg                                                                             >  1:778354/1‑98 (MQ=255)
                          ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCg                                                                           >  1:3321519/1‑100 (MQ=255)
                          ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCg                                                                           <  2:1385896/100‑1 (MQ=255)
                          ttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCg                                                                           <  2:1192354/100‑1 (MQ=255)
                          ttAACCAGACGTTCCGGCGCGAACATGCCTTTCGTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCg                                                                           <  1:2308556/100‑1 (MQ=255)
                           tAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGt                                                                          <  1:2405520/100‑1 (MQ=255)
                            aaCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTa                                                                         >  1:3239125/1‑100 (MQ=255)
                             aCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTAc                                                                        <  2:2179372/100‑1 (MQ=255)
                             aCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTAc                                                                        >  2:1608231/1‑100 (MQ=255)
                              ccAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACt                                                                       <  2:1661797/100‑1 (MQ=255)
                                aGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCa                                                                     <  2:2601936/100‑1 (MQ=255)
                                   cGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTTTTTTTTTTTACCTGCAACAGAGGTTGTTGCAGGGTGGCCTCTTGCTTACTCATa                                                                   >  1:3069178/1‑99 (MQ=255)
                                   cGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAg                                                                  <  2:2116746/100‑1 (MQ=255)
                                   cGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAg                                                                  <  2:2655428/100‑1 (MQ=255)
                                   cGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTATTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAg                                                                  <  1:3088181/100‑1 (MQ=255)
                                    gTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAgt                                                                 <  2:659850/100‑1 (MQ=255)
                                    gTTCCGGCGCGAACATGCCTTGCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAgt                                                                 <  1:2782144/100‑1 (MQ=255)
                                     ttCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAgtg                                                                >  2:2019146/1‑100 (MQ=255)
                                      tCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAgtgt                                                               >  1:4286375/1‑100 (MQ=255)
                                        cGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCg                                                             >  2:4239760/1‑100 (MQ=255)
                                         ggCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCgg                                                            >  1:1924675/1‑100 (MQ=255)
                                          gcgcGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGc                                                           >  2:3553457/1‑100 (MQ=255)
                                           cgcgAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGcc                                                          <  2:4195416/100‑1 (MQ=255)
                                               aaCATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTccc                                                      <  1:3376794/100‑1 (MQ=255)
                                                  aTGCCTTTCTTCACCGGATAATGTTTTTTTTTTTTTTTTTCTTCCAACAGCGGTTGTTTCATTGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGc                                                   >  1:1988503/1‑100 (MQ=255)
                                                  aTGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGc                                                   <  1:1288248/100‑1 (MQ=255)
                                                   tGCCTTTCTTCACCGGATAATGTTTTTTTTTTTTTTTTTCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGc                                                   >  1:1123608/1‑99 (MQ=255)
                                                   tGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGca                                                  <  1:1532980/100‑1 (MQ=255)
                                                   tGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGca                                                  >  1:2104618/1‑100 (MQ=255)
                                                    gCCTTTCTTCACCGGATAATGTTTTTTTTGTTTGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcat                                                 >  1:3136035/1‑100 (MQ=255)
                                                     ccTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatc                                                >  2:2343023/1‑100 (MQ=255)
                                                     ccTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatc                                                <  2:1340023/100‑1 (MQ=255)
                                                      cTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcatca                                               <  2:1839274/100‑1 (MQ=255)
                                                        ttcttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCatcatc                                             <  1:1625903/100‑1 (MQ=255)
                                                          cttcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCaa                                           >  2:3859647/1‑100 (MQ=255)
                                                            tcACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGt                                         >  1:3906971/1‑100 (MQ=255)
                                                               ccGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTgtg                                      <  2:4081463/100‑3 (MQ=255)
                                                                 ggATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTa                                    <  1:2536975/100‑1 (MQ=255)
                                                                    tAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACGCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATg                                 <  1:2062010/100‑1 (MQ=255)
                                                                      aTGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGc                               <  2:3131658/100‑1 (MQ=255)
                                                                       tGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCa                              <  2:224441/100‑1 (MQ=255)
                                                                        gTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCAt                             <  2:904694/100‑1 (MQ=255)
                                                                         gtttttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCAtt                            >  1:3498529/2‑100 (MQ=255)
                                                                           ttttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTg                          <  2:1987640/100‑1 (MQ=255)
                                                                            tttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTgg                         >  2:323155/1‑100 (MQ=255)
                                                                            tttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTg                          >  2:2056200/1‑99 (MQ=255)
                                                                            tttCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTg                          <  2:1288336/99‑1 (MQ=255)
                                                                              tCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGAc                       >  1:3993069/1‑100 (MQ=255)
                                                                               cAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACt                      <  2:2042232/100‑1 (MQ=255)
                                                                               cAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACt                      <  2:100564/100‑1 (MQ=255)
                                                                                 ggTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGa                    >  1:1753268/1‑100 (MQ=255)
                                                                                  gTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGAc                   <  2:3330701/100‑1 (MQ=255)
                                                                                  gTCGATAGCCTGCAACAGCGGTGGTTGCAGGGTGGCCGCTTGCTTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGAc                   <  1:3449936/100‑1 (MQ=255)
                                                                                    ggccacccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGc                 >  2:2515322/7‑100 (MQ=38)
                                                                                       cacccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCg              >  2:889498/4‑100 (MQ=38)
                                                                                         cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc            >  2:1597732/2‑100 (MQ=255)
                                                                                         cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc            >  1:676700/2‑100 (MQ=255)
                                                                                         cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc            >  1:1796493/2‑100 (MQ=255)
                                                                                          ccTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCa           >  1:1952413/1‑100 (MQ=255)
                                                                                          ccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCa           >  1:81989/1‑100 (MQ=255)
                                                                                          ccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCa           >  1:2755516/1‑100 (MQ=255)
                                                                                           cTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAg          <  1:2088799/100‑1 (MQ=255)
                                                                                           cTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCACCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAg          >  2:734311/1‑100 (MQ=255)
                                                                                             gCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCg        <  2:3686564/100‑1 (MQ=255)
                                                                                                 cAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCa    <  1:2432399/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATA  >  NC_000913/3702736‑3702934

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: