Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,901,947 T→G 29.6% intergenic (‑362/‑101) yoaE ← / → manX putative membrane protein/conserved protein/fused mannose‑specific PTS enzymes: IIA component/IIB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,901,9470TG29.6% 97.4 / 79.9 78intergenic (‑362/‑101)yoaE/manXputative membrane protein/conserved protein/fused mannose‑specific PTS enzymes: IIA component/IIB component
Reads supporting (aligned to +/- strand):  ref base T (19/35);  new base G (11/12);  total (31/47)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

TGGCAAAGATTTACCTTCCTTTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTATTGTTATAGGCACACATGGTTGGGCTGCAGAG  >  NC_000913/1901802‑1902092
                                                                                                                                                 |                                                                                                                                                 
tGGCAAAGATTTACCTTCCTTTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGtttt                                                                                                                                                >  2:827324/1‑149 (MQ=255)
  gCAAAGATTTACCTTCCTTTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATACAATAATTCTTTCGATATCTAAAATAAATCGCCAAACGCAGGGGTTTTTg                                                                                                                                              <  4:534846/149‑1 (MQ=255)
   cAAAGATTTACCTTCCTTTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTgg                                                                                                                                             >  1:785171/1‑149 (MQ=255)
    aaaGATTTACCTTCCTTTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGt                                                                                                                                            >  1:367109/1‑149 (MQ=255)
     aaGATTTACCTTCCTTTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGtt                                                                                                                                           <  2:174572/149‑1 (MQ=255)
     aaGATTTACCTTCCTTTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGtt                                                                                                                                           <  2:866165/149‑1 (MQ=255)
        aTTTACCTTCCTTTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTa                                                                                                                                        <  2:367109/149‑1 (MQ=255)
              cttccttTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCtt                                                                                                                                  >  1:797666/1‑149 (MQ=255)
              cttccttTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCtt                                                                                                                                  >  4:645187/1‑149 (MQ=255)
                 ccttTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATc                                                                                                                               >  2:682121/1‑149 (MQ=255)
                 ccttTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATc                                                                                                                               <  1:11595/149‑1 (MQ=255)
                 ccttTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTATTCTTATCGAAAGTTAAATTACGGATCTTCATCACATGAAATAATTTTTTCGATATCTAAAATAAATCGCTAAACGCAGGGGTTTTTGGTTGTAGCCCTTATc                                                                                                                               <  4:13159/149‑1 (MQ=255)
                   tttGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTg                                                                                                                             >  4:470234/1‑149 (MQ=255)
                   tttGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTg                                                                                                                             >  3:627858/1‑149 (MQ=255)
                     tGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGaa                                                                                                                           <  2:507368/149‑1 (MQ=255)
                      gCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAAt                                                                                                                          <  1:865135/149‑1 (MQ=255)
                       cAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAAtc                                                                                                                         <  4:399548/149‑1 (MQ=255)
                        aaaCGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAAtcg                                                                                                                        >  4:197153/1‑149 (MQ=255)
                        aaaCGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAAtcg                                                                                                                        >  4:355434/1‑149 (MQ=255)
                         aaCGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAAtcga                                                                                                                       <  3:380904/149‑1 (MQ=255)
                          aCGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAAtcgat                                                                                                                      <  4:807614/149‑1 (MQ=255)
                              aTGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTCCGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAAtcgattcga                                                                                                                  <  1:665753/149‑1 (MQ=255)
                              aTGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAAtcgattcga                                                                                                                  <  3:721448/149‑1 (MQ=255)
                               tgtgACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTAt                                                                                                                                >  1:171095/1‑134 (MQ=255)
                               tgtgACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTAt                                                                                                                                <  2:171095/134‑1 (MQ=255)
                                  gACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATtgt                                                                                                              >  3:331186/1‑149 (MQ=255)
                                  gACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATtgt                                                                                                              >  2:6310/1‑149 (MQ=255)
                                   aCAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGccc                                                                                                                                    >  3:362828/1‑126 (MQ=255)
                                   aCAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGccc                                                                                                                                    <  4:362828/126‑1 (MQ=255)
                                   aCAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATtgtg                                                                                                             <  2:120369/149‑1 (MQ=255)
                                     aaGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGga                                                                                                           >  4:266462/1‑149 (MQ=255)
                                       ggATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGgacg                                                                                                         >  4:78015/1‑149 (MQ=255)
                                        gATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGgacga                                                                                                        >  2:170912/1‑149 (MQ=255)
                                           atTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGAcgacga                                                                                                     <  2:460220/149‑1 (MQ=255)
                                            ttttACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGAt                                                                                                    <  1:126338/149‑1 (MQ=255)
                                            ttttACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGAt                                                                                                    <  3:371332/149‑1 (MQ=255)
                                            ttttACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGAt                                                                                                    >  3:248808/1‑149 (MQ=255)
                                                   tttCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCaaaaa                                                                                             <  1:234391/149‑1 (MQ=255)
                                                   tttCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCaaaaa                                                                                             <  4:640863/149‑1 (MQ=255)
                                                                     tCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAg                                                                           <  1:827324/149‑1 (MQ=255)
                                                                     tCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAg                                                                           <  1:6310/149‑1 (MQ=255)
                                                                     tCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAg                                                                           >  3:459002/1‑149 (MQ=255)
                                                                         aaGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGgtgt                                                                       >  2:717641/1‑149 (MQ=255)
                                                                          aGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGtt                                                                      <  1:170912/149‑1 (MQ=255)
                                                                          aGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGtt                                                                      <  1:361479/149‑1 (MQ=255)
                                                                            ttAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTaa                                                                    <  1:555064/149‑1 (MQ=255)
                                                                            ttAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTaa                                                                    <  3:645187/149‑1 (MQ=255)
                                                                             tAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAAc                                                                   <  1:124939/149‑1 (MQ=255)
                                                                               aaTTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGa                                                                 <  1:525096/149‑1 (MQ=255)
                                                                                aTTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGat                                                                <  1:652172/149‑1 (MQ=255)
                                                                                 ttACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGata                                                               >  3:328115/1‑149 (MQ=255)
                                                                                   aCGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGataat                                                             <  1:877293/149‑1 (MQ=255)
                                                                                      gATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATaaa                                                          <  1:591388/149‑1 (MQ=255)
                                                                                         cTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAagga                                                       <  2:75623/149‑1 (MQ=255)
                                                                                         cTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAagga                                                       <  3:99004/149‑1 (MQ=255)
                                                                                          ttCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAaggag                                                      <  3:266462/149‑1 (MQ=255)
                                                                                          ttCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAaggag                                                      >  3:696091/1‑149 (MQ=255)
                                                                                                acaTAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAg                                                                           <  3:563543/122‑1 (MQ=255)
                                                                                                acaTAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAg                                                                           >  4:563543/1‑122 (MQ=255)
                                                                                                 caTAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCaa                                               <  1:194589/149‑1 (MQ=255)
                                                                                                 caTAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCaa                                               <  3:568694/149‑1 (MQ=255)
                                                                                                  aTAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAg                                              >  2:784297/1‑149 (MQ=255)
                                                                                                   tAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGt                                             <  1:682121/149‑1 (MQ=255)
                                                                                                         aatTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCAt                                       <  4:627858/149‑1 (MQ=255)
                                                                                                           ttttttCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTg                                     <  1:784297/149‑1 (MQ=255)
                                                                                                              tttCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTa                                  <  4:248808/149‑1 (MQ=255)
                                                                                                                tCGATATCTAAAATAAATCGCGAAACGCAGGGGGTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTAtt                                <  3:653996/149‑1 (MQ=255)
                                                                                                                  gATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTATTGt                              >  1:590329/1‑149 (MQ=255)
                                                                                                                     atCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTATTGTtat                           <  3:643754/149‑1 (MQ=255)
                                                                                                                        tAAAATAAATCGCGAAACGCAGTGGTTTTTGGTTGTAGCCCTTATCTGAATCGCTTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTATTGTTATAgg                        >  4:410394/1‑149 (MQ=255)
                                                                                                                        tAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTATTGTTATAgg                        >  1:101338/1‑149 (MQ=255)
                                                                                                                                atCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTATTGTTATAGGcacaca                  <  3:640059/147‑1 (MQ=255)
                                                                                                                                atCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTATTGTTATAGGcacaca                  >  4:640059/1‑147 (MQ=255)
                                                                                                                                     gAACCGCAGGGGCTCTTGGTTGTAGCCCTTATCTGAATCGACTCGATTGTGGTCTACGATACAAAAATACATCTGACACGCTGAGGGGTCTACGATAACAAAGGAGCTAGCAAGTGACCATTATTATTGTTGTAGGCACACATGGTTgg           >  2:158194/1‑149 (MQ=255)
                                                                                                                                      aaaCGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTATTGTTATAGGCACACATGGTTggg          <  2:590329/149‑1 (MQ=255)
                                                                                                                                            aGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTATTGTCATAGGCACACATGGTTGGGCTGCag    >  3:470875/1‑149 (MQ=255)
                                                                                                                                             ggggTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTATTGTTATAGGCACACATGGTTGGGCTGCaga   <  3:418719/149‑1 (MQ=255)
                                                                                                                                              gggTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTATTGTTATAGGCACACATGGTTGGGCTGCagag  >  1:38527/1‑149 (MQ=255)
                                                                                                                                              gggTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTATTGTTATAGGCACACATGGTTGGGCTGCagag  >  4:214423/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TGGCAAAGATTTACCTTCCTTTGCAAACGAATGTGACAAGGATATTTTACCTTTCGAAATTTCTGCTAATCGAAAGTTAAATTACGGATCTTCATCACATAAAATAATTTTTTCGATATCTAAAATAAATCGCGAAACGCAGGGGTTTTTGGTTGTAGCCCTTATCTGAATCGATTCGATTGTGGACGACGATTCAAAAATACATCTGGCACGTTGAGGTGTTAACGATAATAAAGGAGGTAGCAAGTGACCATTGCTATTGTTATAGGCACACATGGTTGGGCTGCAGAG  >  NC_000913/1901802‑1902092

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: