Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 177,405 G→T 5.8% L74M (CTG→ATG)  yadS ← UPF0126 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913177,4050GT 5.8% 372.1 / 7.3 130L74M (CTG→ATG) yadSUPF0126 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (73/47);  new base T (2/5);  total (76/54)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 134 total aligned reads displayed.

GCCCATACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCATACGCAATT  >  NC_000913/177256‑177554
                                                                                                                                                     |                                                                                                                                                     
gCCCATACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAATCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGGCTTGGCTGGCGCACCAGCACGATGGTCAg                                                                                                                                                       <  1:1926069/150‑1 (MQ=255)
   cATACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCAt                                                                                                                                                    >  1:1404612/1‑150 (MQ=255)
     tACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGc                                                                                                                                                  >  2:174633/1‑150 (MQ=255)
     tACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGc                                                                                                                                                  >  2:235677/1‑150 (MQ=255)
     tACAAACCGCGATTAACGGACAGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACCCCGCCAGACCAACGGCGTCCAGCACCGGCAAAAACCATTTTGGTAAGCTTCTTGTCTGGCGCACCAGCACGATGGTCAGCATGc                                                                                                                                                  >  2:1459905/1‑150 (MQ=255)
       cAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACTATCACGATGGTCAGCATGCTg                                                                                                                                                >  1:1211459/1‑150 (MQ=255)
       cAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTg                                                                                                                                                >  1:724217/1‑150 (MQ=255)
       cAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTg                                                                                                                                                >  2:661239/1‑150 (MQ=255)
          aCCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACTCCAATGCATACAAACACCGCCAGACCAACGGCGTCCAGCACCTGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGGCAGCATGCTGGTg                                                                                                                                             <  2:1207898/150‑1 (MQ=255)
          aCCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTg                                                                                                                                             >  1:800675/1‑150 (MQ=255)
           ccGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCAACGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGa                                                                                                                                            >  2:746250/1‑150 (MQ=255)
              cgATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCATGCTGACCATCGTGCTGGTgcgcc                                                                                                                                         <  2:1856166/150‑5 (MQ=255)
                tttAACGGATCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACTAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCATGCTGACCATCGTGCTGGTgcgccag                                                                                                                                       <  2:1262230/149‑7 (MQ=255)
                aTTAACGGACCGGCTTCCGCATTAAAGGCTTTATTAACGCCAATGCCGACAAACCCCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAATCGTCTTGGCTGGCGCACCAGCATGCTGACCATCGTGCTGGTgcgccag                                                                                                                                       <  1:1511870/150‑7 (MQ=255)
                  tAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGc                                                                                                                                     <  2:1577283/150‑1 (MQ=255)
                   aaCGGACCGGCTTCCGCATTAAAGGCGTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGGAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATCCTGGTGACCATTGcc                                                                                                                                    >  2:195484/1‑149 (MQ=255)
                     cGGACCGGCTTCCGCATTAAAGGCTTTATTCACCCCACTGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCACCAGCATGCTGACCagcgtgct                                                                                                                                  <  1:1157569/150‑8 (MQ=255)
                      ggACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAATCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACg                                                                                                                                 <  2:891726/150‑1 (MQ=255)
                          cGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCa                                                                                                                             <  1:764013/150‑1 (MQ=255)
                          cGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCa                                                                                                                             <  2:900313/150‑1 (MQ=255)
                             cTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGAt                                                                                                                          >  1:1549444/1‑150 (MQ=255)
                                 cGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGt                                                                                                                      >  1:1973382/1‑150 (MQ=255)
                                  gCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTg                                                                                                                     >  1:750892/1‑150 (MQ=255)
                                   cATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTgg                                                                                                                    >  2:608304/1‑150 (MQ=255)
                                    aTTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTggg                                                                                                                   <  2:1873450/150‑1 (MQ=255)
                                      tAAAGGCTTTATTCACGCCAATGCCGACAACCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGCCCATTGCAACGACCAGCTCGGTgggct                                                                                                                 >  2:170459/1‑148 (MQ=255)
                                      tAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCAGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGAt                                                                                                                 <  2:11023/150‑1 (MQ=255)
                                      tAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGAACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGATGGAGACCATTGCAACGACCAGATCCGTGGGAt                                                                                                                 >  2:2034867/1‑150 (MQ=255)
                                         aGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCCCCAGCACGATGGTCAGCATGCTGGTGCCCATTGCACCGACCAGATCGGTGGGATCtt                                                                                                              >  2:1340860/1‑150 (MQ=255)
                                           gCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTc                                                                                                            >  2:941693/1‑150 (MQ=255)
                                            gtttATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGGGGGATCTTTCa                                                                                                           <  2:1572628/149‑1 (MQ=255)
                                            cTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCa                                                                                                           <  1:717152/150‑1 (MQ=255)
                                            cTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGAACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGCCCAGATCGGTGGGATCTTTCa                                                                                                           >  1:231168/1‑150 (MQ=255)
                                               tATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCAccc                                                                                                        <  2:487510/150‑1 (MQ=255)
                                                aTTCACGCCAATGCCGACAAACCCCGCCAGACCAACGGCGTCCAGCACCGGCAACAGCCATTTTGGAAAGCGACTTGGCTGGCGCACCAGCACGAAGCTCAGCCTCCTGGGGACCATTGCTACGACCACATCGGTGGGAGCTTTCACCCa                                                                                                       >  2:301814/1‑150 (MQ=255)
                                                aTTCACGCCAATGCCGACAAACACCGCCAGACCACCGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCa                                                                                                       >  1:1755462/1‑150 (MQ=255)
                                                aTTAACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCAATTTGGTAAGCGTCTTGGCTGGCGCACCAGCCCGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCa                                                                                                       >  2:255341/1‑150 (MQ=255)
                                                  tcccGCCAATGCCGACAAACACCGCCAGACCACCGGCGTCCAGCACCGGCAACATCCATATTGGTAAGCGTCTTGGCTCGCGCACCAGCACGATGGTCCGCATGCTCGTGACCATAGCCACGCCCAGATCGGTGGGATCTTACACCCaaa                                                                                                     >  1:1092079/4‑150 (MQ=255)
                                                   cACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaaaa                                                                                                    >  2:740136/1‑150 (MQ=255)
                                                     cGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCAAGCTGGTGACCATTGCAACGACCAGATCGGTGGGCTCTTTCACCCAAAATa                                                                                                  >  2:1044031/1‑150 (MQ=255)
                                                          aTGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTTAGCGTCGTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCTTTGCAACGACCATATCGGTGGGAGCTTTCACCCAAAATACCggg                                                                                             <  1:241677/150‑1 (MQ=255)
                                                             ccGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCCCCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCg                                                                                          >  1:931396/1‑150 (MQ=255)
                                                               gACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTg                                                                                        >  1:1007969/1‑150 (MQ=255)
                                                                  aaaCACCGACAGACCAACCGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATc                                                                                     <  2:1830872/150‑1 (MQ=255)
                                                                   aaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCCCGATGGTCAGCATGCGGGTGACCATTGCAACGACCAGCTCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGCTcc                                                                                    >  1:1292366/1‑150 (MQ=255)
                                                                      acCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGCTCCAgc                                                                                 >  1:665181/1‑150 (MQ=255)
                                                                        cGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGAGGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAgcgc                                                                               <  1:2042574/150‑1 (MQ=255)
                                                                         gCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGAGGGTCAGCAAGCTGGTGACCATTGCAACGACCAGATCGGGGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGcc                                                                              <  2:1276855/150‑1 (MQ=255)
                                                                             gACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGt                                                                          >  1:435080/1‑150 (MQ=255)
                                                                                cAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTcgc                                                                       >  2:1489255/1‑150 (MQ=255)
                                                                                 aaCGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCCACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCAGGTcgcg                                                                      >  1:643350/1‑150 (MQ=255)
                                                                                 aaCGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGCTCTTTCACCCAAAATACCGGGCCGTGATCCCGCGCCATGTcgcg                                                                      >  2:1863389/1‑150 (MQ=255)
                                                                                 aaCGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTcgcg                                                                      <  1:1438002/150‑1 (MQ=255)
                                                                                    ggCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAAt                                                                   <  2:724217/150‑1 (MQ=255)
                                                                                       gTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGt                                                                >  2:166545/1‑150 (MQ=255)
                                                                                        tCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTc                                                               >  2:1087286/1‑150 (MQ=255)
                                                                                         ccAGCACCTGCAACATCCATTTTGGGACGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGGGACCATTGCAACGACCAGATCGGTTGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTcc                                                              <  2:118754/150‑1 (MQ=255)
                                                                                             cACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGGCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCccgcc                                                          <  2:267397/150‑1 (MQ=255)
                                                                                               ccGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGCCCATTGCAACGCCCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCcgccgc                                                        >  1:1395864/1‑150 (MQ=255)
                                                                                               ccGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCcgccgc                                                        <  2:750892/150‑1 (MQ=255)
                                                                                                   cAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATCGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAAAACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTAc                                                    <  1:911283/150‑1 (MQ=255)
                                                                                                     ccATCCCTTTTGCTAAGCCGCTAGGCAGGCGCACCAGCAAGATGTTAATCATGCTGGTGAGCATTGCAAGGACCAGATCTGTTTGATCTTTCACCCAAAAAACCGGGCAGTGATGCAGCGCCATGTCCCGAATTGTCCCGCCGCCTACTg                                                  <  2:396894/149‑1 (MQ=255)
                                                                                                      cATCCATTTTTGTACGCGTCTTGGCTGGCGCAACAGCACGATGGTCAGCATGCTGGGGACCATTGCAACGCCCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATGCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGc                                                 <  1:1304221/150‑1 (MQ=255)
                                                                                                      cATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGc                                                 >  2:1680139/1‑150 (MQ=255)
                                                                                                      cATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGc                                                 >  2:553850/1‑150 (MQ=255)
                                                                                                      cATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGc                                                 >  1:1806555/1‑150 (MQ=255)
                                                                                                        tCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCgg                                               <  1:80450/150‑1 (MQ=255)
                                                                                                        tCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCgg                                               <  2:1939171/150‑1 (MQ=255)
                                                                                                         ccATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGt                                              <  1:661239/150‑1 (MQ=255)
                                                                                                          cATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGtt                                             >  1:1682290/1‑150 (MQ=255)
                                                                                                             tttGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTAcc                                          >  1:984674/1‑150 (MQ=255)
                                                                                                             tttGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTAcc                                          <  2:1007969/150‑1 (MQ=255)
                                                                                                             tttGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTAcc                                          <  1:746250/150‑1 (MQ=255)
                                                                                                                ggTAAGCGTCTTGGCTGGCGCACCAGCCCGATGGTCAGCATGCTGGTGCCCATTGCACCGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGAGACAGCGCCATGGCGCGAAGTCTCCCGCCGCCTACTGCGGGTagcccg                                       >  2:562691/1‑145 (MQ=255)
                                                                                                                 gTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGc                                      >  1:1702083/1‑150 (MQ=255)
                                                                                                                 gTAAGCGTCTTCGCTGGCGCTCCAGCACGATGGGCAGCATGCTGGGGACCATTGCAACGACCAGAGCGGTGGGAGCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTTCTGCGGTTACCACGc                                      <  1:268159/150‑1 (MQ=255)
                                                                                                                   aaGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCAAGccc                                    >  2:1861118/1‑150 (MQ=255)
                                                                                                                    aGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCa                                   <  1:1695440/150‑1 (MQ=255)
                                                                                                                     gCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAg                                  >  1:6954/1‑150 (MQ=255)
                                                                                                                      cGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGCCCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGt                                 >  2:1477615/1‑150 (MQ=255)
                                                                                                                       gTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTa                                <  2:1796991/150‑1 (MQ=255)
                                                                                                                       gTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCCCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCTGTa                                >  1:178664/1‑150 (MQ=255)
                                                                                                                           tggctggcGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAg                            >  2:21895/1‑150 (MQ=255)
                                                                                                                            ggctggcGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGa                           <  1:2091000/150‑1 (MQ=255)
                                                                                                                               gggcGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAaca                        <  1:1779845/149‑1 (MQ=255)
                                                                                                                                 ggtcACCAGCATGCTGACCATCGTGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACAcc                      >  2:1923024/4‑150 (MQ=255)
                                                                                                                                  cgcACCCGCACGATGGTCAGCATGCTGGTGACCATTGCAACGCCCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTCCTGCGGTTACCACGCCCAGTCCCAGAACACCa                     >  2:1998309/1‑150 (MQ=255)
                                                                                                                                  cgcACCAGCACGATGGTCAGCATTCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCa                     >  2:196700/1‑150 (MQ=255)
                                                                                                                                  cgcACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCa                     >  2:904343/1‑150 (MQ=255)
                                                                                                                                     aCCAGCACGATTGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGACGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCaaaa                  <  2:942428/150‑1 (MQ=255)
                                                                                                                                      ccAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGCTCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACCCCAAAAg                 >  2:358341/1‑150 (MQ=255)
                                                                                                                                         ttgcaaTGGTCACCAGCATGCTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACGAAAAGGGt              >  1:1291908/6‑150 (MQ=255)
                                                                                                                                          cACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGTGCCGTGATCCAGCGCCAGGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGACCACCAAAAGGGTc             <  2:280888/150‑1 (MQ=255)
                                                                                                                                          cACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACAACGCCCAGTACCAGAACACCAAAAGGGTc             >  2:1714287/1‑150 (MQ=255)
                                                                                                                                           aCGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTcc            >  1:1195632/1‑150 (MQ=255)
                                                                                                                                            cGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTCCCAGAACACCAAAAGGGTCCa           >  1:233911/1‑150 (MQ=255)
                                                                                                                                            cGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCa           <  2:1973382/150‑1 (MQ=255)
                                                                                                                                              aTGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCATa         <  1:10399/150‑1 (MQ=255)
                                                                                                                                                 gTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCATACGc      >  2:1572174/1‑150 (MQ=255)
                                                                                                                                                     gCATGCTGGAGACCATTGCAACGACAAGAACGGTGGGATCTTACACCAAAAATACCCGGCCGTGATCCAGCGCCATGTCGCGAAATGTCCAGCCGCCTACTGCGGTAACCACGCCCAGTAACAGAACACCAAAAGGGTCGATAAGCAAtt  >  2:1745401/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCCATACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCATACGCAATT  >  NC_000913/177256‑177554

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