Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,547,319	C→A	100%	R742S (CGT→AGT)	fimD →	fimbrial usher outer membrane porin protein; FimCD chaperone‑usher

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,547,319	0	C	A	100.0%	321.6 / ‑6.7	98	R742S (CGT→AGT)	fimD	fimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base A (41/55);  minor base G (0/1);  total (41/57)

GCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGAT  >  NC_000913/4547222‑4547411                                                                                                  |                                                                                             gCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGt                                                                                            <  2:1607467/100‑1 (MQ=255)   tGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGTCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTgg                                                                                          <  2:293310/100‑1 (MQ=255)   tGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTgg                                                                                          <  2:776054/100‑1 (MQ=255)   tGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTgg                                                                                          <  1:1282941/100‑1 (MQ=255)     gggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGTCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGtt                                                                                        <  2:1638478/100‑1 (MQ=255)         aGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATg                                                                                     <  2:1475483/99‑1 (MQ=255)         aGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGc                                                                                    >  1:197551/1‑100 (MQ=255)         aGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGc                                                                                    <  2:1332128/100‑1 (MQ=255)         aGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGc                                                                                    <  2:1612733/100‑1 (MQ=255)         aGCCGGTAAACGATACGGTCGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAGGTCGAAAACCAGACCGGGCTGCGTACCGACTGGAGTg                                                                                           <  2:1905410/93‑1 (MQ=255)          gCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATg                                                                                     >  2:1870289/1‑98 (MQ=255)          gCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGcc                                                                                   >  1:1447758/1‑100 (MQ=255)           ccGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCg                                                                                  <  1:2307961/100‑1 (MQ=255)           ccGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCg                                                                                  <  2:2103293/100‑1 (MQ=255)            cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGt                                                                                 >  2:873424/1‑100 (MQ=255)            cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGCGTGCGTACCGACTGGAGTGGTTATGCCGt                                                                                 >  2:873305/1‑100 (MQ=255)             gTTAAACGATACGGTGGTGCTTGTTAGAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtg                                                                                >  2:2311805/1‑100 (MQ=255)             gTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtg                                                                                >  2:2433383/1‑100 (MQ=255)               tAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgct                                                                              >  1:1223158/1‑100 (MQ=255)                aaaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgctg                                                                             <  2:505632/100‑1 (MQ=255)                 aaCTATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGAAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtg                                                                                <  2:959575/96‑1 (MQ=255)                 aaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgctgc                                                                            >  1:555811/1‑100 (MQ=255)                  aCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGcc                                                                           <  2:2369831/100‑1 (MQ=255)                  aCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGcc                                                                           <  1:302082/100‑1 (MQ=255)                  aCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGcc                                                                           <  2:2452222/100‑1 (MQ=255)                  aCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGcc                                                                           <  1:2224611/100‑1 (MQ=255)                    gATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCtt                                                                         >  1:1457453/1‑100 (MQ=255)                    gATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCtt                                                                         >  1:2279549/1‑100 (MQ=255)                     aTACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTa                                                                        >  1:1501995/1‑100 (MQ=255)                       aCGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATg                                                                      <  1:2311819/100‑1 (MQ=255)                       aCGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATg                                                                      <  2:611297/100‑1 (MQ=255)                        cGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGc                                                                     >  1:1390064/1‑100 (MQ=255)                         ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGcc                                                                    <  1:2283505/100‑1 (MQ=255)                           tggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCAc                                                                  <  1:1882781/100‑1 (MQ=255)                            ggtgCTAGTTAAAGCGCCTGGCGCACAAGATGCAAAAGCCGAAAACCAGCCGGGGGTGCTTACCGACTGCACTGGATACACCGTGCTGCCTTATGCCAc                                                                  >  2:2248097/1‑99 (MQ=255)                             gtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTg                                                                <  2:850289/100‑1 (MQ=255)                              tgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGa                                                               >  2:2237270/1‑100 (MQ=255)                               gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                              >  2:832878/1‑100 (MQ=255)                               gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                              >  1:1771966/1‑100 (MQ=255)                                cTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAat                                                             <  1:141940/100‑1 (MQ=255)                                cTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAat                                                             <  1:740705/100‑1 (MQ=255)                                  tGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAatat                                                           >  2:983641/1‑100 (MQ=255)                                  tGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAatat                                                           <  1:357188/100‑1 (MQ=255)                                   gTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATc                                                          >  2:1414913/1‑100 (MQ=255)                                    ttAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGACGTGCTGCCTTATGGCACTGAATATCg                                                         >  2:1269778/1‑100 (MQ=255)                                     tAAAGCGGCTGGCGCAAAAGATGCAAAAGTGGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                        <  2:2248200/100‑1 (MQ=255)                                     tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                        <  1:2237284/100‑1 (MQ=255)                                     tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                        >  2:921512/1‑100 (MQ=255)                                     tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                        >  1:1283298/1‑100 (MQ=255)                                      aaaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGTTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                        <  1:1472424/99‑1 (MQ=255)                                      aaaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCggg                                                       >  2:745820/1‑100 (MQ=255)                                        aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaa                                                     <  1:2119380/100‑1 (MQ=255)                                        aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaa                                                     >  1:700617/1‑100 (MQ=255)                                         gcgcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaaa                                                    >  2:819728/1‑100 (MQ=255)                                         gcgcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGAGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaaa                                                    >  1:1486093/1‑100 (MQ=255)                                           gcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAAt                                                  >  2:1255215/1‑100 (MQ=255)                                            ccTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATa                                                 >  1:2431990/1‑100 (MQ=255)                                               ggggcAAAAGAGGCAAAAGGCGAAAGCCAGACGGGGGTGCGGACCGACTGGGGTGGGTATGCCGTGGTGCCTTATGCCACTGAATATCGGGAAAAt                                                  <  2:1740808/93‑1 (MQ=255)                                               ggCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATagag                                              <  1:1418555/100‑1 (MQ=255)                                                gcgcAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGt                                             <  1:393104/100‑1 (MQ=255)                                                  gcAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTgg                                           <  2:1758365/100‑1 (MQ=255)                                                    aaaaGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgcg                                         <  1:599012/100‑1 (MQ=255)                                                         aTGCAAAAGTCGAAACCCAGACGGGGGTGAGTACCGACTGGAGTGGTTATGCCGTGTTGCCGTATGCCACTGAATATCGGGTAAATGGAGCGGCGCTgg                                     >  2:1489703/1‑99 (MQ=255)                                                          tGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGAt                                   >  1:2224442/1‑100 (MQ=255)                                                           gCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATa                                  <  2:135198/100‑1 (MQ=255)                                                            cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGGCTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                 >  2:1693337/1‑100 (MQ=255)                                                            cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                 >  1:1083568/1‑100 (MQ=255)                                                            cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                 >  2:2009050/1‑100 (MQ=255)                                                             aaaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAcc                                <  1:1309126/100‑1 (MQ=255)                                                             aaaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAcc                                <  2:2391882/100‑1 (MQ=255)                                                              aaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCa                               <  2:1559047/100‑1 (MQ=255)                                                                 gTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATa                            >  1:58010/1‑100 (MQ=255)                                                                   cGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAcc                          <  1:545801/100‑1 (MQ=255)                                                                   cGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAcc                          >  1:1247077/1‑100 (MQ=255)                                                                   cGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAcc                          >  1:1745867/1‑100 (MQ=255)                                                                   cGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAc                           <  2:1277394/99‑1 (MQ=255)                                                                        aCCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGc                     >  2:62601/1‑100 (MQ=255)                                                                         ccAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCt                    <  2:1090823/100‑1 (MQ=255)                                                                         ccAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCt                    <  2:327901/100‑1 (MQ=255)                                                                            gACGGGGGGGCGTACCGACTGGCGTGGTTATGGCGTGCTGCCGTATGCCACTGAATAGCGGGAAAAGAGAGTGGCGCTGGATACCAATACCCTg                       <  2:1524294/94‑1 (MQ=255)                                                                              cGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATaa               >  1:1237687/1‑100 (MQ=255)                                                                               gggggTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAAc              <  1:1145800/100‑1 (MQ=255)                                                                                ggggTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACg             <  2:1615445/100‑1 (MQ=255)                                                                                  ggTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGt            <  2:2271762/99‑1 (MQ=255)                                                                                  ggTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTc           <  1:283616/100‑1 (MQ=255)                                                                                   gTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCg          >  1:2183094/1‑100 (MQ=255)                                                                                   gTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACAAATACCCTGGCTGATAACGTCg          <  1:1624244/100‑1 (MQ=255)                                                                                     gCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAt        <  1:2101340/100‑1 (MQ=255)                                                                                      cGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAtt       >  2:77034/1‑100 (MQ=255)                                                                                       gTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt      <  1:1528364/100‑1 (MQ=255)                                                                                        tACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTa     <  1:808889/100‑1 (MQ=255)                                                                                        tACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTa     <  1:1004087/100‑1 (MQ=255)                                                                                          ccGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAGTAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGa   <  1:2398200/100‑1 (MQ=255)                                                                                          ccGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACGCTGGCTGATAACGTCGATTTAg    >  2:924433/1‑99 (MQ=255)                                                                                          ccGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGa   <  1:745827/100‑1 (MQ=255)                                                                                           cGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGAt  <  1:1543721/100‑1 (MQ=255)                                                                                           cGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGAt  <  2:263746/100‑1 (MQ=255)                                                                                           cGACTGGAGTGGTTATGCCGTGCTGCCATATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGAt  <  1:2137306/100‑1 (MQ=255)                                                                                                  |                                                                                             GCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGAT  >  NC_000913/4547222‑4547411

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