Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,544,588	G→A	7.2%	G95E (GGA→GAA)	fimC →	periplasmic chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,588	0	G	A	7.2%	278.1 / 2.1	97	G95E (GGA→GAA)	fimC	periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base G (46/43);  new base A (2/6);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATC  >  NC_000913/4544441‑4544736                                                                                                                                                     |                                                                                                                                                     tttATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGa                                                                                                                                                     >  2:1323110/1‑149 (MQ=255)  ttATCCGGCAGGGCAAAAACAAGAGCAAGTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAGTTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGaa                                                                                                                                                    >  1:878998/1‑149 (MQ=255)   tATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAGG‑GGaaa                                                                                                                                                   >  2:942009/1‑149 (MQ=255)      ccGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGCGTAAGGGATGGTCGTGGTATCGTGACGCCGCCTCTGTTTGCGATGAAG‑GGaaaaaa                                                                                                                                                >  1:1425250/1‑149 (MQ=255)         ggaGGGCAAAAACAAGACAAACTTGCATGGACAAATAATGATGAAAATAGTAACTAGTGAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGGGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAaga                                                                                                                                             <  2:390715/147‑1 (MQ=255)         gCAGGGCAAAAACAAGAGCCACTTGCCGCGACAAATCATTATGAAAAGAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAaga                                                                                                                                             <  1:1168458/149‑1 (MQ=255)         gCAAGGAAAAAACGAGAGCAACGTGCCGGGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAaga                                                                                                                                             <  2:1369488/149‑1 (MQ=255)             ggCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTGTCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAACACAGAGAAt                                                                                                                                         >  2:1194433/1‑149 (MQ=255)             ggCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGGCGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAAAGAAt                                                                                                                                         >  2:33244/1‑149 (MQ=255)             ggCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGCTGAAG‑GGAAAAAAAGAGAAt                                                                                                                                         >  1:1022092/1‑149 (MQ=255)             ggCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGCAAAAAAGAGAAt                                                                                                                                         >  2:1117734/1‑149 (MQ=255)              gCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATa                                                                                                                                        >  2:1307572/1‑149 (MQ=255)                aaaaaCAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATAcc                                                                                                                                      >  2:549460/1‑149 (MQ=255)                 aaaaCAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCt                                                                                                                                     <  1:934250/149‑1 (MQ=255)                 aaaaCAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCt                                                                                                                                     >  1:447381/1‑149 (MQ=255)                   aaCAAGAGAAACTTGCCGTGACAAAGAAAGAAAAAAATAGAACCGAATAAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGAAAAAAAA‑AAAAAAAAAGAGAATACCTTa                                                                                                                                   <  2:375646/149‑1 (MQ=255)                    aCAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCGTAc                                                                                                                                  >  1:83187/1‑149 (MQ=255)                                tGCCGTGACAAATAATGATGAAAATAGTGCCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATg                                                                                                                      <  1:917033/149‑1 (MQ=255)                                     tGACAAATAATGATGAAAATAGTACCTATTTAATTCAAACATGGGTGGAAAATGCCGATGGCGTCAAGGGTGGTCGTTATATCGTGACGCCGCCTCTGTTTGCGATGAAG‑GGAAACAAAGGGGATACCTTACGTATTCTTCATGcaaca                                                                                                                 >  2:320791/1‑149 (MQ=255)                                       aCAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGATCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACaaa                                                                                                               >  2:342391/1‑149 (MQ=255)                                         aaaTAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATa                                                                                                             >  1:594258/1‑149 (MQ=255)                                           ataatGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTCACGATGAAG‑GACAACAACGCGAATGCCATACGTATGCCGGAGGCAACAAATAAc                                                                                                           >  2:1370974/1‑149 (MQ=255)                                                gatgaAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAAtt                                                                                                      <  1:17717/149‑1 (MQ=255)                                                gatgaAAATAGTACCTATTTAATTCAATCAGGGGGGGAAAAGGCCGATGGTGTAAAGAATGGACGTTTTATCGTGACGCCTCTTCTGTTTGAAAAAAAA‑AAAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAAtt                                                                                                      <  2:518009/149‑1 (MQ=255)                                                  tgaAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGc                                                                                                    <  1:1153723/149‑1 (MQ=255)                                                  agaAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGAAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGc                                                                                                    <  1:406585/148‑1 (MQ=255)                                                      aaTAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCcaca                                                                                                <  1:1323110/149‑1 (MQ=255)                                                            aCCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCg                                                                                          <  1:522158/149‑1 (MQ=255)                                                                aTTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGCTGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGac                                                                                      >  1:260620/1‑148 (MQ=255)                                                                aTTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGGCCGGGaa                                                                                      >  1:555126/1‑149 (MQ=255)                                                                aTTTAATTAAATCAGGGGTGGAAAATGCCGAAGGCGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCAAAAAAA‑AAAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGaa                                                                                      <  2:776432/149‑1 (MQ=255)                                                                 tttAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGaaa                                                                                     <  1:475766/149‑1 (MQ=255)                                                                 tttAATTCAATCATGGGGGGAAAATGCCTGTGGTGTAAAGGATGGTCGTTTTAGCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAAGTGCCACAGGACCGGGaaa                                                                                     <  2:194302/149‑1 (MQ=255)                                                                  gcaatagaaTCATGAGTAGAAAAAGCCGAGGCTCAAAACAATGGTCGTTTTACCGTGACGCCCCCTCTGCTTGAAAAAAAGAAAAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGaaa                                                                                     <  2:73137/142‑1 (MQ=255)                                                                   tAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCCCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGCATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGt                                                                                   <  1:784777/149‑1 (MQ=255)                                                                    aaTTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGtt                                                                                  >  1:1381507/1‑149 (MQ=255)                                                                    aaTTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGtt                                                                                  <  1:1285420/149‑1 (MQ=255)                                                                      ttCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTa                                                                                <  1:867383/149‑1 (MQ=255)                                                                       tCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAt                                                                               <  1:843246/149‑1 (MQ=255)                                                                        cAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTGACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTAtt                                                                              <  1:1369688/149‑1 (MQ=255)                                                                         aaTCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATCACCAATTGCCACAGGACCGGGAAAGTTTACTc                                                                             >  2:337000/1‑149 (MQ=255)                                                                         aaTCATGGGGGGAAAATGCCGATGGTGTAAAGGATGGTCGGTTTGTCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATATCCTTTTGCCACAGGACCGGGAAAGTTTATTc                                                                             <  2:1404456/149‑1 (MQ=255)                                                                           tCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTg                                                                           >  1:1163953/1‑149 (MQ=255)                                                                                 gTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGCAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGaa                                                                     >  1:921820/1‑149 (MQ=255)                                                                                 gTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGaa                                                                     <  2:634385/149‑1 (MQ=255)                                                                                  tGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGCATTCTTGATGCAACTGATAGCCAATTTCCACGGGACCGGGAAAGTTTATTCTGGATGAAc                                                                    <  1:938255/149‑1 (MQ=255)                                                                                  tGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCGGTTTGCGAGGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAAc                                                                    >  1:1002465/1‑149 (MQ=255)                                                                                   ggAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGGTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACg                                                                   <  1:678718/149‑1 (MQ=255)                                                                                    gAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAGTAACCAATGGCCACAGGACCAGGAAAGTTTATTTTGGATGAACGt                                                                  >  2:108615/1‑149 (MQ=255)                                                                                    gAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCCCAGGACCGGGAAAGTTTATTCTGGATGAACGt                                                                  >  1:1276172/1‑149 (MQ=255)                                                                                    gAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGt                                                                  <  1:1310728/149‑1 (MQ=255)                                                                                       aaTGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAACCTACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTaa                                                               >  2:522428/1‑149 (MQ=255)                                                                                        aTGCCGATGGTGTAAAGGTTGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTaaa                                                              <  1:549460/149‑1 (MQ=255)                                                                                            cGATGGTGTAAAGGAGGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACATATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGa                                                          <  2:700361/149‑1 (MQ=255)                                                                                             gATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGAt                                                         <  1:1320860/149‑1 (MQ=255)                                                                                               tGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTc                                                       <  1:1325802/149‑1 (MQ=255)                                                                                                    tAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCa                                                  >  2:1415869/1‑149 (MQ=255)                                                                                                      aaGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAAt                                                >  2:603753/1‑149 (MQ=255)                                                                                                      aaGGATGGTCGTTCTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAAt                                                >  2:1012322/1‑149 (MQ=255)                                                                                                           tGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATAGCTTACGTATTCTTGATGCAACAAATAACCAATTGCCGCAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTGAATGGATa                                           >  1:185077/1‑149 (MQ=255)                                                                                                              tCGTTTTATTGTGACGTCTCTTCTGTGTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAAt                                        <  1:58383/149‑1 (MQ=255)                                                                                                              tCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGCAAAAAAGAGAATACCTTACGTATTCTTGATGCAACGAATAACCAATTGCCACAGGCCCGGGACAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATaact                                        >  2:546/1‑147 (MQ=255)                                                                                                              tCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAAt                                        <  1:747024/149‑1 (MQ=255)                                                                                                                 ttttATCGTGACGCCTCCTCTGTTGGCGATGAAG‑TGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCaa                                     >  1:738686/1‑149 (MQ=255)                                                                                                                 tattAACGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATACTGGAAGAACGTTAAAGCGATTCCGTCAATGGAAAAATCaa                                     >  2:97675/3‑149 (MQ=255)                                                                                                                   ttATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCaaaa                                   <  1:1447484/149‑1 (MQ=255)                                                                                                                   ttATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GAAAAGAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATGCTGGATGAACGTTAAAGCGATTCCGTCAATGGAAAAATCaaaa                                   >  2:811312/1‑149 (MQ=255)                                                                                                                     aTCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATGCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAAtt                                 >  1:127268/1‑149 (MQ=255)                                                                                                                      tCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGGCCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTg                                >  2:1157446/1‑149 (MQ=255)                                                                                                                         tGACGCCTCCTCTGTTTGCGATGCAT‑TGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTTTGGATGAACGTTAAACCGATTCCGTCACTGGATAAGTCCAAATTGaat                             >  1:292007/1‑147 (MQ=255)                                                                                                                         tGACGCCTCCTCTGGTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCGGTCAATGGATAAATCAAAATTGACt                             <  1:868156/149‑1 (MQ=255)                                                                                                                           aCGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTga                           >  1:1232781/1‑149 (MQ=255)                                                                                                                           aCGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTga                           <  2:1163953/149‑1 (MQ=255)                                                                                                                           aCGCCTCCTCTGTTTGAGATGAAG‑GGAAAAAAAGAGAATACCTTACGGATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGCCAATGGATAAATCAAAATTGACTga                           <  1:320791/149‑1 (MQ=255)                                                                                                                           aCGCCCCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTga                           <  1:503328/149‑1 (MQ=255)                                                                                                                             gCCTCCTCTGTTTGCAATAAAA‑GAAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTgaga                         <  2:555126/149‑1 (MQ=255)                                                                                                                               ctcctcTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAAt                       >  1:1095128/1‑149 (MQ=255)                                                                                                                                tcctcTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATa                      >  1:1253433/1‑149 (MQ=255)                                                                                                                                tcctcTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATa                      >  2:1285825/1‑149 (MQ=255)                                                                                                                                ccctcTGTTTGCCATTAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTTGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATa                      <  1:647266/148‑1 (MQ=255)                                                                                                                                  ctcgGTTTGCGATGAAG‑GGAAAAAAAGAGAATGCCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACg                    <  2:260620/149‑1 (MQ=255)                                                                                                                                   tctGTTTGCGATGCAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTAACTGAGAATACGc                   <  1:1041581/149‑1 (MQ=255)                                                                                                                                   tctGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGc                   >  2:420797/1‑149 (MQ=255)                                                                                                                                   tctGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGc                   <  1:512501/149‑1 (MQ=255)                                                                                                                                      gTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTAc                <  1:420797/149‑1 (MQ=255)                                                                                                                                       tttGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGCTGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAAGTGACTGAGAGAACGCTACa               >  2:894186/1‑149 (MQ=255)                                                                                                                                          gCGATGAAG‑GAAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCt            <  2:1022092/149‑1 (MQ=255)                                                                                                                                           cGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTTTAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTc           <  1:1370974/149‑1 (MQ=255)                                                                                                                                           cGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTc           <  1:274823/149‑1 (MQ=255)                                                                                                                                           cGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTc           >  1:732461/1‑149 (MQ=255)                                                                                                                                           cGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACAAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAACTGACTGAGAATACGCTACAGcgc           >  2:219246/1‑147 (MQ=255)                                                                                                                                            gATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCg          >  1:1284874/1‑149 (MQ=255)                                                                                                                                               gAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCaa       <  1:797098/149‑1 (MQ=255)                                                                                                                                               gAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCaa       <  1:1452828/149‑1 (MQ=255)                                                                                                                                               gAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTCATTCTGGAGGAACGTTAAAGCGATTCCGACAATGGATAAATCAAAAATGACTGAGACTACCCTACAGCTCGCaa       >  2:1153571/1‑149 (MQ=255)                                                                                                                                                aaG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAAt      >  1:27520/1‑149 (MQ=255)                                                                                                                                                 aG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAGGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAAtt     <  2:1213540/149‑1 (MQ=255)                                                                                                                                                  g‑ggAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGCTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACAGAGAATACGCTACAGCTCGCAATTa    >  2:1107625/1‑149 (MQ=255)                                                                                                                                                     gAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATc  <  2:963599/149‑1 (MQ=255)                                                                                                                                                     gAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGAGAAATCAAAATTGACTGAGAATACGCTCCAGCCCGCAATTATc  >  1:786030/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAG‑GGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATC  >  NC_000913/4544441‑4544736

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