Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,466,562	G→A	5.4%	pseudogene (1171/2513 nt)	ydbA →	pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,466,562	0	G	A	5.4%	297.5 / 5.2	94	pseudogene (1171/2513 nt)	ydbA	pseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (52/36);  new base A (3/2);  total (55/39)
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.42e-01

TCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAA  >  NC_000913/1466416‑1466709                                                                                                                                                   |                                                                                                                                                    tCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTg                                                                                                                                                  >  1:4062985/1‑150 (MQ=255)    aCGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACg                                                                                                                                               <  1:2519750/150‑1 (MQ=255)    aCGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACg                                                                                                                                               <  1:183156/150‑1 (MQ=38)    aCGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACAGTCACCGATCCGCAGTCCATCGGTATCCAGGTTGACg                                                                                                                                               >  1:1183475/1‑150 (MQ=255)     cGGTAATAACGGTAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACgg                                                                                                                                              <  2:909630/150‑1 (MQ=255)        taataaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGGTGACGGCGa                                                                                                                                           >  2:2631404/1‑150 (MQ=255)        taataaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGa                                                                                                                                           <  2:614070/150‑1 (MQ=25)           taaCGGCAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCa                                                                                                                                        <  2:1042501/150‑1 (MQ=37)           taaCGGAAAAGATATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGCTCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCa                                                                                                                                        <  2:1042467/150‑1 (MQ=255)            aaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAg                                                                                                                                       >  1:1457866/1‑150 (MQ=39)            aaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAg                                                                                                                                       >  1:3791871/1‑150 (MQ=21)            aaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAg                                                                                                                                       >  1:502025/1‑150 (MQ=21)            aaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAg                                                                                                                                       >  1:502060/1‑150 (MQ=21)               gggaaaGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCg                                                                                                                                    <  2:2950240/147‑1 (MQ=2)                   aaGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttg                                                                                                                                <  2:3330607/150‑1 (MQ=40)                    aGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgt                                                                                                                               >  2:538521/1‑150 (MQ=40)                     gtgaTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCCACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgtt                                                                                                                              <  2:426550/147‑1 (MQ=25)                     gTTATTCAGGACGGAGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGGCTGCGACCAGGCGgttgtt                                                                                                                              <  2:4113840/150‑1 (MQ=255)                         ttCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaaca                                                                                                                          >  2:1818363/1‑150 (MQ=40)                         ttCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaaca                                                                                                                          >  1:837785/1‑150 (MQ=40)                          tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCTACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTg                                                                                                                                                  <  2:2435609/125‑1 (MQ=17)                          tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaa                                                                                                                         >  2:303908/1‑150 (MQ=255)                          tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaa                                                                                                                         <  2:679709/150‑1 (MQ=255)                          tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaa                                                                                                                         >  1:3545333/1‑150 (MQ=255)                           cAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaac                                                                                                                        >  1:3427104/1‑150 (MQ=255)                             ggACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGa                                                                                                                      >  1:1051715/1‑150 (MQ=255)                              gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGaa                                                                                                                     <  1:251060/150‑1 (MQ=255)                              gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGaa                                                                                                                     <  1:715378/150‑1 (MQ=255)                                  gCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCg                                                                                                                 <  2:1051715/150‑1 (MQ=255)                                   cGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCga                                                                                                                <  2:1494675/150‑1 (MQ=255)                                     accTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaac                                                                                                                           <  2:651378/135‑1 (MQ=255)                                     aacTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaac                                                                                                                           >  1:651378/3‑137 (MQ=255)                                                aGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaac                                                                                                                        <  2:3716544/129‑1 (MQ=255)                                                aGCGGCGGCGGTCACGGTATTGATATCACCGGAGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaac                                                                                                                        >  1:3716544/1‑129 (MQ=255)                                                      ggcggTCACGGTATTGATATCACCGGGGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGt                                                                                             <  2:2018417/150‑1 (MQ=255)                                                      ggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGt                                                                                             <  2:3427104/150‑1 (MQ=255)                                                                 tATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCAc                                                                                  <  2:1457281/150‑1 (MQ=255)                                                                 tATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCAc                                                                                  <  2:342791/150‑1 (MQ=255)                                                                   ttGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGc                                                                                >  2:4188617/1‑150 (MQ=255)                                                                         tCACCGGTGACAGCGCGACGGTGGATAACAAGGGCTCCATTACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTa                                                                          <  2:571302/150‑1 (MQ=255)                                                                         tCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCa                                                                                                                                        >  1:890114/1‑88 (MQ=255)                                                                         tCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCa                                                                                                                                        <  2:890114/88‑1 (MQ=255)                                                                           aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                                                                        <  1:3746548/150‑1 (MQ=255)                                                                           aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                                                                        >  1:3515505/1‑150 (MQ=255)                                                                            ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                       >  1:2797279/1‑150 (MQ=255)                                                                            ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                       >  2:2905955/1‑150 (MQ=255)                                                                            ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                       >  2:3146744/1‑150 (MQ=255)                                                                            ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                       >  1:1275271/1‑150 (MQ=255)                                                                            ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                       >  1:3698482/1‑150 (MQ=255)                                                                             cGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACgg                                                                      >  1:1699774/1‑150 (MQ=255)                                                                              ggCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGt                                                                     <  1:749985/150‑1 (MQ=255)                                                                               gCGACAGCGCACCGGCGGATAACAAGGGCACCATGACCGTCACCGATCCGGATTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTCAACACATATGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtg                                                                    <  2:1853100/150‑1 (MQ=17)                                                                                cgtcAGCGCTACTGTTCATAAAAAGGGCACACTGACCGACTCCTCTACGGAGTTCATCGGTATCCAGGTAGACGGCGACAAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga                                                                   <  2:2205801/147‑1 (MQ=17)                                                                                cGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga                                                                   >  1:3305956/1‑150 (MQ=255)                                                                                cGACAGCCCAACTTTTGATAACAACGGCACCATGACCCTCACCTATCCGGAGTCCATATGTATCCAGGTTGACGGAGACCAGGCGTTAGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga                                                                   <  2:1390982/150‑1 (MQ=17)                                                                                    aGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGAc                                                               <  2:3305956/150‑1 (MQ=255)                                                                                     gcgcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATCaa                                                                         >  1:3888629/1‑139 (MQ=255)                                                                                     gcgcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTCGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATCaa                                                                         <  2:3888629/139‑1 (MQ=255)                                                                                      cgcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGc                                                             >  1:4041389/1‑150 (MQ=255)                                                                                             gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTCAAAAACGAAGGCGAGAGCGCCATCACCAAAGGAGGCACCGGCACTCAGATAAACGGTAAAGACGCCACCGcg                                                      >  2:1627983/1‑149 (MQ=11)                                                                                               ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCaaa                                                    <  2:1659140/150‑1 (MQ=255)                                                                                               ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCaaa                                                    <  2:837785/150‑1 (MQ=255)                                                                                               ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCaaa                                                    <  2:2597293/150‑1 (MQ=255)                                                                                                 aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACg                                                                                                                                               >  2:83227/1‑57 (MQ=35)                                                                                                 aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACg                                                                                                                                               >  2:83176/1‑57 (MQ=35)                                                                                                 aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACg                                                                                                                                               <  1:83176/57‑1 (MQ=35)                                                                                                 aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACg                                                                                                                                               <  1:83227/57‑1 (MQ=35)                                                                                                 aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaaca                                                  >  1:3663214/1‑150 (MQ=255)                                                                                                     cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaacaa                                              >  2:2349594/1‑150 (MQ=255)                                                                                                     cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaacaa                                              >  1:319816/1‑150 (MQ=255)                                                                                                         ggCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACAAGGCGTTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGACACGGCAAACAACAACGGc                                          <  2:3938296/150‑1 (MQ=255)                                                                                                            aCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGc                                                                                                                                     >  1:3013991/1‑56 (MQ=35)                                                                                                            aCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGc                                                                                                                                     >  1:3014038/1‑56 (MQ=35)                                                                                                             ccATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCCCCGGCACGCAGATTAACGGTGATGACGCCACGGCCAACAACAACGGCacaa                                      >  1:2286631/1‑150 (MQ=255)                                                                                                              cATGACCGTCACCGATCCGGAGTCCCTCGGGATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                     >  2:2265065/1‑150 (MQ=255)                                                                                                                      tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                             >  1:524205/1‑150 (MQ=255)                                                                                                                      tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                             >  1:3152364/1‑150 (MQ=255)                                                                                                                      tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGATg                             >  1:3236979/1‑150 (MQ=255)                                                                                                                      tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAAGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                             >  1:3152326/1‑150 (MQ=255)                                                                                                                      tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTAAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                             >  1:3236921/1‑150 (MQ=255)                                                                                                                        aCCGATCCGGAGTCCATCGGTATCCAGGTTTACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGAACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAAAAACGGCAAAACAACCGTTGaa                           >  1:1871422/1‑149 (MQ=255)                                                                                                                        aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAc                           >  1:838746/1‑150 (MQ=255)                                                                                                                          cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                         >  1:2181984/1‑150 (MQ=255)                                                                                                                            ctCCGGAGTTCATCTGTATCCAGTTTGACGGCGACCCGGCGTTTGTGAAGAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCa                       <  2:3697867/149‑1 (MQ=17)                                                                                                                            aTCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCa                       >  2:1922560/1‑150 (MQ=255)                                                                                                                               cGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAAAAACGAAGGCGAGAGCGCAATCAACATATGTGGCACCGTCACTAATAATAACTGATATGAAGAGACGGCAAACAAAAAAGTCAAAAAAAAAGTTGACGGGAAgg                    >  1:3159502/1‑150 (MQ=12)                                                                                                                                ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTCAAGACGCCACCGCGAAAAACAACGGAAAAACCACAGTAGACCGCAAGGa                   >  1:3435160/1‑150 (MQ=18)                                                                                                                                 gAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATAAACAACGGTGTCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAAAAACGGCAAAACCACCGTTGACGGCAAGGAt                  >  1:1947523/1‑150 (MQ=255)                                                                                                                                    tCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGTGGCACCGGCACGCAGATAAACGGTGAAGACGCCACGGCAAAAAACAAAGGCAAAACCACCGTTGAAGGCAATGATTcc               >  1:653962/1‑150 (MQ=255)                                                                                                                                       aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc            >  1:3196369/1‑150 (MQ=255)                                                                                                                                       aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc            >  2:2001181/1‑150 (MQ=255)                                                                                                                                       aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc            >  2:750894/1‑150 (MQ=255)                                                                                                                                              tCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGtggaa     <  2:1757375/150‑5 (MQ=255)                                                                                                                                                 aGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATCAACGGTCAAGACGCCACCGCGAAAAACAGAGGAAAAACCACGGTTGATGGAAAGGATTCCACGGGTAccaaaa  >  2:2113693/1‑145 (MQ=11)                                                                                                                                                   |                                                                                                                                                    TCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAA  >  NC_000913/1466416‑1466709

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