Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 217,044 A→G 21.2% intergenic (‑41/+13) 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,0440AG21.2% 245.6 / 63.5 129intergenic (‑41/+13)yaeF/proSputative lipoprotein/prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (54/47);  new base G (18/9);  total (73/56)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 132 total aligned reads displayed.

GCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTC  >  NC_000913/216945‑217143
                                                                                                   |                                                                                                   
gCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGa                                                                                                     >  2:782472/1‑100 (MQ=255)
  cGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGAcc                                                                                                   <  1:2616628/100‑1 (MQ=255)
   gAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGAccc                                                                                                  >  1:850680/1‑100 (MQ=255)
       gTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTc                                                                                              >  1:1843157/1‑100 (MQ=255)
        tAAAAGAGAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCt                                                                                             >  2:1524381/1‑100 (MQ=255)
        tAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCttttt                                                                                             <  2:611813/100‑3 (MQ=38)
        tAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCt                                                                                             >  2:2346652/1‑100 (MQ=255)
        tAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCt                                                                                             >  2:2696759/1‑100 (MQ=255)
         aaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAGAAAAAGGGGCCTGAAATTCAGGACCCTTTCt                                                                                             >  1:2702043/1‑99 (MQ=255)
          aaagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTgg                                                                                           >  1:1056320/1‑100 (MQ=255)
           aagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTgg                                                                                           <  1:1234627/99‑1 (MQ=255)
           aagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGc                                                                                          <  1:2107177/100‑1 (MQ=255)
            ataaagCTGCGGAGTAACAGGCGGCAGTACGCCTTTGGTTTACCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCttttttt                                                                                           <  1:2656991/96‑5 (MQ=25)
            agaaagCTGGGGAGAAGCAGGCGGCAGTACGCCTTTGGTTTATCCATCTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGc                                                                                          >  1:1675522/1‑99 (MQ=255)
            agaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttac                                                                                         <  1:2532283/100‑7 (MQ=37)
            agaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGCAATTCAGGACCCTTTCTGGCa                                                                                         >  1:891125/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)
                  cTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGcc                                                                                   >  1:2060995/1‑100 (MQ=255)
                  cTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGcc                                                                                   <  2:2565272/100‑1 (MQ=255)
                     cggAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAGAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCATCCttt                                                                                >  1:714322/2‑100 (MQ=255)
                      ggAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTa                                                                               >  1:2319598/1‑100 (MQ=255)
                        agaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAAt                                                                             <  1:2382161/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)
                           aaCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                          >  1:1402030/1‑100 (MQ=255)
                             cAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGtt                                                                        >  2:1861936/1‑100 (MQ=255)
                             cAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGtt                                                                        >  2:1444021/1‑100 (MQ=255)
                              aGGCGGCAGTACGCCTTTGTTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGTACCCTTTCTGGCATCAGCCTTTAATc                                                                            <  1:620961/95‑1 (MQ=255)
                               ggcggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTc                                                                      >  2:2018699/1‑100 (MQ=255)
                               ggcggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTc                                                                      >  1:994971/1‑100 (MQ=255)
                               ggcggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTc                                                                      >  1:253841/1‑100 (MQ=255)
                               ggcggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTc                                                                      >  1:2543098/1‑100 (MQ=255)
                               ggcggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTc                                                                      <  2:2397200/100‑1 (MQ=255)
                                gcggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCa                                                                     <  2:1105934/100‑1 (MQ=255)
                                 cggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAc                                                                    <  2:497040/100‑1 (MQ=255)
                                     aGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGa                                                                <  2:2373180/100‑1 (MQ=255)
                                     aGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGa                                                                <  2:2375244/100‑1 (MQ=255)
                                      gTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGat                                                               >  2:2638986/1‑100 (MQ=255)
                                       tACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGata                                                              >  2:1627303/1‑100 (MQ=255)
                                            cTTTGGTTTATCCATTTTATACAATCCATTTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcg                                                          <  1:464296/99‑1 (MQ=255)
                                            cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcga                                                         >  1:2272588/1‑100 (MQ=255)
                                            cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcga                                                         >  2:1051054/1‑100 (MQ=255)
                                            cTTTGGTGTATACATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGACATCAGCCTTTAATCTGTTTCACCAGATATTcga                                                         >  1:1526047/1‑100 (MQ=255)
                                             tttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgac                                                        <  1:806165/100‑1 (MQ=255)
                                              ttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacg                                                       <  1:171705/100‑1 (MQ=255)
                                              ttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacg                                                       <  2:2913468/100‑1 (MQ=255)
                                               tGGTTTATCCATTTTATACAATCCAAGTAAAAAAAGGGCCCTGAAATTTAGTTGCCTTTTTGGCATCAGCATATAATCTGTTTTAACAGATATTaga                                                         >  1:1273861/1‑97 (MQ=255)
                                                 gTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATg                                                    >  2:1193082/1‑100 (MQ=255)
                                                    tATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCa                                                 <  2:2911586/100‑1 (MQ=255)
                                                    tATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCa                                                 <  1:1706247/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)
                                                         aTTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTCCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGt                                            >  2:1834016/1‑100 (MQ=255)
                                                             tataCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATTTTCGACGATGTCACCAGTCtta                                        >  1:2114286/1‑100 (MQ=255)
                                                             tataCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCtta                                        >  1:2139484/1‑100 (MQ=255)
                                                               taCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaat                                      >  1:1615892/1‑100 (MQ=255)
                                                                  aaTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                                   >  1:574072/1‑100 (MQ=255)
                                                                  aaTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTAACCAGATATTCGACGATGTCACCAGTCttaattaa                                   <  1:1162380/100‑1 (MQ=255)
                                                                   aTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAAc                                  >  2:2367854/1‑100 (MQ=255)
                                                                     ccATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTg                                >  2:2254914/1‑100 (MQ=255)
                                                                      tgatgccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCTACGATGTCACCAGTCTTAATTAACTGt                               >  2:2743946/10‑100 (MQ=25)
                                                                      atgtaaaaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGt                               >  1:2265453/6‑100 (MQ=255)
                                                                       gatgccagaaaGTGTCCTGAATTTCAGGGCCCTTTGTGGCATCCGCCTTTAATCTGTTTCACCCGATATTCGCCGATGTCCCCAGTCTTAATTAACTg                                >  1:1685415/9‑98 (MQ=21)
                                                                       aTGTAAAAAAAGAGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGtt                              >  1:1936981/1‑100 (MQ=255)
                                                                        atgccagaaaGGGTCCTGAATTTCTGGGTCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGt                               <  1:2743928/91‑1 (MQ=37)
                                                                         tgccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtc                            >  2:1255713/7‑100 (MQ=37)
                                                                          tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtct                           >  1:2634020/1‑100 (MQ=255)
                                                                          gccaggaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtct                           >  1:2889866/7‑100 (MQ=38)
                                                                           ccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCTCCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc                          >  1:2394747/5‑100 (MQ=37)
                                                                           ccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc                          >  1:1500498/5‑100 (MQ=38)
                                                                           ccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc                          >  2:2122077/5‑100 (MQ=38)
                                                                           ccagaaaGGGTCCTGAATTTCAGGGCCCGTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc                          >  1:1821260/5‑100 (MQ=38)
                                                                           aaaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc                          >  2:767730/1‑100 (MQ=255)
                                                                                aaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGt                     >  1:1956340/1‑100 (MQ=39)
                                                                                 aGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGtt                    >  2:1856036/1‑100 (MQ=255)
                                                                                 aGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGtt                    >  1:1429058/1‑100 (MQ=255)
                                                                                  gggTCCTGAATTTCAGGGCCCTTTCTGGCGTCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTGCTCGCCGTTg                   >  1:1967073/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)
                                                                                   ggCCCTGAAATTCAGGACCCTTTTTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGc                  <  2:761666/100‑1 (MQ=255)
                                                                                   ggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTCATTAACTGTTTCTCGCCGTTGc                  <  1:1693607/100‑1 (MQ=255)
                                                                                   ggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGc                  <  1:2750395/100‑1 (MQ=255)
                                                                                    gCCCTGACATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCCTAATTAACTGTTTCTCGCCGTTGc                  <  1:428001/99‑1 (MQ=255)
                                                                                      ccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcgac               <  1:1874993/100‑1 (MQ=255)
                                                                                        tGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGCCTTAATTAACTGTTTCTCCCCGTTGcgacga             >  1:1875794/1‑100 (MQ=255)
                                                                                          aaaTTCAGGACACTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCATAATTAACTGTTTCTCGCCGTTGCGACGata           <  1:780866/100‑1 (MQ=255)
                                                                                           taTTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGatat          <  1:1733977/99‑1 (MQ=255)
                                                                                            aTTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGatata         >  2:1733985/1‑99 (MQ=255)
                                                                                            aTTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATAtt         <  2:407755/100‑1 (MQ=255)
                                                                                              tCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTta       >  2:340919/1‑100 (MQ=255)
                                                                                              tCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTta       <  1:184096/100‑1 (MQ=255)
                                                                                                aGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTtata     <  1:845291/100‑1 (MQ=255)
                                                                                                   aCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTc  <  2:131833/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTC  >  NC_000913/216945‑217143

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: