Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 156,683 G→A 5.6% G67G (GGC→GGT yadN ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913156,6830GA 5.6% 285.0 / 3.7 88G67G (GGC→GGTyadNputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (40/43);  new base A (3/2);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTG  >  NC_000913/156535‑156831
                                                                                                                                                    |                                                                                                                                                    
tGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACg                                                                                                                                                      <  8:346933/149‑1 (MQ=255)
tGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACg                                                                                                                                                      >  1:397192/1‑149 (MQ=255)
           aacgTCCCTTTGCTGTTACCACAGAAAACAGACCCGAATGTCATTTTAGCCGTGCTACCTGGATTCGGATTTTCTTAGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGAGCCAACGGTGTCATTCAGTACGCCAGCCTCGAt                                                                                                                                           <  1:113261/146‑1 (MQ=255)
           aaaGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGAt                                                                                                                                           >  2:427106/1‑149 (MQ=255)
            aaGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGAtt                                                                                                                                          <  7:64146/149‑1 (MQ=255)
               tCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCa                                                                                                                                       <  4:118218/149‑1 (MQ=255)
                 ccTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCAc                                                                                                                                      <  6:70360/148‑1 (MQ=255)
                      gCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCGGAAAGGTTTAGCTTTAGCGCCAACGGGGTCATTCAGTACGCCAGCATCGATTTCACCAACg                                                                                                                                 <  6:195534/148‑1 (MQ=255)
                          ttACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGc                                                                                                                            <  1:41220/149‑1 (MQ=255)
                            aCCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCgg                                                                                                                          >  3:258569/1‑149 (MQ=255)
                             ccAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGt                                                                                                                         <  7:324967/149‑1 (MQ=255)
                             ccAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGt                                                                                                                         >  2:225553/1‑149 (MQ=255)
                              cAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCACCGGTTGCGGtt                                                                                                                        >  5:244269/1‑149 (MQ=255)
                               aaaGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGttt                                                                                                                       >  4:265601/1‑149 (MQ=255)
                                 aGAAAACAGAACGGAAGGTCATTTTAGCGGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTTATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTgc                                                                                                                     <  8:170203/149‑1 (MQ=255)
                                 aGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTTCCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTgc                                                                                                                     <  1:430176/149‑1 (MQ=255)
                                 aGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTg                                                                                                                      >  5:412222/1‑148 (MQ=255)
                                   aaaaCAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTgcag                                                                                                                   >  5:236392/1‑149 (MQ=255)
                                      aCAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCAGCGATTTCACCAACGGTTGCGGTTTGcagcag                                                                                                                >  1:402774/1‑149 (MQ=255)
                                       cAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGcagcag                                                                                                                >  6:25306/1‑148 (MQ=255)
                                         gAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAggct                                                                                                             >  2:304774/1‑147 (MQ=255)
                                           aCCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCa                                                                                                           <  7:404079/149‑1 (MQ=255)
                                            ccGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGAt                                                                                                             >  7:234967/1‑146 (MQ=255)
                                            ccGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGAt                                                                                                             <  8:234967/146‑1 (MQ=255)
                                             cGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCCTTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGa                                                                                                         >  4:206087/1‑149 (MQ=255)
                                             cGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGa                                                                                                         <  4:159573/149‑1 (MQ=255)
                                             cGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGa                                                                                                         <  3:7188/149‑1 (MQ=255)
                                                  gTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCGGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATc                                                                                                    <  4:319721/149‑1 (MQ=255)
                                                  gTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGTTTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATc                                                                                                    <  2:397192/149‑1 (MQ=255)
                                                  gTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATc                                                                                                    >  6:256473/1‑149 (MQ=255)
                                                     aTTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGCCCATCCtt                                                                                                 >  8:419274/1‑149 (MQ=255)
                                                      ttttCGCCGTGCTACCTGGATTCGGATTTGCTTTGATGCAGGCATCGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGGCATTCAGTACGCCCGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTg                                                                                                <  1:119670/149‑1 (MQ=255)
                                                         tAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGtt                                                                                              <  5:298814/148‑1 (MQ=255)
                                                          aGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTa                                                                                             <  8:94397/148‑1 (MQ=255)
                                                          aGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTa                                                                                             <  2:128813/148‑1 (MQ=255)
                                                              gTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTAccgcc                                                                                        >  5:240142/1‑149 (MQ=255)
                                                              gTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTAccgcc                                                                                        <  1:225553/149‑1 (MQ=255)
                                                                   aCCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTAGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCaa                                                                                   >  5:125675/1‑149 (MQ=255)
                                                                      tGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAAcgc                                                                                >  4:51768/1‑149 (MQ=255)
                                                                      tGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAAcgc                                                                                >  8:135321/1‑149 (MQ=255)
                                                                            cGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGATTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATAGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGtt                                                                           <  3:330616/148‑1 (MQ=255)
                                                                              gATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTca                                                                        <  5:25306/149‑1 (MQ=255)
                                                                               aTTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTcac                                                                       <  8:248938/149‑1 (MQ=255)
                                                                                    ctttgctGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCg                                                                                                                           <  8:191369/92‑1 (MQ=255)
                                                                                    ctttgctGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCg                                                                                                                           >  7:191369/1‑92 (MQ=255)
                                                                                     tttgctGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGgtgt                                                                 >  6:29603/1‑149 (MQ=255)
                                                                                         ctgcAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATc                                                             <  6:240142/149‑1 (MQ=255)
                                                                                         ctgcAGTAAACGGTGAAGCTGAAAGGTTTATCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCCCCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTAt                                                              <  4:50287/148‑1 (MQ=255)
                                                                                            cAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGGAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATcaa                                                           <  7:297898/148‑1 (MQ=255)
                                                                                                 aaCGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCaacaac                                                       >  1:177351/1‑147 (MQ=255)
                                                                                                 aaCGGTGATGATGAAAGGTTTAGCTTTAGCGACAACGGTGTCATTCAGTACACCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACAGCCGTCAACGCGGGTTTAACAGGTGTTATAaacaac                                                       >  4:250711/1‑147 (MQ=255)
                                                                                                   cGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCCGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGa                                                   <  7:408647/149‑1 (MQ=255)
                                                                                                    ggTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGAc                                                  >  6:53319/1‑149 (MQ=255)
                                                                                                        aTGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACt                                              <  1:427106/149‑1 (MQ=255)
                                                                                                          gCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGa                                            >  3:143435/1‑149 (MQ=255)
                                                                                                          gCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCAATCAGTAAGCAAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGCATCATACCATCCTTGTTACCGCCGTCAACGCCGGTTTCACAGGTGTTATCAACAACCAGACCACTg                                             >  5:225960/1‑148 (MQ=255)
                                                                                                            tGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTGCGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGata                                          >  4:87804/1‑149 (MQ=255)
                                                                                                             gAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGgtgt                                                                 >  1:314574/1‑125 (MQ=255)
                                                                                                             gAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGgtgt                                                                 <  2:314574/125‑1 (MQ=255)
                                                                                                             gAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACAGCCGTCCACGCGGGTTTCACAGGTGTTATCAACAACCAGCCCACTGatat                                         >  3:126648/1‑149 (MQ=255)
                                                                                                              aaaGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATcaacaa                                                        >  8:297792/1‑133 (MQ=255)
                                                                                                              aaaGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATcaacaa                                                        <  7:297792/133‑1 (MQ=255)
                                                                                                              aaaGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGat                                           >  8:133828/1‑146 (MQ=255)
                                                                                                              aaaGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGat                                           <  7:133828/146‑1 (MQ=255)
                                                                                                               aaGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACACCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATAttt                                       <  3:250711/149‑1 (MQ=255)
                                                                                                                   tttagctttagcGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACt                                   >  8:411902/1‑149 (MQ=255)
                                                                                                                         tttagcGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccgc                             >  4:148324/1‑149 (MQ=255)
                                                                                                                           tagcGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGAccgcc                            <  6:125675/148‑1 (MQ=255)
                                                                                                                             gcgcCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGt                          <  6:412222/148‑1 (MQ=255)
                                                                                                                             gcgcCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTc                         >  8:185357/1‑149 (MQ=255)
                                                                                                                              cgcCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCACCGGTTGCGGTTTGCCGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTcc                        <  2:402774/149‑1 (MQ=255)
                                                                                                                              cgcCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTcc                        <  7:6709/149‑1 (MQ=255)
                                                                                                                                   aCGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGCCCACTGATATTTAACTGACCGCCGTCCatat                    >  4:209795/1‑148 (MQ=255)
                                                                                                                                     ggTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCCTCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACCACCAGCCCACTGATATTTAACTGCCCGCCGTCCATATcag                 >  8:319643/1‑149 (MQ=255)
                                                                                                                                     ggTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATcag                 <  3:148324/149‑1 (MQ=255)
                                                                                                                                       tgtgATTCAGTACGCGAGCATAGATTTCACCAACGGTTGCGGTTTGCAGCTGGATCAGTCCATCCTTGTTAACGCCGTGAAGGCGGGTTTCACAGGTGTTATCAACTACCTGACCACTGATATTTAACTGACCGCCGTCCATATcagca               <  6:225960/149‑1 (MQ=255)
                                                                                                                                       tgtCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATcagca               <  4:143435/149‑1 (MQ=255)
                                                                                                                                       tgtCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATcagca               <  3:265601/149‑1 (MQ=255)
                                                                                                                                        gtCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATcag                 >  3:177250/1‑146 (MQ=255)
                                                                                                                                         gcATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGAAATTTAACTGACCGCCGTCCATATcagcag              >  7:332160/2‑148 (MQ=255)
                                                                                                                                           aTTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGaa            <  4:392197/148‑1 (MQ=255)
                                                                                                                                           aTTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAg           >  2:169757/1‑149 (MQ=255)
                                                                                                                                           aTTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAg           >  1:17470/1‑149 (MQ=255)
                                                                                                                                             tCAGTACACCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAg           >  8:171486/1‑147 (MQ=255)
                                                                                                                                              cAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTATCTGACCGCCGTCCATATCAGCAGAAGCa         <  7:419274/148‑1 (MQ=255)
                                                                                                                                                gTACACCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGt                                                                                      >  2:257279/1‑69 (MQ=255)
                                                                                                                                                gTACACCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGt                                                                                      <  1:257279/69‑1 (MQ=255)
                                                                                                                                                 tACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACAACTTATATTTAACTGACCGCCTTCCATATAAGCAGAAGCAGTAc     >  2:212864/1‑149 (MQ=255)
                                                                                                                                                   cGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTAc                                                                                            <  3:224907/60‑1 (MQ=255)
                                                                                                                                                   cGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTAc                                                                                            >  4:224907/1‑60 (MQ=255)
                                                                                                                                                    gCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTg  >  6:187302/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTAGCGCCAACGGTGTCATTCAGTACGCCAGCATCGATTTCACCAACGGTTGCGGTTTGCAGCAGGATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTG  >  NC_000913/156535‑156831

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: