Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_pgi	3,993,512	Δ1 bp	coding (2360/2547 nt)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_pgi	3,993,512	0	A	.	87.9%	353.8 / 49.9	99	coding (2360/2547 nt)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (6/6);  new base . (43/44);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGC  >  NC_000913_3_pae_pgi/3993372‑3993653                                                                                                                                             |                                                                                                                                              gATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCg                                                                                                                                     >  2:707030/1‑150 (MQ=255)           ttCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGtt                                                                                                                           <  3:819449/150‑1 (MQ=255)           ttCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGACGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGtt                                                                                                                           >  1:1273626/1‑150 (MQ=255)               ttttCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTAc                                                                                                                       <  2:288925/150‑1 (MQ=255)                  tCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTa                                                                                                                    >  1:319231/1‑150 (MQ=255)                   cGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTAc                                                                                                                   >  1:993341/1‑150 (MQ=255)                    gaagaaACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACg                                                                                                                  >  1:61982/1‑150 (MQ=255)                    gaagaaACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACg                                                                                                                  >  1:61329/1‑150 (MQ=255)                     aagaaACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACgg                                                                                                                 >  2:1420043/1‑150 (MQ=255)                        aaaCGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTc                                                                                                              >  2:1028021/1‑150 (MQ=255)                         aaCGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTc                                                                                                              <  2:660953/149‑1 (MQ=255)                         aaCGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTc                                                                                                              >  1:660953/1‑149 (MQ=255)                           cGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAg                                                                                                           <  4:418342/150‑1 (MQ=255)                            gCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAg                                                                                                           >  1:190149/1‑150 (MQ=255)                            gCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGc                                                                                                          >  3:81282/1‑150 (MQ=255)                            gCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGc                                                                                                          >  1:321051/1‑150 (MQ=255)                             cAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCt                                                                                                         >  3:123057/1‑150 (MQ=255)                             cAAGACGAGAAGGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCt                                                                                                         <  1:1028021/150‑1 (MQ=255)                                  cGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtcatc                                                                                                    <  1:240939/150‑1 (MQ=255)                                    agaATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTGTTACCTACGGCTTAACCATCATCaa                                                                                                  >  2:1185880/1‑150 (MQ=255)                                      aaTGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGc                                                                                                          >  2:202228/1‑141 (MQ=255)                                      aaTGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGc                                                                                                          <  1:202228/141‑1 (MQ=255)                                           cTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAAc                                                                                           <  2:319231/150‑1 (MQ=255)                                            tttAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAAcc                                                                                          >  4:487144/1‑150 (MQ=255)                                            tttAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAAcc                                                                                          >  3:140732/1‑150 (MQ=255)                                                atatCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGcc                                                                                      >  3:703816/1‑150 (MQ=255)                                                 tatCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCg                                                                                     <  2:1478538/150‑1 (MQ=255)                                                 tatCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCg                                                                                     <  2:101665/150‑1 (MQ=255)                                                      aCATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGt                                                                                 <  4:736283/150‑1 (MQ=255)                                                       cATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTg                                                                                        >  1:564235/1‑141 (MQ=255)                                                       cATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTg                                                                                        <  2:564235/141‑1 (MQ=255)                                                       cATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTc                                                                               <  2:958533/150‑1 (MQ=255)                                                        aTTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCt                                                                              <  4:63013/150‑1 (MQ=255)                                                          tctcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtc                                                                                                       <  3:122816/123‑1 (MQ=255)                                                          tctcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtc                                                                                                       >  4:122816/1‑123 (MQ=255)                                                           ctcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGCTCTATc                                                                           <  4:703816/150‑1 (MQ=255)                                                            tcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtcatc                                                                                                    <  2:1428126/124‑1 (MQ=255)                                                            tcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtcatc                                                                                                    >  1:1428126/1‑124 (MQ=255)                                                            tcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCa                                                                          >  1:1227648/1‑150 (MQ=255)                                                             cgacgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAg                                                                         <  1:707030/150‑1 (MQ=255)                                                              gacgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGa                                                                        >  4:654644/1‑150 (MQ=255)                                                              gacgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGa                                                                        <  1:1336893/150‑1 (MQ=255)                                                               acgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGa                                                                        >  1:1455321/1‑150 (MQ=255)                                                               acgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCa                                                                                   >  2:551134/1‑138 (MQ=255)                                                               acgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCa                                                                                   <  1:551134/138‑1 (MQ=255)                                                                 gaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATtg                                                                     >  3:518535/1‑150 (MQ=255)                                                                 gaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATtg                                                                     >  1:605817/1‑150 (MQ=255)                                                                    aGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGa                                                                  <  2:1098904/150‑1 (MQ=255)                                                                         ccGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtcat                                                                                                     >  2:391956/1‑110 (MQ=255)                                                                         ccGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtcat                                                                                                     <  1:391956/110‑1 (MQ=255)                                                                            ggTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGc                                                                                                                                    <  2:463739/76‑1 (MQ=255)                                                                            ggTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGc                                                                                                                                    >  1:463739/1‑76 (MQ=255)                                                                            ggTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGAt                                                          >  2:174715/1‑150 (MQ=255)                                                                            ggTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGAt                                                          >  2:174748/1‑150 (MQ=255)                                                                                gAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTc                                                      <  1:174715/150‑1 (MQ=255)                                                                                gAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTc                                                      <  1:174748/150‑1 (MQ=255)                                                                                gAGGTATATCAACACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTAATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTc                                                      <  2:758016/150‑1 (MQ=255)                                                                                 aGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGc                                                                                                          <  1:585596/97‑1 (MQ=255)                                                                                 aGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGc                                                                                                          >  2:585596/1‑97 (MQ=255)                                                                                  ggTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGt                                                    <  2:1233423/150‑1 (MQ=255)                                                                                  ggTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGt                                                    <  4:136176/150‑1 (MQ=255)                                                                                   gTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCtt                                                                                                        >  3:462814/1‑97 (MQ=255)                                                                                   gTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCtt                                                                                                        <  4:462814/97‑1 (MQ=255)                                                                                     atatCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGaa                                                 <  2:61982/150‑1 (MQ=255)                                                                                     atatCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGaa                                                 >  1:545381/1‑150 (MQ=255)                                                                                     atatCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGaa                                                 <  2:61329/150‑1 (MQ=255)                                                                                       atCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACa                                               <  2:710675/150‑1 (MQ=255)                                                                                       atCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACa                                               <  1:1420043/150‑1 (MQ=255)                                                                                        tCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACa                                               >  1:1107538/1‑150 (MQ=255)                                                                                         caccacTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAg                                              <  2:190149/150‑1 (MQ=255)                                                                                           ccacTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTg                                           >  4:254620/1‑150 (MQ=255)                                                                                            cacTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGa                                          <  3:487144/150‑1 (MQ=255)                                                                                             acTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtc                                                                                                       >  1:1145552/1‑88 (MQ=255)                                                                                             acTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtc                                                                                                       <  2:1145552/88‑1 (MQ=255)                                                                                              cTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAtt                                        <  2:1376007/150‑1 (MQ=255)                                                                                               tGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAttc                                       >  1:441985/1‑150 (MQ=255)                                                                                               tGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAttc                                       >  1:442067/1‑150 (MQ=255)                                                                                                   aaGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAGCCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAttccgtt                                   <  4:318120/150‑1 (MQ=255)                                                                                                   aaGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAttccgtt                                   >  2:192584/1‑150 (MQ=255)                                                                                                       cagcaAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGc                               <  4:528024/150‑1 (MQ=255)                                                                                                                 gagCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCt                      <  2:1107538/150‑1 (MQ=255)                                                                                                                 gagCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCt                      >  2:219980/1‑150 (MQ=255)                                                                                                                 gagCTGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTa                     >  1:1248896/1‑150 (MQ=255)                                                                                                                    cTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATc                   <  2:1455321/150‑1 (MQ=255)                                                                                                                     tGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGt                                                              <  2:534140/105‑1 (MQ=255)                                                                                                                     tGGTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGt                                                              >  1:534140/1‑105 (MQ=255)                                                                                                                      ggTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCCAGCTTCATCAACTTCGACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGt                <  1:1185880/150‑1 (MQ=255)                                                                                                                       gTACGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTa               >  1:1382165/1‑150 (MQ=255)                                                                                                                          cGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtc                                                                                                       >  2:259627/1‑59 (MQ=255)                                                                                                                          cGTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTtc                                                                                                       <  1:259627/59‑1 (MQ=255)                                                                                                                           gTGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACAt           >  1:1245869/1‑150 (MQ=255)                                                                                                                            tGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTg                                                            <  1:307648/101‑1 (MQ=255)                                                                                                                            tGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTg                                                            >  2:307648/1‑101 (MQ=255)                                                                                                                            tGACGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATg          <  2:883474/150‑1 (MQ=255)                                                                                                                              aCGTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATgcc        >  1:963140/1‑150 (MQ=255)                                                                                                                                gTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATgccgc      <  2:963140/150‑1 (MQ=255)                                                                                                                                gTCAGTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATgccgc      >  3:556225/1‑150 (MQ=255)                                                                                                                                    gTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGc  <  2:1273626/150‑1 (MQ=255)                                                                                                                                    gTCGCTTCT‑CTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTATATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGc  <  2:321051/150‑1 (MQ=255)                                                                                                                                             |                                                                                                                                              GATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGC  >  NC_000913_3_pae_pgi/3993372‑3993653

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