Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,542,645 G→T intergenic (+12/‑470) fimE → / → fimA regulator for fimA/type 1 fimbriae major subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,542,6450GT98.8% 301.5 / ‑3.1 85intergenic (+12/‑470)fimE/fimAregulator for fimA/type 1 fimbriae major subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (41/43);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTA  >  NC_000913/4542513‑4542780
                                                                                                                                    |                                                                                                                                       
tctcGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTAtaataa                                                                                                                                  <  1:2184647/140‑1 (MQ=255)
tctcGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTAtaataa                                                                                                                                  >  2:827074/1‑140 (MQ=255)
tctcGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTAtaataa                                                                                                                                  >  2:529042/1‑140 (MQ=255)
tctcGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGCTTGATTTAACTTATTTAtaataa                                                                                                                                  >  1:1602099/1‑140 (MQ=255)
     ggCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTa                                                                                                                             >  1:2038258/1‑140 (MQ=255)
     ggCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTa                                                                                                                             >  1:1375520/1‑140 (MQ=255)
     ggCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTa                                                                                                                             >  1:1559549/1‑140 (MQ=255)
         tCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTaaaaa                                                                                                                         <  1:726547/140‑1 (MQ=255)
          cGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTaaaaaa                                                                                                                        >  2:930425/1‑140 (MQ=255)
          cGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTaaaaaa                                                                                                                        >  2:1775610/1‑140 (MQ=255)
           gAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTaaaaaaa                                                                                                                       <  1:1055687/140‑1 (MQ=255)
           gAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTaaaaaaa                                                                                                                       <  2:2116319/140‑1 (MQ=255)
           gAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTaaaaaaa                                                                                                                       <  2:67395/140‑1 (MQ=255)
                atTCGCCATATTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACaaat                                                                                                                  <  1:935824/140‑1 (MQ=255)
                atTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACaaat                                                                                                                  <  1:394618/140‑1 (MQ=255)
                atTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACaaat                                                                                                                  <  2:2004161/140‑1 (MQ=255)
                     ccATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAaataaata                                                                                                             <  1:529042/140‑1 (MQ=255)
                     ccATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAaataaata                                                                                                             <  1:930425/140‑1 (MQ=255)
                             tgCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACaaat                                                                                                                  >  2:1606502/1‑127 (MQ=255)
                             tgCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACaaat                                                                                                                  <  1:1606502/127‑1 (MQ=255)
                             tgCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACaa                                                                                                     >  2:2169473/1‑140 (MQ=255)
                             tgCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACaa                                                                                                     >  2:187520/1‑140 (MQ=255)
                             tgCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACaa                                                                                                     >  1:1339324/1‑140 (MQ=255)
                             tgCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACaa                                                                                                     >  2:1259683/1‑140 (MQ=255)
                             tgCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGATTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACaa                                                                                                     >  1:529068/1‑140 (MQ=255)
                                  tataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACaaataa                                                                                                                <  2:915603/124‑1 (MQ=255)
                                  tataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACaaataa                                                                                                                >  1:915603/1‑124 (MQ=255)
                                  tataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTggg                                                                                                >  2:1010493/1‑140 (MQ=255)
                                  tataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTggg                                                                                                >  1:2054700/1‑140 (MQ=255)
                                  tataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTggg                                                                                                >  1:275838/1‑140 (MQ=255)
                                  tataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTggg                                                                                                >  1:339561/1‑140 (MQ=255)
                                  tataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTggg                                                                                                >  2:1602411/1‑140 (MQ=255)
                                   ataCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTgggg                                                                                               <  2:2239781/140‑1 (MQ=255)
                                     cccGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGcc                                                                                             <  1:1529825/139‑1 (MQ=255)
                                       cGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                                                           <  2:1184468/140‑1 (MQ=255)
                                       cGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                                                           >  2:979969/1‑140 (MQ=255)
                                       cGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                                                           <  1:2169473/140‑1 (MQ=255)
                                        gCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaaa                                                                                          <  1:1090051/140‑1 (MQ=255)
                                        gCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaaa                                                                                          <  1:1667246/140‑1 (MQ=255)
                                        gCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaaa                                                                                          <  2:2038258/140‑1 (MQ=255)
                                        gCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaaa                                                                                          <  2:2054700/140‑1 (MQ=255)
                                        gCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaaa                                                                                          <  2:1339324/140‑1 (MQ=255)
                                        gCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaaa                                                                                          <  1:1602411/140‑1 (MQ=255)
                                        gCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaaa                                                                                          <  2:339561/140‑1 (MQ=255)
                                                 gctgctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCata                                                                                 <  2:2072397/140‑1 (MQ=255)
                                                  ctgctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCatat                                                                                <  1:2227023/140‑1 (MQ=255)
                                                    gctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCa                                                                              <  2:105354/140‑1 (MQ=255)
                                                    gctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCa                                                                              >  1:1914859/1‑140 (MQ=255)
                                                    gctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCa                                                                              >  2:529336/1‑140 (MQ=255)
                                                     ctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCat                                                                             >  1:2241089/1‑140 (MQ=255)
                                                       cGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataa                                                                           >  2:1896264/1‑140 (MQ=255)
                                                        gTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaa                                                                          <  2:2101234/140‑1 (MQ=255)
                                                              ccGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGtt                                                                    <  2:1559549/140‑1 (MQ=255)
                                                               cGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTa                                                                   >  1:1585791/1‑140 (MQ=255)
                                                                  aTTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGt                                                                >  2:1850138/1‑140 (MQ=255)
                                                                       gggAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAtttt                                                           >  1:544026/1‑140 (MQ=255)
                                                                           aagaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtct                                                       <  1:924074/140‑1 (MQ=255)
                                                                           aagaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtct                                                       >  1:1747561/1‑140 (MQ=255)
                                                                            agaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                      <  1:529336/140‑1 (MQ=255)
                                                                            agaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                      >  1:1893615/1‑140 (MQ=255)
                                                                            agaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                      >  1:179907/1‑140 (MQ=255)
                                                                                ataatCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGc                                                  >  2:834658/1‑140 (MQ=255)
                                                                                ataatCTCATAAACGAAAAATTAAAAAGAGAAGAGGGTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGc                                                  >  1:498768/1‑140 (MQ=255)
                                                                                    tctcATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAACAAATAGATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa                                             <  1:1914599/140‑1 (MQ=255)
                                                                                        aTAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAat                                          >  1:6746/1‑140 (MQ=255)
                                                                                        aTAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAat                                          >  2:2328710/1‑140 (MQ=255)
                                                                                         tAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACCTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAata                                         >  2:1932860/1‑140 (MQ=255)
                                                                                                aaaaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaaa                                  <  1:1956119/140‑1 (MQ=255)
                                                                                                  aaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAc                                >  2:53362/1‑140 (MQ=255)
                                                                                                        aaaGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaa                          <  1:1850138/140‑1 (MQ=255)
                                                                                                          agagaagagGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaaa                        <  1:1259683/140‑1 (MQ=255)
                                                                                                             gaagagGTTTGATTTAACCTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCa                     <  1:1932860/140‑1 (MQ=255)
                                                                                                              aagagGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCAt                    >  2:15612/1‑140 (MQ=255)
                                                                                                              aagagGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCAt                    >  1:1547212/1‑140 (MQ=255)
                                                                                                                gagGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCt                  <  2:2241089/140‑1 (MQ=255)
                                                                                                                 agGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa                 <  2:544026/140‑1 (MQ=255)
                                                                                                                 agGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa                 <  1:1524385/140‑1 (MQ=255)
                                                                                                                 agGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa                 <  1:1010493/140‑1 (MQ=255)
                                                                                                                  ggTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTaa                <  1:15612/140‑1 (MQ=255)
                                                                                                                  ggTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTaa                <  2:1585791/140‑1 (MQ=255)
                                                                                                                   gcttGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAAc               <  2:1602099/138‑1 (MQ=255)
                                                                                                                         tttAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTg         >  2:1517896/1‑140 (MQ=255)
                                                                                                                           tAACTTATTTATAATAAAGTTAAAAAACCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGat       >  2:1115214/1‑140 (MQ=255)
                                                                                                                             aCTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGatat     >  1:1348898/1‑140 (MQ=255)
                                                                                                                                tATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTa  <  2:1914859/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
TCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTA  >  NC_000913/4542513‑4542780

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