Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,002,463 A→G F113L (TTC→CTC)  yigF ← DUF2628 domain‑containing protein YigF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,002,4630AG98.9% 321.2 / ‑5.3 91F113L (TTC→CTC) yigFDUF2628 domain‑containing protein YigF
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (47/43);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCATAACCCCAGGAC  >  NC_000913/4002326‑4002598
                                                                                                                                         |                                                                                                                                       
cTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGcc                                                                                                                                       >  2:36719/1‑140 (MQ=255)
cTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGcc                                                                                                                                       >  1:1376598/1‑140 (MQ=255)
 tCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGATAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCACGTATTTAACCCTCTTACATCTTTCTCATAGAGCCa                                                                                                                                      >  2:696805/1‑140 (MQ=255)
 tCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCAAGGTATTAAACACTGTTACCTATTTCTCATAGAGCCa                                                                                                                                      >  1:646961/1‑140 (MQ=255)
   gcgGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCAtt                                                                                                                                    <  1:956255/140‑1 (MQ=255)
               ccGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTaa                                                                                                                        <  2:1376598/140‑1 (MQ=255)
                 gTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACa                                                                                                                      <  1:36719/140‑1 (MQ=255)
                 gTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACa                                                                                                                      >  2:828021/1‑140 (MQ=255)
                  tATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAg                                                                                                                     >  2:849933/1‑140 (MQ=255)
                     cTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCaa                                                                                                                  >  2:768147/1‑140 (MQ=255)
                       gTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACa                                                                                                                <  1:334312/140‑1 (MQ=255)
                            ttAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGgc                                                                                                           >  1:167854/1‑140 (MQ=255)
                             tAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGgct                                                                                                          >  1:1160016/1‑140 (MQ=255)
                             tAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGgct                                                                                                          >  1:914787/1‑140 (MQ=255)
                             tAAGCGGAAGAACTAGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGgct                                                                                                          >  2:967109/1‑140 (MQ=255)
                               aGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAAc                                                                                                                       <  2:1185993/125‑1 (MQ=255)
                               aGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAAc                                                                                                                       >  1:1185993/1‑125 (MQ=255)
                               aGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGgctgc                                                                                                        >  2:879281/1‑140 (MQ=255)
                               aGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCCGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGgctgc                                                                                                        >  2:823143/1‑140 (MQ=255)
                                  ggAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTAc                                                                                                     >  2:175651/1‑140 (MQ=255)
                                    aagaaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGGCATTTGTTCGTATTAACAGCAACATAGGCTGCTACgg                                                                                                   >  1:495994/1‑140 (MQ=255)
                                     agaaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGc                                                                                                  <  1:879281/140‑1 (MQ=255)
                                           ggTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTcc                                                                                            <  1:988035/140‑1 (MQ=255)
                                             ttGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCa                                                                                          <  1:828021/140‑1 (MQ=255)
                                              tGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAg                                                                                         <  1:801008/140‑1 (MQ=255)
                                               gTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGa                                                                                        <  1:1040227/140‑1 (MQ=255)
                                                 aaaaaaCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAAt                                                                                      >  1:268092/1‑140 (MQ=255)
                                                           agCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACGCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCa                                                                            <  2:766231/140‑1 (MQ=255)
                                                             cGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAAt                                                                          >  2:1349094/1‑140 (MQ=255)
                                                             cGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAAt                                                                          <  1:1169758/140‑1 (MQ=255)
                                                                    ttAGTCTGGCTAAGCCACCTGTCATTCAAAGTCTCCAGGCATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCa                                                                   <  2:495994/140‑1 (MQ=255)
                                                                       gTCTGGCTAAGCCACCTGGTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACca                                                                >  2:1302646/1‑140 (MQ=255)
                                                                        tCTGGCTAAGCCACCTGTTATTTAAAGGCTCTAGGAATTCAACCCTTTCACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaa                                                               <  2:405105/140‑1 (MQ=255)
                                                                        tCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCTCAGTATATTAACGACAATTAGTGCACcaa                                                               >  1:405105/1‑140 (MQ=255)
                                                                          tGGCTAAGCCACCTGTTATTCAAAGGCTCCATGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaaca                                                             <  2:646961/140‑1 (MQ=255)
                                                                          tGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTACTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaaca                                                             <  1:768147/140‑1 (MQ=255)
                                                                            gCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACCTAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACaaa                                                           >  2:212810/1‑140 (MQ=255)
                                                                            gCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCATTTCCCAGAATGTTAACGCCAATTAGTGCACCAACaaa                                                           >  1:296393/1‑140 (MQ=255)
                                                                            gCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACaaa                                                           >  1:822480/1‑140 (MQ=255)
                                                                            gCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACaaa                                                           >  2:1135593/1‑140 (MQ=255)
                                                                             cTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAAt                                                          <  2:167854/140‑1 (MQ=255)
                                                                               aaGCCACCTGTTATTGAAAGGCACCTGGTATTTAACCTTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATag                                                        <  1:604715/140‑1 (MQ=255)
                                                                                 gCCTCCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaa                                                               <  2:99081/131‑1 (MQ=255)
                                                                                 gCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaa                                                               >  1:99081/1‑131 (MQ=255)
                                                                                 gCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACaaa                                                           >  1:669887/1‑135 (MQ=255)
                                                                                 gCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACaaa                                                           <  2:669887/135‑1 (MQ=255)
                                                                                  ccACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCATTTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAgaga                                                     <  2:296393/140‑1 (MQ=255)
                                                                                      cTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaa                                                               >  2:303926/1‑126 (MQ=255)
                                                                                      cTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaa                                                               <  1:303926/126‑1 (MQ=255)
                                                                                      cTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGa                                                 >  1:1195956/1‑140 (MQ=255)
                                                                                         ttATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATa                                              <  1:1014401/140‑1 (MQ=255)
                                                                                         ttATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATa                                              <  2:1450798/140‑1 (MQ=255)
                                                                                         ttATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATa                                              <  2:1195956/140‑1 (MQ=255)
                                                                                         ttATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATa                                              <  2:268092/140‑1 (MQ=255)
                                                                                          tATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATAc                                             <  2:822480/140‑1 (MQ=255)
                                                                                           aTTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATAcc                                            <  2:1240214/140‑1 (MQ=255)
                                                                                            ttCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCa                                           <  1:1349094/140‑1 (MQ=255)
                                                                                                 aGGCTCCAGGTATTTATCCCTTGTACCTCTTCCGCCTAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTaa                                      <  2:696989/140‑1 (MQ=255)
                                                                                                 aGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAGTAGAGATAGAATACCAAGTaa                                      <  1:849933/140‑1 (MQ=255)
                                                                                                 aGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATa                                              >  1:1126497/1‑132 (MQ=255)
                                                                                                 aGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATa                                              <  2:1126497/132‑1 (MQ=255)
                                                                                                 aGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTaa                                      <  1:3305/140‑1 (MQ=255)
                                                                                                 aGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTaa                                      >  1:834943/1‑140 (MQ=255)
                                                                                                 aGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTaa                                      <  1:1309442/140‑1 (MQ=255)
                                                                                                  ggCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAAt                                     <  2:834943/140‑1 (MQ=255)
                                                                                                  ggCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAAt                                     >  2:1341090/1‑140 (MQ=255)
                                                                                                  ggCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAAt                                     <  1:823143/140‑1 (MQ=255)
                                                                                                   gCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATAcc                                            >  2:1278363/1‑132 (MQ=255)
                                                                                                   gCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATAcc                                            <  1:1278363/132‑1 (MQ=255)
                                                                                                   gCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATa                                    >  2:307107/1‑140 (MQ=255)
                                                                                                     tCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGa                                  <  1:1135593/140‑1 (MQ=255)
                                                                                                      ccAGGTATTTAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGAGCACCAAAAAACAGGGATAGAATACCAAGTAATAGaa                                 >  1:696989/1‑140 (MQ=255)
                                                                                                             tttAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCa                                           <  2:122634/123‑1 (MQ=255)
                                                                                                             tttAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCa                                           >  1:122634/1‑123 (MQ=255)
                                                                                                              ttAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAAtt                         >  2:1184671/1‑140 (MQ=255)
                                                                                                               tAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCACCAAGTAGTGATAGAATACCAATAAATAGAACTATAATTg                        >  2:724717/1‑140 (MQ=255)
                                                                                                               tAACCCTTTTACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTg                        >  2:122487/1‑140 (MQ=255)
                                                                                                                      tttACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCtttttt                 >  2:219994/1‑140 (MQ=255)
                                                                                                                      tttACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCtttttt                 >  2:1185646/1‑140 (MQ=255)
                                                                                                                      tttACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAAATAGAATGCCATGTGATAGAAATATAATTGCtttttt                 >  1:633880/1‑140 (MQ=255)
                                                                                                                       ttACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGGTAGAATATCAATTATTAGAACTATAATTGCTTTTTTc                >  1:461059/1‑140 (MQ=255)
                                                                                                                       ttACCTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTc                >  2:42023/1‑140 (MQ=255)
                                                                                                                          ccTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCAt             <  1:175651/140‑1 (MQ=255)
                                                                                                                          ccTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCAt             <  1:307107/140‑1 (MQ=255)
                                                                                                                          ccTCTTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCAt             <  1:42023/140‑1 (MQ=255)
                                                                                                                          ccTCTTTCTCATAGAGCCATTTATTAGTGTTAACAGCCACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCAt             <  1:724717/140‑1 (MQ=255)
                                                                                                                             ctTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCAt             >  1:1482367/1‑137 (MQ=255)
                                                                                                                             ctTTCTCATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCAt             <  2:1482367/137‑1 (MQ=255)
                                                                                                                                tctcATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCATAAcccc       >  2:823703/1‑140 (MQ=255)
                                                                                                                                tctcATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCATAAcccc       >  1:1404038/1‑140 (MQ=255)
                                                                                                                                 ctcATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCtt                     <  1:1548127/125‑1 (MQ=255)
                                                                                                                                 ctcATAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCtt                     >  2:1548127/1‑125 (MQ=255)
                                                                                                                                     tAGAGCCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCATAACCCCAGGAc  >  1:954022/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                       
CTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCATAACCCCAGGAC  >  NC_000913/4002326‑4002598

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