Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 778,380 T→A 5.1% L214Q (CTG→CAG)  tolB → periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913778,3800TA 5.1% 246.8 / 3.4 78L214Q (CTG→CAG) tolBperiplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base A (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

CTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTGGCATTCGCCTTGTC  >  NC_000913/778233‑778527
                                                                                                                                                   |                                                                                                                                                   
cTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTg                                                                                                                                                    <  5:47893/149‑1 (MQ=255)
   gTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTgg                                                                                                                                                   <  1:236878/147‑1 (MQ=255)
         ccTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTAtgtg                                                                                                                                           <  1:193071/149‑1 (MQ=255)
         ccTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTAtgt                                                                                                                                            >  7:106773/1‑148 (MQ=255)
         ccTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTAtgt                                                                                                                                            <  1:362089/148‑1 (MQ=255)
         ccTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTAtgt                                                                                                                                            >  5:244957/1‑148 (MQ=255)
           tACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGAc                                                                                                                                         >  8:280834/1‑149 (MQ=255)
             cGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGAc                                                                                                                                         <  6:244957/147‑1 (MQ=255)
                     aGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAg                                                                                                                                <  5:190556/148‑1 (MQ=255)
                     aGACCAACGGCGGTCAGTTCCCGTATCAACTGCGCGTATCTGACAATGACGGTAACAACCAGTTTGTCGTTAACCGATCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGCTTCTAAACTGGCTTATGTGACCTTCGaaa                                                                                                                                 >  5:138050/1‑147 (MQ=255)
                      gACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGACACCAGACGGTTCTAAACTGGc                                                                                                                                                  <  4:190071/129‑1 (MQ=255)
                      gACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGACACCAGACGGTTCTAAACTGGc                                                                                                                                                  >  3:190071/1‑129 (MQ=255)
                      gACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTGCTCCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCg                                                                                                                              <  1:170957/149‑1 (MQ=255)
                              gcggTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGATTGTCGTTCACAGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTcc                                                                                                                      >  4:77359/1‑149 (MQ=255)
                                     gTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTgg                                                                                                               <  3:310895/149‑1 (MQ=255)
                                         ccGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTAt                                                                                                           >  4:300609/1‑149 (MQ=255)
                                           gTATGAACTGCGCGTCTCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTc                                                                                                         <  4:5522/149‑1 (MQ=255)
                                           gTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGt                                                                                                                            >  4:321581/1‑130 (MQ=255)
                                           gTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGt                                                                                                                            <  3:321581/130‑1 (MQ=255)
                                            tATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTAtt                                                                                                          >  4:94637/1‑147 (MQ=255)
                                              tGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGa                                                                                                      <  3:94637/149‑1 (MQ=255)
                                                        gTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCaaa                                                                                            <  7:280834/149‑1 (MQ=255)
                                                         tATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAAt                                                                                           >  8:329973/1‑149 (MQ=255)
                                                            cTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGgc                                                                                        >  6:324664/1‑149 (MQ=255)
                                                               aCTATGACGGTTACAACCAGTGTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCTTTGTCACCAGACGGTTCTAAACTTGCTTATGTGACATTCGAAAGCGGTCGTTCCTCGCTGGTTATTCATACCCTGGCAAATGGCGCt                                                                                     <  6:260452/149‑1 (MQ=255)
                                                                   tGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTg                                                                                    >  3:145796/1‑146 (MQ=255)
                                                                   tGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTg                                                                                    <  4:145796/146‑1 (MQ=255)
                                                                       ggTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCa                                                                             <  4:181781/149‑1 (MQ=255)
                                                                         ttACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAgg                                                                           >  3:304452/1‑149 (MQ=255)
                                                                         ttACAACCAGTTTGTCGTTCACAGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTAATGTGACCTTCGAAAGCGGTCGTTCCGAGCTGGTTATTCAGACGCTGGCAAATGGCGATGTACGTCAgt                                                                           >  6:47803/1‑148 (MQ=255)
                                                                             aaCCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCa                                                                             <  6:297464/143‑1 (MQ=255)
                                                                             aaCCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCa                                                                             >  5:297464/1‑143 (MQ=255)
                                                                             aaCCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGc                                                                       >  8:211617/1‑149 (MQ=255)
                                                                             aaCCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGc                                                                       <  5:324664/149‑1 (MQ=255)
                                                                             aaCCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGc                                                                       >  7:249295/1‑149 (MQ=255)
                                                                                cAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCTAATGGCGCTGTACGTCAGGTGGCTTc                                                                    <  8:106773/149‑1 (MQ=255)
                                                                                 aGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCa                                                                   <  7:329973/149‑1 (MQ=255)
                                                                                         gttcaccgttcaccGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCcgc                                                            <  4:304452/148‑1 (MQ=255)
                                                                                             accgttcaccgCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGc                                                                       >  2:86429/1‑133 (MQ=255)
                                                                                             accgttcaccgCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGc                                                                       <  1:86429/133‑1 (MQ=255)
                                                                                             accgttcaccgCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTca                                                        >  8:219494/1‑148 (MQ=255)
                                                                                               cgttcacgtCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGGCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACaa                                                     <  5:9011/149‑1 (MQ=255)
                                                                                                 ttcaccgGAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTATACTGGTTTTTGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTTTTCAGACGCTTGCAATTGGCGCTGTTCTTCAGGTGGCTTCATTCCCGCGTCACAACg                                                   <  6:138050/149‑1 (MQ=255)
                                                                                                    accgCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGATAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTg                                                <  8:249295/149‑1 (MQ=255)
                                                                                                      cgCAGCCGCTGATGTCTCCGGCGTGGTCACAAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCATACGCTGGCAAATGGCTCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCa                                              >  8:340710/1‑149 (MQ=255)
                                                                                                             gCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCAt                                       >  6:38404/1‑149 (MQ=255)
                                                                                                             gCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACAGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTc                                                                              >  4:252731/1‑110 (MQ=255)
                                                                                                             gCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACAGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTc                                                                              >  4:248168/1‑110 (MQ=255)
                                                                                                             gCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACAGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTc                                                                              <  3:252731/110‑1 (MQ=255)
                                                                                                             gCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACAGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTc                                                                              <  3:248168/110‑1 (MQ=255)
                                                                                                               tGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGt                                                                          <  1:28624/112‑1 (MQ=255)
                                                                                                               tGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGt                                                                          >  2:28624/1‑112 (MQ=255)
                                                                                                                   gTCTCCGGCGTGGTCACGAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGc                                 <  8:70858/149‑1 (MQ=255)
                                                                                                                      tcCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGc                                 >  3:271572/1‑146 (MQ=255)
                                                                                                                      tcCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAg                              >  5:60994/1‑149 (MQ=255)
                                                                                                                      tcCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGACACAACGGTGCACCTGCATTCTCGc                                 >  8:68270/1‑146 (MQ=255)
                                                                                                                        cGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCa                               <  3:219606/146‑1 (MQ=255)
                                                                                                                             tGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGgcag                       >  3:13296/1‑149 (MQ=255)
                                                                                                                             tGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGgcag                       <  7:211617/149‑1 (MQ=255)
                                                                                                                               gTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGgcagc                      >  2:106659/1‑148 (MQ=255)
                                                                                                                                  aCCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCaaa                                                                                            >  2:143727/1‑75 (MQ=255)
                                                                                                                                    cAGACTGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGGCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCATGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTg                <  5:47803/149‑1 (MQ=255)
                                                                                                                                    cAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCt                                           >  6:113490/1‑122 (MQ=255)
                                                                                                                                    cAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCt                                           <  5:113490/122‑1 (MQ=255)
                                                                                                                                    cAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGATATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACACCGGTTCACCTGCATTCTCGCCAGAAGGCAGCAAACTg                >  2:265450/1‑149 (MQ=255)
                                                                                                                                     aGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTgg               >  7:160832/1‑149 (MQ=255)
                                                                                                                                        cGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAACTGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTGGCAt            <  3:77359/149‑1 (MQ=255)
                                                                                                                                           ttCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTGGCATTCg         >  5:5406/1‑149 (MQ=255)
                                                                                                                                             cTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGTAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGgcagc                      >  4:298264/1‑134 (MQ=255)
                                                                                                                                             cTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGTAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGgcagc                      <  3:298264/134‑1 (MQ=255)
                                                                                                                                             cTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTGGCATTCGcc       >  4:152584/1‑149 (MQ=255)
                                                                                                                                              tAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTGGCATTCGCCt      >  2:27060/1‑149 (MQ=255)
                                                                                                                                              tAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTGGCATTCGCCt      >  1:327817/1‑149 (MQ=255)
                                                                                                                                                aaCTGGCTTATGTGACCTTCTAAATCGGTCGTTCCGCTCTGGTTATTCAGACGATGTCAAATGGCGCTTTACGTCAGGTGGCTTCTTTCCCGCGTAACAACGGTTCACCTGCATTCTAGCCAGACGGCAGCAAACTGGCATTCGCCtt     <  6:245897/148‑1 (MQ=255)
                                                                                                                                                aaCTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTGGCATTCGcc       <  2:23861/146‑1 (MQ=255)
                                                                                                                                                aaCTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTGGCATTCGcc       >  1:23861/1‑146 (MQ=255)
                                                                                                                                                aaCTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTGGCATTCGCCt      <  4:13296/147‑1 (MQ=255)
                                                                                                                                                  cTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTGGCATTCGCCTTGTc  >  1:31796/1‑149 (MQ=255)
                                                                                                                                                   tGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTGGCATTCGCCTTGTc  >  2:118578/1‑148 (MQ=255)
                                                                                                                                                   tGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTGGCATTCGCCTTGTc  >  6:199482/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCGGTCGTTCCGCGCTGGTTATTCAGACGCTGGCAAATGGCGCTGTACGTCAGGTGGCTTCATTCCCGCGTCACAACGGTGCACCTGCATTCTCGCCAGACGGCAGCAAACTGGCATTCGCCTTGTC  >  NC_000913/778233‑778527

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