Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 217,030 C→T 19.0% intergenic (‑27/+27) yaeF ← / ← proS putative lipoprotein/prolyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913217,0300CT19.0% 315.7 / 54.8 142intergenic (‑27/+27)yaeF/proSputative lipoprotein/prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (62/51);  new base T (18/9);  total (82/60)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 142 total aligned reads displayed.

CAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGC  >  NC_000913/216934‑217127
                                                                                                |                                                                                                 
cAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCt                                                                                                >  1:2286955/1‑100 (MQ=255)
  aCGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGa                                                                                              >  1:2164649/1‑100 (MQ=255)
  aCGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGa                                                                                              >  2:1319319/1‑100 (MQ=255)
  aCGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAACGGCCCtca                                                                                              >  1:1634461/1‑98 (MQ=255)
  aCGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGGTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGa                                                                                              >  1:2163492/1‑100 (MQ=255)
     gTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAAt                                                                                           <  2:1800499/100‑1 (MQ=255)
     gTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAAt                                                                                           >  2:1091743/1‑100 (MQ=255)
       aCAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAAt                                                                                           <  1:2113230/98‑1 (MQ=255)
        cAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTAAGCATTTTATACAATCCATGCAAAACAAGGGACCTAGAATTCa                                                                                        >  1:230429/1‑100 (MQ=255)
         aaGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGGAGTACGCCTTTGGTTTATCCATTTTATACAATACATGTAAAAAAAGGGCCCTGAAATTCAg                                                                                       <  1:1316301/100‑1 (MQ=255)
         aaGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAg                                                                                       <  2:1839083/100‑1 (MQ=255)
         aaGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAg                                                                                       <  2:2014751/100‑1 (MQ=255)
         aaGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAg                                                                                       <  2:2328583/100‑1 (MQ=255)
         aaGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAg                                                                                       <  1:1718843/100‑1 (MQ=255)
          aGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAgg                                                                                      >  2:2359431/1‑100 (MQ=255)
          aGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAgg                                                                                      >  1:2341234/1‑100 (MQ=255)
           gCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGa                                                                                     >  1:2077402/1‑100 (MQ=255)
           gCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGa                                                                                     >  2:782472/1‑100 (MQ=255)
              gAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGAccc                                                                                  >  1:850680/1‑100 (MQ=255)
                 aGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCttt                                                                               >  2:738040/1‑100 (MQ=255)
                  gTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTc                                                                              >  1:1843157/1‑100 (MQ=255)
                   tAAAAGAGAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCt                                                                             >  2:1524381/1‑100 (MQ=255)
                   tAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCt                                                                             >  2:2346652/1‑100 (MQ=255)
                   tAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCt                                                                             >  2:2696759/1‑100 (MQ=255)
                     aaagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTgg                                                                           >  1:1056320/1‑100 (MQ=255)
                      aagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTgg                                                                           <  1:1234627/99‑1 (MQ=255)
                      aagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGc                                                                          >  1:2591729/1‑100 (MQ=255)
                       ataaagCTGCGGAGTAACAGGCGGCAGTACGCCTTTGGTTTACCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCttttttt                                                                           <  1:2656991/96‑5 (MQ=25)
                       agaaagCTGGGGAGAATCAGGCGGCAGGACGCCTTTGGTTTATGCATTTTATACAATCCAGGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttac                                                                         <  1:1867240/100‑7 (MQ=25)
                       agaaagCTGGGGAGAAGCAGGCGGCAGTACGCCTTTGGTTTATCCATCTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGc                                                                          >  1:1675522/1‑99 (MQ=255)
                       agaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttac                                                                         <  1:2532283/100‑7 (MQ=37)
                       agaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCa                                                                         >  2:1764542/1‑100 (MQ=255)
                          aagCTGGTGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttacatg                                                                      <  1:2869217/100‑10 (MQ=25)
                          aagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttacatg                                                                      <  2:1993550/100‑10 (MQ=25)
                           agCTTGGGAGTAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttacat                                                                       <  1:2379431/98‑9 (MQ=25)
                           agCTGTGGTGAAACAGGCGGCAGTACGCCTTTGTTTTGTCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttaca                                                                        <  1:2811081/97‑8 (MQ=21)
                            gCTGGGGAGAAACAGGTGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGAGCCTTTCTGGCATCAGc                                                                    >  1:1700068/1‑100 (MQ=255)
                            gCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGc                                                                    >  1:1484705/1‑100 (MQ=255)
                             cTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGcc                                                                   >  1:2060995/1‑100 (MQ=255)
                             cTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGcc                                                                   <  2:2565272/100‑1 (MQ=255)
                                cggAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAGAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCATCCttt                                                                >  1:714322/2‑100 (MQ=255)
                                   agaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAAt                                                             <  2:1289835/100‑1 (MQ=255)
                                   agaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAAt                                                             <  1:2247282/100‑1 (MQ=255)
                                    gaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATc                                                            >  2:64492/1‑100 (MQ=255)
                                    gaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATc                                                            >  1:1512777/1‑100 (MQ=255)
                                     aaaCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATGCAATCCATGTAAAAAAAGGTCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCt                                                           <  1:1588301/100‑1 (MQ=255)
                                     aaaCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACTATCCATGTAAAAAATGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAAt                                                             <  1:1764693/98‑1 (MQ=255)
                                        cAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGtt                                                        >  2:1861936/1‑100 (MQ=255)
                                         aGGCGGCAGTACGCCTTTGTTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGTACCCTTTCTGGCATCAGCCTTTAATc                                                            <  1:620961/95‑1 (MQ=255)
                                          ggcggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTc                                                      <  2:2397200/100‑1 (MQ=255)
                                          ggcggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTc                                                      >  1:994971/1‑100 (MQ=255)
                                          ggcggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTc                                                      >  1:2543098/1‑100 (MQ=255)
                                           gcggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCa                                                     <  2:1105934/100‑1 (MQ=255)
                                            cggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAc                                                    <  2:497040/100‑1 (MQ=255)
                                                aGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGa                                                <  2:2373180/100‑1 (MQ=255)
                                                     gCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATAtt                                            >  2:2368315/1‑99 (MQ=255)
                                                      ccTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcg                                          >  1:837537/1‑100 (MQ=255)
                                                       cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcga                                         >  2:1051054/1‑100 (MQ=255)
                                                       cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcga                                         >  1:2272588/1‑100 (MQ=255)
                                                       cTTTGGTGTATACATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGACATCAGCCTTTAATCTGTTTCACCAGATATTcga                                         >  1:1526047/1‑100 (MQ=255)
                                                        tttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgac                                        <  1:806165/100‑1 (MQ=255)
                                                         ttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacg                                       <  2:2913468/100‑1 (MQ=255)
                                                         ttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacg                                       <  1:171705/100‑1 (MQ=255)
                                                            gTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATg                                    >  2:1193082/1‑100 (MQ=255)
                                                               tATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCa                                 <  2:2911586/100‑1 (MQ=255)
                                                                aTCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCAc                                <  1:2018688/100‑1 (MQ=255)
                                                                aTCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCAc                                <  2:354363/100‑1 (MQ=255)
                                                                  ccATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCa                              >  2:2023657/1‑100 (MQ=255)
                                                                  ccATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCa                              <  2:1321170/100‑1 (MQ=255)
                                                                      tttATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGTCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCt                          <  1:29422/100‑1 (MQ=255)
                                                                        tataCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATTTTCGACGATGTCACCAGTCtta                        >  1:2114286/1‑100 (MQ=255)
                                                                           aCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaatt                     <  2:2741466/100‑1 (MQ=255)
                                                                             aaTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                   >  1:574072/1‑100 (MQ=255)
                                                                             aaTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTAACCAGATATTCGACGATGTCACCAGTCttaattaa                   <  1:1162380/100‑1 (MQ=255)
                                                                              aaccAAGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAAc                  <  1:1324089/98‑1 (MQ=255)
                                                                                 cATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGt               <  2:419340/100‑1 (MQ=255)
                                                                                 atgtaaaaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGt               >  1:2265453/6‑100 (MQ=255)
                                                                                  aTGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGtt              >  1:1942117/1‑100 (MQ=255)
                                                                                   tGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGttt             >  1:677245/1‑100 (MQ=255)
                                                                                   atgccagaaaGGGTCCTGAATTTCTGGGTCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGt               <  1:2743928/91‑1 (MQ=37)
                                                                                     tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtct           >  1:2634020/1‑100 (MQ=255)
                                                                                     gccaggaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtct           >  1:2889866/7‑100 (MQ=38)
                                                                                      ccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCTCCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc          >  1:2394747/5‑100 (MQ=37)
                                                                                      ccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc          >  2:2122077/5‑100 (MQ=38)
                                                                                      aaaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc          >  2:767730/1‑100 (MQ=255)
                                                                                       aaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCg         <  1:2359416/100‑1 (MQ=255)
                                                                                       aaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCg         <  1:2696743/100‑1 (MQ=255)
                                                                                           aaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGt     >  1:1956340/1‑100 (MQ=39)
                                                                                           aaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGt     >  1:1640725/1‑100 (MQ=39)
                                                                                            aGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGtt    >  1:1429058/1‑100 (MQ=255)
                                                                                            aGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGtt    >  2:1856036/1‑100 (MQ=255)
                                                                                            aGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGtt    >  2:1523614/1‑100 (MQ=255)
                                                                                             gggTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGGTATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTg   >  2:1239326/1‑100 (MQ=255)
                                                                                             gggTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTg   >  1:1469389/1‑100 (MQ=255)
                                                                                             gggCCCTGAAATGCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATACGACGATGTCACCAGTCTTAATTAACTGGTTCTCGCCGTTg   <  1:1592453/100‑1 (MQ=255)
                                                                                              ggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTCATTAACTGTTTCTCGCCGTTGc  <  1:1693607/100‑1 (MQ=255)
                                                                                              ggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGc  <  2:2430723/100‑1 (MQ=255)
                                                                                              ggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGc  <  1:2750395/100‑1 (MQ=255)
                                                                                              ggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGc  <  1:907011/100‑1 (MQ=255)
                                                                                               gCCCTGACATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCCTAATTAACTGTTTCTCGCCGTTGc  <  1:428001/99‑1 (MQ=255)
                                                                                                |                                                                                                 
CAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGC  >  NC_000913/216934‑217127

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: